Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,649,2950TG41.2% 145.3 / 21.9 80intergenic (+254/‑314)ydfC/dicBuncharacterized protein, Qin prophage/Qin prophage; cell division inhibition protein
Reads supporting (aligned to +/- strand):  ref base T (8/39);  new base G (33/0);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 3.10e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGT  >  NC_000913/1649145‑1649444
                                                                                                                                                       |                                                                                                                                                         
gAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGt                                                                                                                                                           <  1:174238/151‑1 (MQ=255)
 aaaGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGtt                                                                                                                                                          <  1:751360/151‑1 (MQ=255)
  aaGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTa                                                                                                                                                         <  2:221349/151‑1 (MQ=255)
   aGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑t                                                                                                                                                       <  2:1454198/151‑1 (MQ=255)
   aGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑t                                                                                                                                                       <  2:832179/151‑1 (MQ=255)
   aGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑t                                                                                                                                                       <  2:1281478/151‑1 (MQ=255)
    ggATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑tt                                                                                                                                                      <  2:387529/151‑1 (MQ=255)
    ggATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑tt                                                                                                                                                      >  1:348507/1‑151 (MQ=255)
           tCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTNNNNNNNTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑ttcctt‑cc                                                                                                                                               <  1:375813/150‑1 (MQ=255)
           tCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑ccc                                                                                                                                              >  2:313264/1‑151 (MQ=255)
            cccTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑cccc                                                                                                                                             >  2:147133/1‑151 (MQ=255)
                gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGt                                                                                                                                           <  2:953933/149‑1 (MQ=255)
                gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTg                                                                                                                                         <  1:1295691/151‑1 (MQ=255)
                gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTg                                                                                                                                         <  1:82225/151‑1 (MQ=255)
                 aaaaGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCCTTgg                                                                                                                                        >  1:1182668/1‑150 (MQ=255)
                  aaaGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑TCCCtttagg                                                                                                                                       >  2:476505/1‑148 (MQ=255)
                       cAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGacac                                                                                                                                  <  1:588228/151‑1 (MQ=255)
                            cATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCCTTGGGGGCACC‑cggg                                                                                                                            >  1:436868/1‑150 (MQ=255)
                              tAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGGATGGC‑GGGGTTT‑ATCCTT‑CCCCGCTGTGGAGACC‑CGGGTg                                                                                                                          >  1:631807/1‑151 (MQ=255)
                                   cgAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGGT‑TTTCCT‑CCCCGTTTGGGGCACC‑ggggcggtagg                                                                                                                     >  2:965017/1‑143 (MQ=255)
                                    gAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTTCTT‑CCCCG‑TGAGGGGACC‑ggggcgggcggg                                                                                                                    >  2:1248774/1‑142 (MQ=255)
                                    gAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTCCTT‑CCCCCCTGAGGGGACC‑GGGGTGGCCGGt                                                                                                                    >  2:536043/1‑151 (MQ=255)
                                      aaGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTT‑CCCCCTTTAGGACAAC‑CGGGTGTCAGGGTg                                                                                                                  >  2:671264/1‑151 (MQ=255)
                                          gcATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTT‑TTTCTT‑TCCCCTTGGGGGCACC‑GGGGTGTCCGGGTGaac                                                                                                               >  1:543270/1‑148 (MQ=255)
                                          gcATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCa                                                                                                              <  1:1245227/151‑1 (MQ=255)
                                           cATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCAt                                                                                                             <  1:1404902/151‑1 (MQ=255)
                                             tGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTTTT‑CCCCCGTGAGGGGACC‑GGGGTGTCAGGTTgtgcaaa                                                                                                            >  1:854907/1‑144 (MQ=255)
                                              ggCGAGGTCGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGGGGCATCCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTACTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGAGCATACg                                                                                                          <  2:544758/151‑1 (MQ=255)
                                              ggCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATAc                                                                                                           <  1:633110/150‑1 (MQ=255)
                                                cGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCt                                                                                                        <  2:814050/151‑1 (MQ=255)
                                                cGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTT‑TCCCCCTGAGGACACC‑CGGGTGTCAGGGTGACCATACGcc                                                                                                        >  2:110313/1‑150 (MQ=255)
                                                 gAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCt‑t                                                                                                      <  1:508721/151‑1 (MQ=255)
                                                  aGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCGCTGGGGGCACC‑CGGGTGTCAGGGTGGCCATACGcctt                                                                                                      >  1:644329/1‑148 (MQ=255)
                                                   ggTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑Taa                                                                                                    <  2:1367603/151‑1 (MQ=255)
                                                    gTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑TCCCCGTGAGGGCACC‑GGGGTGTCAGGGTGACCCTACGCC‑Taa                                                                                                    >  1:216808/1‑150 (MQ=255)
                                                      aGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGt                                                                                                  <  2:734368/150‑1 (MQ=255)
                                                      aGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTg                                                                                                 <  1:89007/151‑1 (MQ=255)
                                                          ggTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACa                                                                                              <  2:678284/150‑1 (MQ=255)
                                                           gTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAAc                                                                                            <  1:210/151‑1 (MQ=255)
                                                             tCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCCTTGAGGACACC‑CGGGTGTCAGGGTGACCATACGCT‑TAAGTGACAAacc                                                                                          >  1:851635/1‑151 (MQ=255)
                                                              cATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAAccc                                                                                          <  1:1000939/150‑1 (MQ=255)
                                                              cATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTT‑CCCCCCTGGGGGCACC‑GGGGTGTCAGGGTGAACATAAGCT‑TTAGTGGCAAaccc                                                                                         >  2:783143/1‑151 (MQ=255)
                                                              cATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑TCCCCTTGAGGACAAC‑GGGGTGGCAGGGTGACCATACGCC‑TAAGTGACAAaccc                                                                                         >  1:1102026/1‑151 (MQ=255)
                                                               atagatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTT‑TCCCCCTGAGGGGACC‑CGGGTGTTTGGGTGACCCTACGCT‑TAAATGaaaaaaccc                                                                                        >  2:21164/1‑146 (MQ=255)
                                                                tagatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTCCTT‑TCCCCCTGAGGGCACC‑CGGGTGGCAGGGTGACCCTACGCT‑TTAGTGACAAaccccc                                                                                       >  2:616098/1‑149 (MQ=255)
                                                                 agatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTT‑TCCCCTTGGGGGGCCC‑GGGGTGTCAGGGTGACCCTTCCCT‑TTAGTGACAAACCCGcc                                                                                      >  1:178046/1‑150 (MQ=255)
                                                                 agatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCGGTGGGGGCACC‑GGGGTGTCAGGGTGACCCTACGCT‑TAAATGACAAacccccc                                                                                      >  1:845879/1‑148 (MQ=255)
                                                                 agatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCCTTGAGGACACC‑CGGGTGTTCGGTTGGTCATAAACT‑TAAGTGAGAAacccccc                                                                                      >  2:1465509/1‑148 (MQ=255)
                                                                 agatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCCT‑CCCCGCTGAGGACACC‑CGGGTGTTAGGGTGAACATACGCC‑TAAGGGACAACCCCCCt                                                                                      >  1:901450/1‑151 (MQ=255)
                                                                  gatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCCTTGGGGGCACC‑GGGGTGTCAGGGTGACCATACGCC‑TAAATGACAAACCCCCTg                                                                                     >  1:927698/1‑151 (MQ=255)
                                                                     agCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCGGTGAGGGCACC‑CGGGTGTCAGGTTGTACCTACGCT‑TAAATGACAACCCCCCTGcc                                                                                   >  2:486812/1‑149 (MQ=255)
                                                                      gCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAAc                                                                                 <  2:283180/151‑1 (MQ=255)
                                                                           tCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTT‑TTCTTT‑CCCCCTT‑AGGAGAAC‑CGGGTGGCAGGGTGACCCTTACCT‑TTAGTGAGCACCCCCCTCCAACGgccc                                                                            >  1:710615/1‑150 (MQ=255)
                                                                             gTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑TCCCCTTTTAGAAACC‑CGGGTGTCAGGGTGACCATACTCTCTAAAGAACAACCCC‑CTGCAACCCCctct                                                                          >  1:615918/1‑151 (MQ=255)
                                                                                aaaTTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGtt                                                                       <  1:526309/151‑1 (MQ=255)
                                                                                 aaTTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTa                                                                      <  2:348512/151‑1 (MQ=255)
                                                                                    tttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATc                                                                    <  1:913936/150‑1 (MQ=255)
                                                                                     ttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCGGTGAGGGCACC‑CGGCTGTTAGGGTGGACATACGCT‑TAAGTGACAAACCCGCTGCCACGCCCCCTGTTATCa                                                                   >  1:1274866/1‑150 (MQ=255)
                                                                                        gtcgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTTT‑CCCCCTTG‑GGGCACC‑GGGGTGTTCGGTTGGCCCTACCCC‑TTAGGGGCAACCCCCCCGCGAACCCCCCTGttttaacttt                                                               >  2:709273/1‑143 (MQ=255)
                                                                                             ccGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCGTTGAGGGCACC‑CGGGTGTCAGGTTGTCCATACCCT‑TAAATGACAAACCCCCTGCAACGCCCTCCGTTTTCAATTTTCTgg                                                          >  1:1416633/1‑151 (MQ=255)
                                                                                               gTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTg                                                        <  2:745049/151‑1 (MQ=255)
                                                                                               gTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTg                                                        <  2:529253/151‑1 (MQ=255)
                                                                                                tGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCCGTGAGGGCACC‑CGGGTGTTAGGGTGGCCATACGCC‑TTAGTGAAAAACCCCCCGCAACGCCCCCCGTTTTTCAtttttttgttg                                                       >  2:1345745/1‑144 (MQ=255)
                                                                                                 gcgcTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTT‑CCCCCCTGGGGGCACCGGGGTGGTCGGGTTGACC‑TACCCC‑TTAGTGACAAACCCCCCGCAACGCCCCCCGTTTTCAAtttttttgttac                                                      >  2:910061/1‑143 (MQ=255)
                                                                                                   gcTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGt                                                    <  2:225012/151‑1 (MQ=255)
                                                                                                     ttCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGGAAGTATGGC‑GGGGGGT‑GTCCTTCCCCCGGGGAGGA‑ACC‑CCGGTGGCAAGGTGAGCCAACCCG‑TAAATGAAAACCCCGCCGCAACGCCCCCTGCTTTTAATTTTGTTGTGACGgtt                                                  >  2:647049/1‑151 (MQ=255)
                                                                                                            tgtgGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGt                                            <  2:963160/150‑1 (MQ=255)
                                                                                                             gtgGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTa                                           <  2:227162/150‑1 (MQ=255)
                                                                                                             gtgGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTT‑TCCCCCTGTGGGCACC‑CGGGTGTCCGGTTGACCATAACCT‑TTAGTGACAAACCCCCCGCCACCCCCCCTTTTTTTAATTTTTTTGTGGCGttttggggggt                                          >  2:1392429/1‑143 (MQ=255)
                                                                                                               ggCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCa                                        <  2:1434803/151‑1 (MQ=255)
                                                                                                                 cAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGt                                      <  2:523138/151‑1 (MQ=255)
                                                                                                                   aCCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGttt                                    <  2:1088677/151‑1 (MQ=255)
                                                                                                                     ccGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTa                                  <  2:1107451/151‑1 (MQ=255)
                                                                                                                          aaaTGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTCCTT‑TCCCGCTGGGGGCACC‑GGGGTGGTAGGGTGAACATACCCC‑TAAGTGACAACCCCCCCGCAACGCCCCCCGTTTTCAATTTTTTGGGGGCGTTTTGCGGTATCAGTTTTTCTccc                             >  1:154328/1‑150 (MQ=255)
                                                                                                                                 gcggcgGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGc                      <  1:468176/151‑1 (MQ=255)
                                                                                                                                    gcggTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGctct                   <  2:300616/151‑1 (MQ=255)
                                                                                                                                               ggC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCtttt        >  1:261426/1‑151 (MQ=255)
                                                                                                                                               ggC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCtttt        >  2:866958/1‑151 (MQ=255)
                                                                                                                                               ggC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCtttt        <  1:1020517/151‑1 (MQ=255)
                                                                                                                                                 c‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTa      <  2:221634/151‑1 (MQ=255)
                                                                                                                                                   ggggTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTa      >  1:893344/1‑150 (MQ=255)
                                                                                                                                                      gTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GTGTTTTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGGGACTGCTCTTCCGCCCTTTTTAAAGt  >  2:1057954/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                         
GAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCT‑TAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGT  >  NC_000913/1649145‑1649444

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑