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

GAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGT  >  NC_000913/1649145‑1649444
                                                                                                                                                       |                                                                                                                                                         
gAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGt                                                                                                                                                           <  1:174238/151‑1 (MQ=255)
 aaaGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGtt                                                                                                                                                          <  1:751360/151‑1 (MQ=255)
  aaGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTa                                                                                                                                                         <  2:221349/151‑1 (MQ=255)
   aGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTAt                                                                                                                                                       <  2:1454198/151‑1 (MQ=255)
   aGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTAt                                                                                                                                                       <  2:832179/151‑1 (MQ=255)
   aGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTAt                                                                                                                                                       <  2:1281478/151‑1 (MQ=255)
    ggATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTAtt                                                                                                                                                      <  2:387529/151‑1 (MQ=255)
    ggATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTAtt                                                                                                                                                      >  1:348507/1‑151 (MQ=255)
           tCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTNNNNNNNTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTAttccttcc                                                                                                                                               <  1:375813/150‑1 (MQ=255)
           tCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTccc                                                                                                                                              >  2:313264/1‑151 (MQ=255)
            cccTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTcccc                                                                                                                                             >  2:147133/1‑151 (MQ=255)
                gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGt                                                                                                                                           <  2:953933/149‑1 (MQ=255)
                gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTg                                                                                                                                         <  1:1295691/151‑1 (MQ=255)
                gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTg                                                                                                                                         <  1:82225/151‑1 (MQ=255)
                 aaaaGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCTTgg                                                                                                                                        >  1:1182668/1‑150 (MQ=255)
                  aaaGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCtttagg                                                                                                                                       >  2:476505/1‑148 (MQ=255)
                       cAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGacac                                                                                                                                  <  1:588228/151‑1 (MQ=255)
                            cATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCTTGGGGGCACCcggg                                                                                                                            >  1:436868/1‑150 (MQ=255)
                              tAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGGATGGCGGGGTTTATCCTTCCCCGCTGTGGAGACCCGGGTg                                                                                                                          >  1:631807/1‑151 (MQ=255)
                                   cgAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGGTTTTCCTCCCCGTTTGGGGCACCggggcggtagg                                                                                                                     >  2:965017/1‑143 (MQ=255)
                                    gAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTTCTTCCCCGTGAGGGGACCggggcgggcggg                                                                                                                    >  2:1248774/1‑142 (MQ=255)
                                    gAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTCCTTCCCCCCTGAGGGGACCGGGGTGGCCGGt                                                                                                                    >  2:536043/1‑151 (MQ=255)
                                      aaGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTCCCCCTTTAGGACAACCGGGTGTCAGGGTg                                                                                                                  >  2:671264/1‑151 (MQ=255)
                                          gcATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTTTTTCTTTCCCCTTGGGGGCACCGGGGTGTCCGGGTGaac                                                                                                               >  1:543270/1‑148 (MQ=255)
                                          gcATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCa                                                                                                              <  1:1245227/151‑1 (MQ=255)
                                           cATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCAt                                                                                                             <  1:1404902/151‑1 (MQ=255)
                                             tGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTTTTCCCCCGTGAGGGGACCGGGGTGTCAGGTTgtgcaaa                                                                                                            >  1:854907/1‑144 (MQ=255)
                                              ggCGAGGTCGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGGGGCATCCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTACTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGAGCATACg                                                                                                          <  2:544758/151‑1 (MQ=255)
                                              ggCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATAc                                                                                                           <  1:633110/150‑1 (MQ=255)
                                                cGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCt                                                                                                        <  2:814050/151‑1 (MQ=255)
                                                cGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTTCCCCCTGAGGACACCCGGGTGTCAGGGTGACCATACGcc                                                                                                        >  2:110313/1‑150 (MQ=255)
                                                 gAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCtt                                                                                                      <  1:508721/151‑1 (MQ=255)
                                                  aGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCGCTGGGGGCACCCGGGTGTCAGGGTGGCCATACGcctt                                                                                                      >  1:644329/1‑148 (MQ=255)
                                                   ggTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTaa                                                                                                    <  2:1367603/151‑1 (MQ=255)
                                                    gTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCCGTGAGGGCACCGGGGTGTCAGGGTGACCCTACGCCTaa                                                                                                    >  1:216808/1‑150 (MQ=255)
                                                      aGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGt                                                                                                  <  2:734368/150‑1 (MQ=255)
                                                      aGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTg                                                                                                 <  1:89007/151‑1 (MQ=255)
                                                          ggTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACa                                                                                              <  2:678284/150‑1 (MQ=255)
                                                           gTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAAc                                                                                            <  1:210/151‑1 (MQ=255)
                                                             tCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCTTGAGGACACCCGGGTGTCAGGGTGACCATACGCTTAAGTGACAAacc                                                                                          >  1:851635/1‑151 (MQ=255)
                                                              cATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAAccc                                                                                          <  1:1000939/150‑1 (MQ=255)
                                                              cATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTCCCCCCTGGGGGCACCGGGGTGTCAGGGTGAACATAAGCTTTAGTGGCAAaccc                                                                                         >  2:783143/1‑151 (MQ=255)
                                                              cATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCCTTGAGGACAACGGGGTGGCAGGGTGACCATACGCCTAAGTGACAAaccc                                                                                         >  1:1102026/1‑151 (MQ=255)
                                                               atagatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTTCCCCCTGAGGGGACCCGGGTGTTTGGGTGACCCTACGCTTAAATGaaaaaaccc                                                                                        >  2:21164/1‑146 (MQ=255)
                                                                tagatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTCCTTTCCCCCTGAGGGCACCCGGGTGGCAGGGTGACCCTACGCTTTAGTGACAAaccccc                                                                                       >  2:616098/1‑149 (MQ=255)
                                                                 agatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTTCCCCTTGGGGGGCCCGGGGTGTCAGGGTGACCCTTCCCTTTAGTGACAAACCCGcc                                                                                      >  1:178046/1‑150 (MQ=255)
                                                                 agatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCGGTGGGGGCACCGGGGTGTCAGGGTGACCCTACGCTTAAATGACAAacccccc                                                                                      >  1:845879/1‑148 (MQ=255)
                                                                 agatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCTTGAGGACACCCGGGTGTTCGGTTGGTCATAAACTTAAGTGAGAAacccccc                                                                                      >  2:1465509/1‑148 (MQ=255)
                                                                 agatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCCTCCCCGCTGAGGACACCCGGGTGTTAGGGTGAACATACGCCTAAGGGACAACCCCCCt                                                                                      >  1:901450/1‑151 (MQ=255)
                                                                  gatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCTTGGGGGCACCGGGGTGTCAGGGTGACCATACGCCTAAATGACAAACCCCCTg                                                                                     >  1:927698/1‑151 (MQ=255)
                                                                     agCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCGGTGAGGGCACCCGGGTGTCAGGTTGTACCTACGCTTAAATGACAACCCCCCTGcc                                                                                   >  2:486812/1‑149 (MQ=255)
                                                                      gCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAAc                                                                                 <  2:283180/151‑1 (MQ=255)
                                                                           tCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTTCCCCCTTAGGAGAACCGGGTGGCAGGGTGACCCTTACCTTTAGTGAGCACCCCCCTCCAACGgccc                                                                            >  1:710615/1‑150 (MQ=255)
                                                                             gTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCCTTTTAGAAACCCGGGTGTCAGGGTGACCATACTCTCTAAAGAACAACCCCCTGCAACCCCctct                                                                          >  1:615918/1‑151 (MQ=255)
                                                                                aaaTTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGtt                                                                       <  1:526309/151‑1 (MQ=255)
                                                                                 aaTTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTa                                                                      <  2:348512/151‑1 (MQ=255)
                                                                                    tttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATc                                                                    <  1:913936/150‑1 (MQ=255)
                                                                                     ttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCGGTGAGGGCACCCGGCTGTTAGGGTGGACATACGCTTAAGTGACAAACCCGCTGCCACGCCCCCTGTTATCa                                                                   >  1:1274866/1‑150 (MQ=255)
                                                                                        gtcgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCCCTTGGGGCACCGGGGTGTTCGGTTGGCCCTACCCCTTAGGGGCAACCCCCCCGCGAACCCCCCTGttttaacttt                                                               >  2:709273/1‑143 (MQ=255)
                                                                                             ccGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCGTTGAGGGCACCCGGGTGTCAGGTTGTCCATACCCTTAAATGACAAACCCCCTGCAACGCCCTCCGTTTTCAATTTTCTgg                                                          >  1:1416633/1‑151 (MQ=255)
                                                                                               gTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTg                                                        <  2:745049/151‑1 (MQ=255)
                                                                                               gTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTg                                                        <  2:529253/151‑1 (MQ=255)
                                                                                                tGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCGTGAGGGCACCCGGGTGTTAGGGTGGCCATACGCCTTAGTGAAAAACCCCCCGCAACGCCCCCCGTTTTTCAtttttttgttg                                                       >  2:1345745/1‑144 (MQ=255)
                                                                                                 gcgcTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCCTGGGGGCACCGGGGTGGTCGGGTTGACCTACCCCTTAGTGACAAACCCCCCGCAACGCCCCCCGTTTTCAAtttttttgttac                                                      >  2:910061/1‑143 (MQ=255)
                                                                                                   gcTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGt                                                    <  2:225012/151‑1 (MQ=255)
                                                                                                     ttCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGGAAGTATGGCGGGGGGTGTCCTTCCCCCGGGGAGGAACCCCGGTGGCAAGGTGAGCCAACCCGTAAATGAAAACCCCGCCGCAACGCCCCCTGCTTTTAATTTTGTTGTGACGgtt                                                  >  2:647049/1‑151 (MQ=255)
                                                                                                            tgtgGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGt                                            <  2:963160/150‑1 (MQ=255)
                                                                                                             gtgGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTa                                           <  2:227162/150‑1 (MQ=255)
                                                                                                             gtgGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTTCCCCCTGTGGGCACCCGGGTGTCCGGTTGACCATAACCTTTAGTGACAAACCCCCCGCCACCCCCCCTTTTTTTAATTTTTTTGTGGCGttttggggggt                                          >  2:1392429/1‑143 (MQ=255)
                                                                                                               ggCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCa                                        <  2:1434803/151‑1 (MQ=255)
                                                                                                                 cAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGt                                      <  2:523138/151‑1 (MQ=255)
                                                                                                                   aCCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGttt                                    <  2:1088677/151‑1 (MQ=255)
                                                                                                                     ccGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTa                                  <  2:1107451/151‑1 (MQ=255)
                                                                                                                          aaaTGGCGCGGCGGTAAGTATGGCGGGGGTATTCCTTTCCCGCTGGGGGCACCGGGGTGGTAGGGTGAACATACCCCTAAGTGACAACCCCCCCGCAACGCCCCCCGTTTTCAATTTTTTGGGGGCGTTTTGCGGTATCAGTTTTTCTccc                             >  1:154328/1‑150 (MQ=255)
                                                                                                                                 gcggcgGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGc                      <  1:468176/151‑1 (MQ=255)
                                                                                                                                    gcggTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGctct                   <  2:300616/151‑1 (MQ=255)
                                                                                                                                               ggCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCtttt        >  1:261426/1‑151 (MQ=255)
                                                                                                                                               ggCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCtttt        >  2:866958/1‑151 (MQ=255)
                                                                                                                                               ggCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCtttt        <  1:1020517/151‑1 (MQ=255)
                                                                                                                                                 cGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTa      <  2:221634/151‑1 (MQ=255)
                                                                                                                                                   ggggTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTa      >  1:893344/1‑150 (MQ=255)
                                                                                                                                                      gTTATTCCTTCCCCGTTGAGGACACCGTGTTTTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGGGACTGCTCTTCCGCCCTTTTTAAAGt  >  2:1057954/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                         
GAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGT  >  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: