Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA37.0% 109.5 / 26.8 76E384D (GAA→GAT) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (10/36);  new base A (26/1);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.32e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508516‑3508805
                                                                                                                                                  |                                                                                                                                                    
gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAttta                                                                                                                                                   >  2:984930/1‑148 (MQ=255)
  tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcaca                                                                                                                                                 <  2:695365/150‑1 (MQ=255)
    tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTca                                                                                                                                                   >  2:1387661/1‑146 (MQ=255)
        tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTACTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGCCGGGGGGAAAAAATCaca‑ca‑cag                                                                                                                                          >  2:401354/1‑147 (MQ=255)
         aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTaca‑ca‑cag                                                                                                                                          >  2:795446/1‑146 (MQ=255)
         aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAg                                                                                                                                       <  1:428212/151‑1 (MQ=255)
             ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACA‑CA‑CGGGAG‑cccc                                                                                                                                  >  2:943343/1‑151 (MQ=255)
             ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑cccc                                                                                                                                  <  2:1305176/151‑1 (MQ=255)
              cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑ccccc                                                                                                                                 <  2:751629/151‑1 (MQ=255)
                gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCcac                                                                                                                               <  2:927470/151‑1 (MQ=255)
                   agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAA‑AA‑CGGAAG‑cccc                                                                                                                                  >  1:1042511/1‑145 (MQ=255)
                    gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACA‑CA‑CGGGAG‑ccccccc                                                                                                                               >  2:479439/1‑145 (MQ=255)
                    gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCAAA‑CA‑CGGGAG‑ccccc                                                                                                                                 >  1:672805/1‑145 (MQ=255)
                    gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTAAA‑CA‑CCGGAG‑gcccccc                                                                                                                               >  2:1205601/1‑145 (MQ=255)
                       cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGTGCGGCCGGGGAAATTCACA‑CA‑CGGGTA‑CCCCCCCATc                                                                                                                            >  1:1000968/1‑147 (MQ=255)
                           atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCg                                                                                                                    <  1:1154445/151‑1 (MQ=255)
                           atGCAGATACTCATGGGCGAGGGTGTAACCCGTCGTACCAAAACTCATTAGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAATAGCGACCATTATGTGGCGGTGGCGCGCGGAAAAATTGACA‑CA‑CGGTAC‑CCCCCCCATCAGAAAGc                                                                                                                     <  1:651547/150‑1 (MQ=255)
                              cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACA‑CA‑CGGGAA‑CCCCCCCACCAGAAAGCgggg                                                                                                                 >  2:976031/1‑151 (MQ=255)
                               aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGCGCGGGAAAAATCACA‑CA‑AGGGAG‑CCCCCACATCAGAAAACgagg                                                                                                                 >  1:960423/1‑150 (MQ=255)
                                gATGCTCACGGGCGAGGGTGTAACCCGTCGGATGAAAACTCATAGGGTGTCTCCCGCAACCCGGAGGGAGGAAAGAAAAGATCGCCCATGTTGTGGCGGTGGCGCGCGGAAAAAGTCACA‑CA‑CCGTAG‑CCACCACATCAGAAAGCgtggtg                                                                                                               <  1:863923/151‑1 (MQ=255)
                                   gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTgcg                                                                                                             <  2:383778/150‑1 (MQ=255)
                                   gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTgcg                                                                                                             <  1:128152/150‑1 (MQ=255)
                                      cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGcaca                                                                                                         <  2:342696/151‑1 (MQ=255)
                                      cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTgcg                                                                                                             <  1:1156695/147‑1 (MQ=255)
                                      cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCACA‑CA‑CGGTCG‑GCCCCCCCTCCGAAAGCGCGGGgcg                                                                                                             >  2:807942/1‑147 (MQ=255)
                                       aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCGCGGGAAAAATCACA‑CA‑CGGGAG‑ACCCCCCACCCGAAAAAGTGGTGGGGaca‑c                                                                                                       >  1:1313209/1‑151 (MQ=255)
                                       aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGCCGCCGGAAAAATCATC‑AA‑CAGTGG‑CCCCCCCACTAGAAAGcagggcgggggcac                                                                                                        >  2:1227994/1‑138 (MQ=255)
                                           gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGGGCGCGGGAAAAATCACA‑CA‑CGGCAG‑CCCCCCCAACAGAAAAGGGGGTGCGGACA‑cc                                                                                                      >  2:1129894/1‑148 (MQ=255)
                                             gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGGGCGGGGGAAAAATTACA‑CA‑CGGGAG‑CCCCCCCCACCGAAAGCGTGGGGCGCACA‑cc                                                                                                      >  2:641379/1‑146 (MQ=255)
                                             gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCTCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGc                                                                                                  <  1:538872/150‑1 (MQ=255)
                                             gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACA‑CA‑CCGGAG‑CCCCCACAACAAAAAGCGTGGGGCGCACA‑cc                                                                                                      >  2:1108095/1‑146 (MQ=255)
                                                     aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGc                                                                                          <  2:499868/150‑1 (MQ=255)
                                                          gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATAACA‑CA‑CGGGTG‑CCCCCCCATGCGAAAACGAGGGGCGGCCC‑CCGAGcccc                                                                                               >  1:35145/1‑137 (MQ=255)
                                                             ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATccc                                                                                 <  2:228893/151‑1 (MQ=255)
                                                               aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCtt                                                                               <  1:41077/151‑1 (MQ=255)
                                                                 cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTg                                                                              <  2:1241206/150‑1 (MQ=255)
                                                                 cAAAACTCATAGGGTGTCGCCTACAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTGGTGGCGGTGGCGAGAGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTTCGTAGCGATCATCCCTTg                                                                              <  1:947913/150‑1 (MQ=255)
                                                                       tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACC‑CA‑CGGTAG‑CCCCCCCCTCAGAAAGCGTGGGGGGCACA‑CCGAGCCGCGTAGCGGCCAACCCCTTCTGATa                                                                        >  1:1418722/1‑150 (MQ=255)
                                                                           aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                    <  2:1328683/150‑1 (MQ=255)
                                                                            gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGTCCATTTTGTGGCGGTGGCGCTCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTTCGCACA‑CCGAGCTGCTTAGCGATCATCCCTTACTGATAAACCa                                                                   <  1:1320691/150‑1 (MQ=255)
                                                                            gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTCACA‑CA‑CGGTAG‑CCCCCCCCTCAGAAAGGGTGGGGGGGACA‑CCGCGCGGCGGAGCGATCATACCTTGCTGTTAAACCa                                                                   >  1:1044039/1‑150 (MQ=255)
                                                                              gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCa                                                               <  2:1374369/151‑1 (MQ=255)
                                                                                gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACA‑CA‑CGGTAG‑CCCCCCCCCCAAAAAACGTGGTGCGCACA‑CCCAGCGGCGGAGCGATAATCCCTTGCTGTTAAAaca‑aca                                                               >  2:53768/1‑149 (MQ=255)
                                                                                 tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGGGCGCACA‑CCCAGCTGCGGAGGGATCATCCCCTCCTGGTAAAcc                                                                    >  2:144651/1‑144 (MQ=255)
                                                                                    ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGca                                                         <  1:455213/151‑1 (MQ=255)
                                                                                    ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGca                                                         <  1:512781/151‑1 (MQ=255)
                                                                                       gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTAGAGCGCGGGAAAACTCAAACCA‑CGG‑AA‑ACCCCACACCCGCAAGCCGGGGGCGCACACCCGGACTG‑GTAGATGTCATCCCCTGTTGATAAACCA‑GCACTGGcgcc                                                       >  2:675182/1‑150 (MQ=255)
                                                                                        cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCCGGGGAAAATTAACA‑CA‑CAGGAG‑CCCCCCCCACAGAAAAAGCGGGGGGCCCA‑CCCAGCCGCGGAGCGGTCCTCCCCTCCTGATTAAACA‑GCCCCGGcgccac                                                     >  1:1000272/1‑151 (MQ=255)
                                                                                        cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACA‑AA‑CAGGAG‑GCCCCACATCAGAAAACGGGGTGGCCACA‑CCGAGCCGCGGATCGATCATCCCTTGTTTATAAAACA‑ACACAGGcacca                                                      >  2:247018/1‑150 (MQ=255)
                                                                                        cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGcaccac                                                     <  2:121986/151‑1 (MQ=255)
                                                                                         aaaCCGGAGGGAGGAAAGAAAAGAGCGCCCATTTTGTGGCGGTTGCGCGCGGAAAAATGCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCCCA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGcaccac                                                     <  1:415446/150‑1 (MQ=255)
                                                                                             cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAAATCACA‑CA‑AGGACG‑ACCCCCAACCCCAAAAAGTGGTGAGCGCA‑CCCAGCTGCGTAGCGAGAATCCCCTTCTGCTAAACAA‑GCACCGGcccccc                                                     >  1:1307374/1‑144 (MQ=255)
                                                                                               gagggaggTAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCTCA‑CC‑CGGTAG‑CACCCGCATCAGAAAGCGGGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTcgcg                                               <  1:1187970/150‑1 (MQ=255)
                                                                                               gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGGGGGAAAAATTACA‑AA‑CGGGAA‑CCCCCACCTCAAAAAGCGGGGGGGGCACA‑CCCAGCTTCCGAGCGATCATCCCTTGCTTTTAAAACA‑ACACCTGGACCACCTCGCgg                                              >  2:1123416/1‑151 (MQ=255)
                                                                                                agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCCCGGAAAAAATCACA‑CA‑CCGTAA‑CCCCCCCCTCCGAAAAAGTGGGGGGCCCA‑CCCGACCGCGTCGCCCACATCCCCTGCTGATAAACCACCCACGGGC‑CCACCTCGCggg                                             >  2:1412252/1‑151 (MQ=255)
                                                                                                agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAACA‑CACCGGAAG‑CCCCC‑CCTTAAAAAACGTGGGGCGGACA‑CCCAAGAGCGGAGCGATTCTCCCCTTTTGTTAAAACA‑ACACTGGCGCCACATCCCggg                                             >  1:287253/1‑151 (MQ=255)
                                                                                                 gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGt                                            <  1:221843/151‑1 (MQ=255)
                                                                                                  ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACA‑CA‑CCGGAG‑GCCCCCCATCCGAAAACGGGGGGGGGACA‑CCGCGGCGCCGAGGGATCATCCCTTGCTGATAAAACA‑ACACCGGGCCCACCTTGCgggg                                            >  2:245621/1‑149 (MQ=255)
                                                                                                  ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTa                                           <  2:532652/151‑1 (MQ=255)
                                                                                                       aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTAAA‑CA‑CCGTAG‑CCCCCCCATCAAAAAGAGGGGGGCGCACA‑CCCCGCTGCGGAGCGATCATCCCCTGCTGATAAAACA‑ACACTGGCACCACTTCCCGGGGAttg                                        >  2:809944/1‑147 (MQ=255)
                                                                                                       aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAA‑CA‑CAGGAG‑CCCCCACCACCGAAAGCGGGGGGGCCACA‑CCCAGCTGCGTTGCGGTCATTCCCTTCTGGTAAACCA‑ACACCGGGACCCCATCGCGGGGAtgg                                        >  1:1149192/1‑149 (MQ=255)
                                                                                                          gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCt                                   <  2:1391757/151‑1 (MQ=255)
                                                                                                             aagaGCGACCATTTTGTGGCGGGGGGGCCCGGAGAAAATCACA‑CA‑CCGGGG‑ACCCCCCCACAGAAAGCCTGGTGGGGGCC‑CCCAGCTGCGGAGGGATCCACCCCTGCTGCAAAAACA‑ACACTGGCACCACCTAGCCGGGATGGTGGAtggg                                >  1:1300690/1‑149 (MQ=255)
                                                                                                              agagCGACCATTTTGTGGCGGGGTCGCCCGGAAAAAATCACA‑CA‑CCGGAGCCCCCCACA‑CTGCAAACGTGGGGGGGACA‑ACGAGACGCGGAGCGATGATCCCCCGCTGGTAAAACA‑GCACCGGCACCCCCTCCCgggggtggggggtg                                  >  1:417267/1‑137 (MQ=255)
                                                                                                                agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                             <  2:222994/151‑1 (MQ=255)
                                                                                                                 gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                            <  1:210384/151‑1 (MQ=255)
                                                                                                                     ccAGTTTGAGGCGGTGGCGCGCGGAAATACTCCCA‑CA‑CGGTAG‑CCCGCACATCAGAAAGAGGGGTGCGCACA‑CCGCGCTGTGTAGCGTTCATTCCTCGCTGATAAACCA‑GCGCTGGCACCAGGTCGCGGGTATGGTGGCTGTGACCAATgg                        <  1:775435/151‑1 (MQ=255)
                                                                                                                              ggccgtggcgCGCGGAAAAATTCACG‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑ACGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTTGCTGTGACCAATGGTCGGATCGt               <  1:1255091/151‑1 (MQ=255)
                                                                                                                               gcggtggcgCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt              <  1:902166/151‑1 (MQ=255)
                                                                                                                               gcggtggcgCCCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCTTGGTGCGCACA‑CCGAGCTGCGTAGCCCTCATCCCTTGCTGATAATCCA‑GCACTGTCACCTCTTCGCGGGTATGGTGGCTGTGCCCCATGTGCGGATCGtt              <  1:1203010/151‑1 (MQ=255)
                                                                                                                                cggtggcgCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg             <  2:541604/151‑1 (MQ=255)
                                                                                                                                 ggtggcgCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg             <  1:1018014/150‑1 (MQ=255)
                                                                                                                                       gcgcgGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt      >  1:249509/1‑151 (MQ=255)
                                                                                                                                        cgcgGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc     <  1:859330/151‑1 (MQ=255)
                                                                                                                                          cgGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg   <  2:696686/151‑1 (MQ=255)
                                                                                                                                           ggAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc  >  2:563581/1‑151 (MQ=255)
                                                                                                                                              aaaaTTCACT‑CA‑CGGTAG‑GCCCCAAATCAGAGGGCGTGGTGCGAACA‑CCGAGCTGCGTAGCGTTCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGGGGCTGGGACAAATGGTCGGATCGTGGCCGTGATCCGc  >  1:915655/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACA‑CA‑CGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACA‑CCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCA‑GCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508516‑3508805

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