Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6570AG38.2% 118.5 / 14.8 76F386S (TTC→TCC) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (8/39);  new base G (29/0);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 4.68e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508510‑3508805
                                                                                                                                                   |                                                                                                                                                         
gTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGaaaa                                                                                                                                                        <  1:1004399/151‑1 (MQ=255)
gTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGaaaa                                                                                                                                                        <  2:1141703/151‑1 (MQ=255)
gTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGaaaa                                                                                                                                                        <  2:229552/151‑1 (MQ=255)
      gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAttta                                                                                                                                                   >  2:984930/1‑148 (MQ=255)
        tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcaca                                                                                                                                                 <  2:695365/150‑1 (MQ=255)
          tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTca                                                                                                                                                   >  2:1387661/1‑146 (MQ=255)
              tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTACTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGCCGGGGGGAAAAAATCacacacag                                                                                                                                          >  2:401354/1‑147 (MQ=255)
               aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTacacacag                                                                                                                                          >  2:795446/1‑146 (MQ=255)
               aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAg                                                                                                                                       <  1:428212/151‑1 (MQ=255)
                   ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGcccc                                                                                                                                  >  2:943343/1‑151 (MQ=255)
                   ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGcccc                                                                                                                                  <  2:1305176/151‑1 (MQ=255)
                    cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGccccc                                                                                                                                 <  2:751629/151‑1 (MQ=255)
                      gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCcac                                                                                                                               <  2:927470/151‑1 (MQ=255)
                         agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAAAACGGAAGcccc                                                                                                                                  >  1:1042511/1‑145 (MQ=255)
                          gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACACACGGGAGccccccc                                                                                                                               >  2:479439/1‑145 (MQ=255)
                          gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCAAACACGGGAGccccc                                                                                                                                 >  1:672805/1‑145 (MQ=255)
                          gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTAAACACCGGAGgcccccc                                                                                                                               >  2:1205601/1‑145 (MQ=255)
                             cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGTGCGGCCGGGGAAATTCACACACGGGTACCCCCCCATc                                                                                                                            >  1:1000968/1‑147 (MQ=255)
                                 atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCg                                                                                                                    <  1:1154445/151‑1 (MQ=255)
                                 atGCAGATACTCATGGGCGAGGGTGTAACCCGTCGTACCAAAACTCATTAGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAATAGCGACCATTATGTGGCGGTGGCGCGCGGAAAAATTGACACACGGTACCCCCCCCATCAGAAAGc                                                                                                                     <  1:651547/150‑1 (MQ=255)
                                    cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAACCCCCCCACCAGAAAGCgggg                                                                                                                 >  2:976031/1‑151 (MQ=255)
                                     aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGCGCGGGAAAAATCACACAAGGGAGCCCCCACATCAGAAAACgagg                                                                                                                 >  1:960423/1‑150 (MQ=255)
                                      gATGCTCACGGGCGAGGGTGTAACCCGTCGGATGAAAACTCATAGGGTGTCTCCCGCAACCCGGAGGGAGGAAAGAAAAGATCGCCCATGTTGTGGCGGTGGCGCGCGGAAAAAGTCACACACCGTAGCCACCACATCAGAAAGCgtggtg                                                                                                               <  1:863923/151‑1 (MQ=255)
                                         gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                             <  2:383778/150‑1 (MQ=255)
                                         gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                             <  1:128152/150‑1 (MQ=255)
                                            cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGcaca                                                                                                         <  2:342696/151‑1 (MQ=255)
                                            cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                             <  1:1156695/147‑1 (MQ=255)
                                            cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCACACACGGTCGGCCCCCCCTCCGAAAGCGCGGGgcg                                                                                                             >  2:807942/1‑147 (MQ=255)
                                             aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCGCGGGAAAAATCACACACGGGAGACCCCCCACCCGAAAAAGTGGTGGGGacac                                                                                                       >  1:1313209/1‑151 (MQ=255)
                                             aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGCCGCCGGAAAAATCATCAACAGTGGCCCCCCCACTAGAAAGcagggcgggggcac                                                                                                        >  2:1227994/1‑138 (MQ=255)
                                                 gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGGGCGCGGGAAAAATCACACACGGCAGCCCCCCCAACAGAAAAGGGGGTGCGGACAcc                                                                                                      >  2:1129894/1‑148 (MQ=255)
                                                   gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGGGCGGGGGAAAAATTACACACGGGAGCCCCCCCCACCGAAAGCGTGGGGCGCACAcc                                                                                                      >  2:641379/1‑146 (MQ=255)
                                                   gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCTCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                                                                  <  1:538872/150‑1 (MQ=255)
                                                   gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCACAACAAAAAGCGTGGGGCGCACAcc                                                                                                      >  2:1108095/1‑146 (MQ=255)
                                                           aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                                                                          <  2:499868/150‑1 (MQ=255)
                                                                gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATAACACACGGGTGCCCCCCCATGCGAAAACGAGGGGCGGCCCCCGAGcccc                                                                                               >  1:35145/1‑137 (MQ=255)
                                                                   ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                                                                 <  2:228893/151‑1 (MQ=255)
                                                                     aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                               <  1:41077/151‑1 (MQ=255)
                                                                       cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                                                              <  2:1241206/150‑1 (MQ=255)
                                                                       cAAAACTCATAGGGTGTCGCCTACAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTGGTGGCGGTGGCGAGAGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTTCGTAGCGATCATCCCTTg                                                                              <  1:947913/150‑1 (MQ=255)
                                                                             tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACCCACGGTAGCCCCCCCCTCAGAAAGCGTGGGGGGCACACCGAGCCGCGTAGCGGCCAACCCCTTCTGATa                                                                        >  1:1418722/1‑150 (MQ=255)
                                                                                 aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                    <  2:1328683/150‑1 (MQ=255)
                                                                                  gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGTCCATTTTGTGGCGGTGGCGCTCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTTCGCACACCGAGCTGCTTAGCGATCATCCCTTACTGATAAACCa                                                                   <  1:1320691/150‑1 (MQ=255)
                                                                                  gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTCACACACGGTAGCCCCCCCCTCAGAAAGGGTGGGGGGGACACCGCGCGGCGGAGCGATCATACCTTGCTGTTAAACCa                                                                   >  1:1044039/1‑150 (MQ=255)
                                                                                    gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                                                               <  2:1374369/151‑1 (MQ=255)
                                                                                      gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACGGTAGCCCCCCCCCCAAAAAACGTGGTGCGCACACCCAGCGGCGGAGCGATAATCCCTTGCTGTTAAAacaaca                                                               >  2:53768/1‑149 (MQ=255)
                                                                                       tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGGAGGGATCATCCCCTCCTGGTAAAcc                                                                    >  2:144651/1‑144 (MQ=255)
                                                                                          ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                                                         <  1:512781/151‑1 (MQ=255)
                                                                                          ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                                                         <  1:455213/151‑1 (MQ=255)
                                                                                             gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTAGAGCGCGGGAAAACTCAAACCACGGAAACCCCACACCCGCAAGCCGGGGGCGCACACCCGGACTGGTAGATGTCATCCCCTGTTGATAAACCAGCACTGGcgcc                                                       >  2:675182/1‑150 (MQ=255)
                                                                                              cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCCGGGGAAAATTAACACACAGGAGCCCCCCCCACAGAAAAAGCGGGGGGCCCACCCAGCCGCGGAGCGGTCCTCCCCTCCTGATTAAACAGCCCCGGcgccac                                                     >  1:1000272/1‑151 (MQ=255)
                                                                                              cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACAGGAGGCCCCACATCAGAAAACGGGGTGGCCACACCGAGCCGCGGATCGATCATCCCTTGTTTATAAAACAACACAGGcacca                                                      >  2:247018/1‑150 (MQ=255)
                                                                                              cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcaccac                                                     <  2:121986/151‑1 (MQ=255)
                                                                                               aaaCCGGAGGGAGGAAAGAAAAGAGCGCCCATTTTGTGGCGGTTGCGCGCGGAAAAATGCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCCCACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcaccac                                                     <  1:415446/150‑1 (MQ=255)
                                                                                                   cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAAATCACACAAGGACGACCCCCAACCCCAAAAAGTGGTGAGCGCACCCAGCTGCGTAGCGAGAATCCCCTTCTGCTAAACAAGCACCGGcccccc                                                     >  1:1307374/1‑144 (MQ=255)
                                                                                                     gagggaggTAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCTCACCCGGTAGCACCCGCATCAGAAAGCGGGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                               <  1:1187970/150‑1 (MQ=255)
                                                                                                     gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGGGGGAAAAATTACAAACGGGAACCCCCACCTCAAAAAGCGGGGGGGGCACACCCAGCTTCCGAGCGATCATCCCTTGCTTTTAAAACAACACCTGGACCACCTCGCgg                                              >  2:1123416/1‑151 (MQ=255)
                                                                                                      agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCCCGGAAAAAATCACACACCGTAACCCCCCCCTCCGAAAAAGTGGGGGGCCCACCCGACCGCGTCGCCCACATCCCCTGCTGATAAACCACCCACGGGCCCACCTCGCggg                                             >  2:1412252/1‑151 (MQ=255)
                                                                                                      agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAACACACCGGAAGCCCCCCCTTAAAAAACGTGGGGCGGACACCCAAGAGCGGAGCGATTCTCCCCTTTTGTTAAAACAACACTGGCGCCACATCCCggg                                             >  1:287253/1‑151 (MQ=255)
                                                                                                       gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                            <  1:221843/151‑1 (MQ=255)
                                                                                                        ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACCGGAGGCCCCCCATCCGAAAACGGGGGGGGGACACCGCGGCGCCGAGGGATCATCCCTTGCTGATAAAACAACACCGGGCCCACCTTGCgggg                                            >  2:245621/1‑149 (MQ=255)
                                                                                                        ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTa                                           <  2:532652/151‑1 (MQ=255)
                                                                                                             aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTAAACACCGTAGCCCCCCCATCAAAAAGAGGGGGGCGCACACCCCGCTGCGGAGCGATCATCCCCTGCTGATAAAACAACACTGGCACCACTTCCCGGGGAttg                                        >  2:809944/1‑147 (MQ=255)
                                                                                                             aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAACACAGGAGCCCCCACCACCGAAAGCGGGGGGGCCACACCCAGCTGCGTTGCGGTCATTCCCTTCTGGTAAACCAACACCGGGACCCCATCGCGGGGAtgg                                        >  1:1149192/1‑149 (MQ=255)
                                                                                                                gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCt                                   <  2:1391757/151‑1 (MQ=255)
                                                                                                                   aagaGCGACCATTTTGTGGCGGGGGGGCCCGGAGAAAATCACACACCGGGGACCCCCCCACAGAAAGCCTGGTGGGGGCCCCCAGCTGCGGAGGGATCCACCCCTGCTGCAAAAACAACACTGGCACCACCTAGCCGGGATGGTGGAtggg                                >  1:1300690/1‑149 (MQ=255)
                                                                                                                    agagCGACCATTTTGTGGCGGGGTCGCCCGGAAAAAATCACACACCGGAGCCCCCCACACTGCAAACGTGGGGGGGACAACGAGACGCGGAGCGATGATCCCCCGCTGGTAAAACAGCACCGGCACCCCCTCCCgggggtggggggtg                                  >  1:417267/1‑137 (MQ=255)
                                                                                                                      agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                             <  2:222994/151‑1 (MQ=255)
                                                                                                                       gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                            <  1:210384/151‑1 (MQ=255)
                                                                                                                           ccAGTTTGAGGCGGTGGCGCGCGGAAATACTCCCACACGGTAGCCCGCACATCAGAAAGAGGGGTGCGCACACCGCGCTGTGTAGCGTTCATTCCTCGCTGATAAACCAGCGCTGGCACCAGGTCGCGGGTATGGTGGCTGTGACCAATgg                        <  1:775435/151‑1 (MQ=255)
                                                                                                                                    ggccgtggcgCGCGGAAAAATTCACGCACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACAACGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTTGCTGTGACCAATGGTCGGATCGt               <  1:1255091/151‑1 (MQ=255)
                                                                                                                                     gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt              <  1:902166/151‑1 (MQ=255)
                                                                                                                                     gcggtggcgCCCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCTTGGTGCGCACACCGAGCTGCGTAGCCCTCATCCCTTGCTGATAATCCAGCACTGTCACCTCTTCGCGGGTATGGTGGCTGTGCCCCATGTGCGGATCGtt              <  1:1203010/151‑1 (MQ=255)
                                                                                                                                      cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg             <  2:541604/151‑1 (MQ=255)
                                                                                                                                       ggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg             <  1:1018014/150‑1 (MQ=255)
                                                                                                                                             gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt      >  1:249509/1‑151 (MQ=255)
                                                                                                                                              cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc     <  1:859330/151‑1 (MQ=255)
                                                                                                                                                cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg   <  2:696686/151‑1 (MQ=255)
                                                                                                                                                 ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc  >  2:563581/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                         
GTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508510‑3508805

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