Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913213,3410AG43.8% 125.1 / 32.4 90G337G (GGA→GGGtilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (13/37);  new base G (39/0);  total (53/37)
Fisher's exact test for biased strand distribution p-value = 2.46e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGCC  >  NC_000913/213193‑213490
                                                                                                                                                     |                                                                                                                                                      
aTGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAg                                                                                                                                                      >  3:484184/1‑151 (MQ=255)
aTGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAg                                                                                                                                                      <  2:153110/151‑1 (MQ=255)
     tCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTAc                                                                                                                                                 <  1:268664/151‑1 (MQ=255)
     tCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTAc                                                                                                                                                 <  1:508326/151‑1 (MQ=255)
      cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTAc                                                                                                                                                 >  4:501202/1‑150 (MQ=255)
      cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTAc                                                                                                                                                 <  2:1331886/150‑1 (MQ=255)
      cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACa                                                                                                                                                <  4:788839/151‑1 (MQ=255)
           tGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGGACAGCtt                                                                                                                                            >  2:22164/1‑150 (MQ=255)
            gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTaa                                                                                                                                          >  3:308218/1‑151 (MQ=255)
            gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTa                                                                                                                                           <  2:971966/150‑1 (MQ=255)
            gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTa                                                                                                                                           <  4:867599/150‑1 (MQ=255)
              ttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAAt                                                                                                                                         <  4:861648/150‑1 (MQ=255)
               tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAACTTAATGc                                                                                                                                       >  2:878492/1‑151 (MQ=255)
               tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATg                                                                                                                                        <  2:1028748/150‑1 (MQ=255)
                  gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCggg                                                                                                                                    >  2:1113385/1‑151 (MQ=255)
                    ttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGAGGACAGCTTAATGCggggg                                                                                                                                  >  2:841739/1‑149 (MQ=255)
                     tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATACgggggg                                                                                                                                 >  2:902291/1‑151 (MQ=255)
                      gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAgg                                                                                                                                 <  2:144979/150‑1 (MQ=255)
                         cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGat                                                                                                                             <  4:1192634/151‑1 (MQ=255)
                               tGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCGGGGGAGGGGAAAGCTTAAATCGGGAGGGGGTAATCGc                                                                                                                       >  4:491811/1‑151 (MQ=255)
                                gAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGGGGTATTTGcc                                                                                                                      >  2:2830/1‑151 (MQ=255)
                                gAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGc                                                                                                                       <  2:281802/150‑1 (MQ=255)
                                 aaaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGccc                                                                                                                     <  4:93243/151‑1 (MQ=255)
                                  aaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGAGCACTTAATGAGGGAGGGGGTATTCccccc                                                                                                                    >  4:1269550/1‑150 (MQ=255)
                                     ccGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGGGGGGAGGGTGAAGCTTAAAGCGGGGGGCGATATTTGccccccc                                                                                                                 >  4:1019519/1‑150 (MQ=255)
                                         cGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGAAAGCTTTATGCGGGGGGGGGTATTCGCGCTCcccggg                                                                                                             >  4:524615/1‑149 (MQ=255)
                                           cTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTTATACGGGAGGGGGTATTCGCCCTCCGCGTGc                                                                                                            >  4:452734/1‑150 (MQ=255)
                                            tATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATTCGGGAGGCGATATTCTCCCCCCCCGTGc                                                                                                            >  4:554147/1‑149 (MQ=255)
                                             aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCa                                                                                                           <  4:967777/149‑1 (MQ=255)
                                              aaaGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGAc                                                                                                        <  2:947861/149‑1 (MQ=255)
                                                aGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGa                                                                                                      <  4:956412/151‑1 (MQ=255)
                                                 gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGaa                                                                                                     >  1:260925/1‑151 (MQ=255)
                                                 gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCTGGCGGAGGGCGGCTTCATGGGGGAGGGGATTTTTGTCCTCCGGGTGCCGAcaga                                                                                                     >  1:539874/1‑148 (MQ=255)
                                                       aCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGACAGCTTAAAGCGGGGGGGGATATTTGCCCCCCCCGGGCAGACGAAGCgggc                                                                                               >  4:144335/1‑149 (MQ=255)
                                                          gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAg                                                                                             <  4:1256139/150‑1 (MQ=255)
                                                              tggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATGCGGGGGGGGATATTTGCCCTCCGCGTGCAGACGAAGCGGGCAgcgtgg                                                                                        >  2:578864/1‑150 (MQ=255)
                                                               ggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcgtgcg                                                                                       <  4:384521/151‑1 (MQ=255)
                                                               ggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcgtgc                                                                                        <  2:1114818/150‑1 (MQ=255)
                                                                gATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcgtgcgt                                                                                      <  2:801475/151‑1 (MQ=255)
                                                                 aTTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGCGGTATTCGCCCTCCCCGTGCAGACGAAGCGGGCAgcgggcggt                                                                                     >  2:1060323/1‑149 (MQ=255)
                                                                  ttAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGTGTACAGCTTAATGCGGGGGGCGGTATTTGCCCTCCGCGGGCAGACGAAGCGGGCAGCGGGCGgtt                                                                                    >  4:1429940/1‑151 (MQ=255)
                                                                   tAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGAGGGTACAGTTTAATGCGGGTGGCGGAATTCGCCCTCCGCGCGCAGACGAAACgggtaacggggggttt                                                                                   >  4:287253/1‑137 (MQ=255)
                                                                   tAAATCCGTCACCGGGCAAAGCGAAAACATTGTGACGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAACGGAATAGTTAAATCGGGTGGGGGTATTTAGCCTCAGGCTGCAGACGAAGCgggctcggagcggttc                                                                                   >  2:374536/1‑137 (MQ=255)
                                                                    aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGCGTTTc                                                                                   >  2:1237842/1‑150 (MQ=255)
                                                                    aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGAGGCGATATTTGCCCTCCGCGTGCAGACGAAGCGGGCCGCGGGCGtttt                                                                                   >  2:801293/1‑149 (MQ=255)
                                                                    aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCa                                                                                  <  2:766877/151‑1 (MQ=255)
                                                                     aaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCaa                                                                                 <  4:420081/151‑1 (MQ=255)
                                                                      aTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATTCGGGAGGCGATATTCGCCCTCCGCGGGCAGACGAAGCGGGCCGCGGGCGTTTCaaa                                                                                >  1:568316/1‑151 (MQ=255)
                                                                      aTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCCCCTGCAGACGAAGCGGTCAGCGGGCGGTTCaa                                                                                 >  2:1311411/1‑150 (MQ=255)
                                                                           tCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGCGATATTTGCCCTCCCCGTGCAGACGAAACGGTCAGCGGGCGGTTCAAAACGcc                                                                           >  4:381454/1‑151 (MQ=255)
                                                                              ccGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTCCCCCTCCCCGTGCAGACGAAGCGGTCAGCGTGCGGTTTAAAACGCCCgg                                                                        >  2:1326153/1‑151 (MQ=255)
                                                                              ccGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAgg                                                                        <  4:1193114/151‑1 (MQ=255)
                                                                                gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGGGCGGTTCAAAACGCCAGGAt                                                                      >  2:706193/1‑151 (MQ=255)
                                                                                 ggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGAtt                                                                     <  4:319268/151‑1 (MQ=255)
                                                                                   cAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGCGGTTCAAAACGCCAGGGTtgc                                                                   >  2:656098/1‑151 (MQ=255)
                                                                                      aGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGATTAATGGGGGGGGGGGTATTTGCCCTCCGCGTGCAGACGGAACGGGCCGCGGGCGGTTTAAAACGCCCGGATtgctgg                                                                >  2:1435078/1‑150 (MQ=255)
                                                                                          aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCatat                                                            <  3:275283/151‑1 (MQ=255)
                                                                                          aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCatat                                                            <  2:580679/151‑1 (MQ=255)
                                                                                           aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATAtt                                                           <  4:194253/151‑1 (MQ=255)
                                                                                            aaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGAGGGGACAGCTTAAAGCGGGGGGCGATATTCGCCCCCCGCGTGCAGACGAAGCGGGCAGCGTGCGGTTCCAAGCGCCACGATTGCTGCCTAttt                                                          >  4:1431149/1‑150 (MQ=255)
                                                                                             aCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGTCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGAATATTGt                                                         <  3:137046/151‑1 (MQ=255)
                                                                                             aCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGGACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGGGCAGACGAAGCGGGCAGCGGGCGGTTCAAAGCGCCAGGGTTGCTGCATATTGt                                                         >  1:159636/1‑151 (MQ=255)
                                                                                                    gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGt                                                 <  2:1402605/151‑1 (MQ=255)
                                                                                                       gTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCCTAATGCGGGAGGCGATATTCCCCCCCCGCGTGCAGACGAAGCGGTCAGCGTGCGGTTCAAAACGCCAGGATTGCTGCATATTGTCGGGCGGAAc                                              >  2:1431651/1‑151 (MQ=255)
                                                                                                       gTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAAAGCGGGGGGCGATATTTGCCCTCCGCGTGCAGACGAAGCGGTCAGCGGGCGGTTCAAAGCGCCAGGATTGCTGCATATTTTCggggggaa                                               >  4:504477/1‑145 (MQ=255)
                                                                                                         ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATTCGGGAGGGGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGGGGGGTTTAAAACGCCAGGATTGCTGCATATTGTCGGGCGGAAcgg                                            >  2:1055133/1‑151 (MQ=255)
                                                                                                         ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGCGGTATTCGCCCCCCGCGTGCAGACGAAGCGGGCAGCGTGCGTTTCAAAACGCCCGGGTTGCTGCATATTTTCGGGCGGAAcgg                                            >  2:520873/1‑151 (MQ=255)
                                                                                                         ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGCGGTGGGGAGTGGACAGATTAATGAGGGGGGGGGGATTTGCCCCCCGCCGGGCGACGAAACGGGGAGCGGGGGGTTCAAAGCGCCCGGATTGGTTCCTATTGTCCGGCGTAAcgg                                            >  4:877320/1‑151 (MQ=255)
                                                                                                          gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggc                                           <  2:1366491/151‑1 (MQ=255)
                                                                                                            aGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGGGGTATTTGCCCTCCCCGTGCAGACGAAGCGGGCAGCGGGCGGTTTAAAACGCCAGGAGTGGTGCATATTGTCGGGCGGAAaggcgg                                         >  4:1301942/1‑151 (MQ=255)
                                                                                                            aGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCCTAATGCGGGGGGGGGTATTTGCCCTCCCCGTGCCGACGAAGCGGGGAGCGGGCGGTTCAAAACGCCAGGATTGCTGCATATTGTCGGGCGTAAaggcgg                                         >  4:179182/1‑151 (MQ=255)
                                                                                                               cGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGAGTAAAACCGGGGGGGATATTTGACCCCCCCGCGCCGGCGGAACGGGCAGAGCGCGTTTCCAAACGCCCGGATTGCTGCATATTGTAGGGGGTAACGGGGGGCg                                      >  4:1392904/1‑151 (MQ=255)
                                                                                                                gTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGt                                     <  4:296420/151‑1 (MQ=255)
                                                                                                                  ggCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTa                                    <  2:813179/150‑1 (MQ=255)
                                                                                                                   gCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGGGCGTTTCAAAACGCCCGGGTTGCTGCATATTTTCGGGGGGAACGGCGGGCGGAAg                                  >  3:461998/1‑151 (MQ=255)
                                                                                                                   gCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAg                                  <  2:165351/151‑1 (MQ=255)
                                                                                                                      tCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTAAAGCTTAATTCGGGGGGCGATATTTGCCCTCCCCGTGCAGACGAAGCGGTCCGCGGGCGTTTCAAAACGCCAGGATTGGTGCATATTGTCGGGCGTAACGGGGGACGTAAGCtt                               >  2:1233633/1‑150 (MQ=255)
                                                                                                                      tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCCCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGGTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGGCGTAAggt                                >  2:1402513/1‑148 (MQ=255)
                                                                                                                              tGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAAt                       <  1:263826/151‑1 (MQ=255)
                                                                                                                               ggAATTACCGGCGGGGCTGGGGAGTGGAAAGCTTAATGCGGGGGGGGATATTCGCCCTCCCCGTGCAGACGAAGCGGTCAGCGGGCGGTTCAAAACCCCCGGATTGCTGCATATTTTCGGGGGGCAACGGGGGCGTAAGGTAAAAAAAAtt                      >  4:237275/1‑150 (MQ=255)
                                                                                                                                  aTTACCGGCGGGGGTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCCCCGCGGGCAGACGAAGCGGTCAGCGTGCGGTTTAAAGCGCCAGGATTGCTGCATAATTTCGGGCGTAACGGCGGGCGGAAGCTAAAAAAAATCTgg                   >  4:1013686/1‑151 (MQ=255)
                                                                                                                                          cggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCt           >  4:276911/1‑151 (MQ=255)
                                                                                                                                           ggggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTg          <  2:491807/151‑1 (MQ=255)
                                                                                                                                               cTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCg      >  4:477660/1‑151 (MQ=255)
                                                                                                                                               cTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCg      <  3:42908/151‑1 (MQ=255)
                                                                                                                                                tGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGGTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCgg     >  2:1152100/1‑150 (MQ=255)
                                                                                                                                                    gAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgcc  <  4:1326152/151‑1 (MQ=255)
                                                                                                                                                    gAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgcc  >  4:1239778/1‑151 (MQ=255)
                                                                                                                                                    gAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgcc  >  4:1032407/1‑151 (MQ=255)
                                                                                                                                                    gAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgc   >  1:526413/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGCC  >  NC_000913/213193‑213490

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