Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA40.0% 139.7 / 40.7 90E384D (GAA→GAT) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (12/42);  new base A (36/0);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 4.10e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACG  >  NC_000913/3508518‑3508812
                                                                                                                                                 |                                                                                                                                                         
tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATtac‑a                                                                                                                                                     >  2:780530/1‑150 (MQ=255)
  tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATcac‑aca                                                                                                                                                   >  1:95155/1‑150 (MQ=255)
  tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGG‑GGAAAAAATcac                                                                                                                                                       >  2:854359/1‑147 (MQ=255)
   ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCtc‑tcac                                                                                                                                                  <  2:809922/150‑1 (MQ=255)
   ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGC‑GGGAAAAATtac‑aca                                                                                                                                                   >  1:355248/1‑149 (MQ=255)
       aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTa                                                                                                                                              <  1:1079462/150‑1 (MQ=255)
       aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTa                                                                                                                                              <  2:1256283/150‑1 (MQ=255)
       aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGC‑GGGAAAAATAAC‑ACACCGt                                                                                                                                               >  2:298727/1‑149 (MQ=255)
         aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGcc                                                                                                                                           <  1:149180/151‑1 (MQ=255)
          gCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGccc                                                                                                                                          <  1:136379/151‑1 (MQ=255)
           ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTaac‑aacaggagtccc                                                                                                                                          >  2:121251/1‑140 (MQ=255)
           ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGCGC‑GGGAAAAAACAC‑AAACAGTCGgccc                                                                                                                                         >  2:1031769/1‑151 (MQ=255)
             ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCca                                                                                                                                       <  1:448068/151‑1 (MQ=255)
             ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGC‑GGCAAAAATCAT‑ACACgggataacgc                                                                                                                                        >  1:628775/1‑142 (MQ=255)
               agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCca                                                                                                                                       <  1:820456/149‑1 (MQ=255)
                 agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAC‑ACACCGGAGccccccc                                                                                                                                      >  1:467506/1‑146 (MQ=255)
                  gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGC‑CGAAAAAATTAC‑ACCCCGGAGCCCCcac                                                                                                                                      >  1:1108984/1‑147 (MQ=255)
                     cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAC‑ACACCGGAACCCCCCCAACa                                                                                                                                  >  1:82960/1‑148 (MQ=255)
                     aaaTATGAAGATGCTCACGGGCGAGGGAGTAAGCCGTCGGATCAAAACTCATAGGCTGTCTCCTGCAAACCGGAGGGAGGAAAGACAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑CGAAAAATTtac‑acaaaataccccc                                                                                                                                         >  2:958150/2‑131 (MQ=255)
                          tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAC‑ACACGGTAGgccccccccccgaaagacg                                                                                                                          >  1:1205403/1‑137 (MQ=255)
                           gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAA‑ACACCGGAGCCCCCACATCACA‑AAACggg                                                                                                                        >  1:328438/1‑149 (MQ=255)
                           gCAGATGCTCATGGGCGAGGCTGTATACCGTCGCATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCACTTTGCGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACAGGGTAGCCCCCACATCAGA‑AAGCgg                                                                                                                         <  2:444751/150‑2 (MQ=255)
                            cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCCC‑GGGAAAATTTAC‑CCACCGGAGGCCCCCCAACAGG‑AAGCgtgg                                                                                                                       >  1:41968/1‑151 (MQ=255)
                               aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGC‑GGGAAAAATTAA‑AAACCGGAGCCCCCACACCAAA‑AAACgggg                                                                                                                       >  2:1056320/1‑148 (MQ=255)
                                tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTAAA‑AAAC‑GTAGCCCCCACCTCAGA‑AAACggggtggg                                                                                                                   >  2:198311/1‑149 (MQ=255)
                                           gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGG‑GGGAAAAATTAA‑ACACAGTAGCCCCCCCATCAGA‑AAAGGTGGTGCGCACCCCGAGc                                                                                                         >  1:632141/1‑150 (MQ=255)
                                           gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAC‑ACACGGGAGACCCCCCAACACA‑AAACGGGGGGGGcac                                                                                                                >  1:868421/1‑143 (MQ=255)
                                           gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAC‑ACACCGTAGCCCCCCCCTCAGA‑AAGCGGGGGGGGCACAcc                                                                                                             >  1:555901/1‑146 (MQ=255)
                                           gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAA‑AAACGGGAGGCCCCCCCACACA‑AAAAGTGGGGGGGACACCCAGc                                                                                                         >  2:531267/1‑150 (MQ=255)
                                            aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGAAAAAATAACAACCCGGGCG‑CCCCACAACACA‑AAGCGTGGTGTGCGCAcc                                                                                                             >  2:250340/1‑145 (MQ=255)
                                                    aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGGGG‑GGGAAAAATCAC‑ACCCGGGAGGCCCCCCATCACA‑AAAACCGGGGCGGACACCGAGCCGCCGAGCg                                                                                                >  2:381131/1‑150 (MQ=255)
                                                     cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                                                                              <  1:217856/151‑1 (MQ=255)
                                                     cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                                                                              <  1:1068594/151‑1 (MQ=255)
                                                          cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAC‑ACACCGTAGCCCCCCCAACAGA‑AAACGTGGTGGGCACACCGAGCCGCGTAGCGATCA‑Tcc                                                                                        >  1:130672/1‑151 (MQ=255)
                                                               cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGc                                                                                   <  2:144316/151‑1 (MQ=255)
                                                                aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAC‑ACACGGGAGCCCCCCCATCAGA‑AAACGTGGGGCGCACACCGAGCTGCGTAGCGATCC‑Tccc                                                                                       >  1:308019/1‑146 (MQ=255)
                                                                 aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAC‑ACACCGTAGCCCCCCCACCAGA‑AAACGGGGTGGGCACACCCAGCCGCGTAGCGAACA‑Tcccct                                                                                     >  1:840392/1‑145 (MQ=255)
                                                                 aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTg                                                                                 <  2:869596/151‑1 (MQ=255)
                                                                  aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGC‑GGGAAAAATCAC‑ACAAAGGAACCCCCCCCTCCGA‑AAACGAGGTGGGGACACCGAGCTGCGTAGCGAGCA‑ACCCTTGc                                                                                   >  2:1276045/1‑148 (MQ=255)
                                                                      cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATaa                                                                             <  2:996842/150‑1 (MQ=255)
                                                                      cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATaa                                                                             <  1:1227065/150‑1 (MQ=255)
                                                                      cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGGGC‑CGGAAAAATCAC‑ACACGGGAGGCCCCCCCCCAGA‑AAAAGGGGGGCGCACACCGAAGGGCGGAGGGGGCA‑TCCCCTGCTGGTaa                                                                             >  1:47200/1‑150 (MQ=255)
                                                                        tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGC‑GGGGAAAATCAC‑ACACCGGAGCCCCCCCATCAGA‑AAACGGGGGGGGCACACCGAGCTGCGGAGCGCTCA‑TCCCCTGCTGATAAAcc                                                                          >  1:921133/1‑151 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAC‑ACACCGGAGCCCCCCCATCAGA‑AAGCGGGGGGCGCCCACCCAGCTGCGTAGCGGACC‑TCCCCTGCTGATAAAcc                                                                          >  1:238447/1‑150 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAC‑ACACCGGAGCCCCCCCATCAGA‑AAACGGGGGGCGCACACCGAGCCGCGGAGCGCTCA‑TCCCCTTCTGATaaaa                                                                           >  1:140052/1‑148 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGC‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCa                                                                         <  2:410452/151‑1 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCa                                                                         <  1:239086/151‑1 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGGGCGG‑GGAAAAAATCAC‑ACACGGGAGCCCCCCCATACGAGAAGAG‑GGGGGGGACACCGCGCTGCGGAGGGGTCA‑TCCCGTGCTGaaaa                                                                             >  2:506447/1‑147 (MQ=255)
                                                                          gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAC‑AAAACGTAGCCCCCACAACAGA‑AAAGGGGGGGGGCCCACCCAGCTGCGGAGCGATCA‑TCCCCTGCTGATaaa                                                                            >  1:838256/1‑147 (MQ=255)
                                                                           ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAA‑AAACGGTAGCCCCCCCCTCAGA‑AAACGGGGGGGGGACACCCAGCTGCGTAGCGATCA‑TCCCCTTCTGTTAAACCa                                                                         >  1:1220896/1‑149 (MQ=255)
                                                                              gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCAc                                                                     <  1:854390/150‑1 (MQ=255)
                                                                              gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGG‑GGGAAAATTTTG‑ACAAGGGAACCCCCACCTCACA‑AAAAGCGGTGCGGAGCCCGAGCAGCGGGGCCAACA‑TCCCTTGCTTGCAAACCAGc                                                                       >  1:109397/1‑148 (MQ=255)
                                                                                ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATAAC‑ACACAGGAACCCCCACCACAGA‑AAAAGGGGGGGGCAACCCCAGCCGCGGAGCGATAA‑TCCCCTGCTTATAAAACAGCACttg                                                                  >  2:579677/1‑149 (MQ=255)
                                                                                        aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACtt                                                          <  2:578439/151‑1 (MQ=255)
                                                                                         aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                         <  1:1082601/151‑1 (MQ=255)
                                                                                            ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                                      <  1:687657/151‑1 (MQ=255)
                                                                                              agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAA‑ACACCGTAACCCCCCCATCAGA‑AAGCGGGGGGGGCACACCCAGCTGCGGAGCGATCA‑TCCCCTTCTTATAAAACAACACCGGCACCACCTCCCggg                                                    >  1:1165285/1‑151 (MQ=255)
                                                                                              agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                                    <  2:625876/151‑1 (MQ=255)
                                                                                              agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                                    <  2:26859/151‑1 (MQ=255)
                                                                                              agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGG‑GGGAAAAATAAC‑ACACCGGAACCCCCACCTCAGA‑AAACGGGGGGCGCCCACCGAGCTGCGGAGAGGTCA‑TCCCTTTCTGGTAAAAAACCACCGGCACCACTTCGCggg                                                    >  2:1034310/1‑151 (MQ=255)
                                                                                               gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                                    <  2:446997/150‑1 (MQ=255)
                                                                                               gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                                   <  1:601951/151‑1 (MQ=255)
                                                                                                    gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggt                                              <  1:530125/151‑1 (MQ=255)
                                                                                                     aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAAACAC‑ACACGGTAGCCCCCCCATCCGA‑AAACGGGGGGGGCACACCCAGCCGGGGGAGGAGCA‑TCCCCTTCTTGTAAAACAACACTGGCACCCCTTCCCGGGGAtgggg                                             >  2:247897/1‑149 (MQ=255)
                                                                                                      aagaaaagaGCGACCATTTTGTGGCGGGGGCGCGC‑GGGAAAAATAAA‑ACACCGGAGCCCCCACATCAGA‑AAGCGGGGGGCGCACCCCGAGCTGCGTAGCGATCA‑TCCCTTTCTTATAAAACAGCACCGGGACCACCTCCCGGGGGtggtgg                                            >  2:686552/1‑151 (MQ=255)
                                                                                                       agaaaagaGCGACCATTTTGTGGCGGAGGCGCGC‑GGAAAAATTAAA‑AAACGGGTACCCCCCCAACAGC‑AAGAGCGGGGGGCACACCGCGCCCCGGAGGGGTCATTCCCTTCTTGTAAAACCA‑CACCGGGACCACCTCGGGGGtgttg                                               >  2:57126/1‑143 (MQ=255)
                                                                                                        gaaaagaGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATAAC‑AAACGGCAGCCCCCCCATCAGA‑AAACGGGGGGGGCACACCCCGCCGCGGAGCGAGCA‑TCCCTTGCTGGTAAAACAACACCGGCACCCCATCCCGGGGAtggtgg                                            >  2:405569/1‑149 (MQ=255)
                                                                                                            agagCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGa                                      <  1:780558/151‑1 (MQ=255)
                                                                                                              agCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                                    <  1:1132264/151‑1 (MQ=255)
                                                                                                               gCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                                   <  1:1028995/151‑1 (MQ=255)
                                                                                                                cGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATTCAC‑ACACGGTAGCCCCCCCATCAGA‑AAACGGGGGGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGATGGTGGGTGTGAcc                                    >  1:29383/1‑149 (MQ=255)
                                                                                                                 gACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAAt                                 <  2:767892/151‑1 (MQ=255)
                                                                                                                  aCCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                                <  1:198319/151‑1 (MQ=255)
                                                                                                                     aTTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                             <  1:1082972/151‑1 (MQ=255)
                                                                                                                       tttGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                            <  1:1159495/150‑1 (MQ=255)
                                                                                                                       tttGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                            <  2:921095/150‑1 (MQ=255)
                                                                                                                         tgtgGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGa                          <  2:411425/150‑1 (MQ=255)
                                                                                                                             gcggtggcgCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt                     <  1:1019408/151‑1 (MQ=255)
                                                                                                                             gcggtggcgCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt                     <  1:1135350/151‑1 (MQ=255)
                                                                                                                              cggtggcgCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGCAGCGTCCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                       >  1:411088/1‑148 (MQ=255)
                                                                                                                                  ggcgCGC‑GGAAAACTTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTCGCTGTGACCAATGGTCGGATCGTTGCCGt                >  2:429092/1‑151 (MQ=255)
                                                                                                                                  ggcgCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg                 <  1:29894/150‑1 (MQ=255)
                                                                                                                                  ggcgCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt                <  1:1225080/151‑1 (MQ=255)
                                                                                                                                   gcgcgc‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg               <  2:856242/151‑1 (MQ=255)
                                                                                                                                    cgcgc‑gGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg               >  2:354857/1‑150 (MQ=255)
                                                                                                                                            aaaaaTTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc        >  2:208705/1‑150 (MQ=255)
                                                                                                                                             aaaaTTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑GCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt      <  2:944075/151‑1 (MQ=255)
                                                                                                                                             aaaaTTCAC‑ACACAGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt      <  2:628751/151‑1 (MQ=255)
                                                                                                                                               aaTTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCa    <  1:79328/151‑1 (MQ=255)
                                                                                                                                                aTTCAC‑ACACGGTAACCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGCATGGCGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAacac   >  2:879325/1‑151 (MQ=255)
                                                                                                                                                 ttCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACg  >  1:926832/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                         
TGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACG  >  NC_000913/3508518‑3508812

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: ‑