Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG31.3% 107.1 / 15.9 68I441V (ATC→GTC) mplUDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase
Reads supporting (aligned to +/- strand):  ref base A (3/43);  new base G (21/0);  total (25/43)
Fisher's exact test for biased strand distribution p-value = 1.56e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTG  >  NC_000913/4456865‑4457349
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                            
cgcTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCCTATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑GGGGGTCCATCAGa                                                                                                                                                                                                                                                     >  1:296889/1‑250 (MQ=255)
  cTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGG‑‑GGGGGCCATCAGA‑AA‑c                                                                                                                                                                                                                                                >  1:140335/1‑251 (MQ=255)
   tCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGCAAAACGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTG‑GGGGGTCCAACAGA‑AA‑c                                                                                                                                                                                                                                                >  1:56776/1‑250 (MQ=255)
      aaTACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGTGATCCGTCAGA‑AA‑ctgct                                                                                                                                                                                                                                            <  1:329247/251‑1 (MQ=255)
        tACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGTGGG‑TCCATCCGA‑AA‑CTGCTgg                                                                                                                                                                                                                                          >  1:142140/1‑251 (MQ=255)
         aCCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTgg                                                                                                                                                                                                                                          <  2:322776/250‑1 (MQ=255)
          ccATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGGCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTTTGCCGGGCGACCATATGCTGGTGATGACCAGCGGCGG‑TTTTGG‑GGGGATCCTGCAGT‑AA‑CTGCTGG‑At                                                                                                                                                                                                                                       <  2:424165/251‑1 (MQ=255)
            aTGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATgg                                                                                                                                                                                                                                     <  2:101067/251‑1 (MQ=255)
              gAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGAGATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCAt                                                                                                                                                                                                                                                         <  1:190429/232‑1 (MQ=255)
               aaaaTGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTTG‑TGGGGTCCCTCCGa‑aaacggtgggtggggg                                                                                                                                                                                                                                    >  1:296107/1‑237 (MQ=255)
                   tGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑CTTTTG‑TGGGGACCATCAAACAA‑AT‑CTGG‑ATGGGCTgggg                                                                                                                                                                                                                              >  1:423542/1‑249 (MQ=255)
                            gCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGaa‑gaag                                                                                                                                                                                                                       <  1:426432/248‑1 (MQ=255)
                                aGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAg                                                                                                                                                                                                                <  1:307645/251‑1 (MQ=255)
                                         ggCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAg                                                                                                                                                                                                        <  2:14030/250‑1 (MQ=255)
                                           cACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTaa                                                                                                                                                                                                     <  2:318803/251‑1 (MQ=255)
                                             ccTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGG‑GGGGGTCC‑TCACA‑AA‑CTGATGG‑GCGGGGTGGGGAAAGAAG‑CGGAAGCCGCGCAGTAAtt                                                                                                                                                                                                   >  1:267571/1‑251 (MQ=255)
                                                  aTTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCCGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCCTATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑GGGGGTCCATCCGC‑AA‑CTGCTGG‑TTGGGCTGGCGGA‑GAAGGGGGAAGGCGCGCCGTAATTCGGcc                                                                                                                                                                                              >  1:254723/1‑251 (MQ=255)
                                                     aGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTTGATACGCTGGCAGAGATGGTGGTGAAAACCGCTCAGCCTGGCGACCATAGTCTGGTGAGGATCAACGGCTG‑TTGTGG‑TTGGATCCATCAGA‑AA‑CTGCTGT‑ATGGGCGGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCtt                                                                                                                                                                                            <  1:189549/250‑1 (MQ=255)
                                                     aGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCtt                                                                                                                                                                                            <  1:180786/250‑1 (MQ=255)
                                                          gTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTG‑TGGGGTCCATTAAA‑AA‑ATGCTGG‑ATGGGGTGGGGAA‑GAAGGGGGAAGCCGCGCAGGAAATCgggccttgcc                                                                                                                                                                                        >  1:241422/1‑241 (MQ=255)
                                                             ccGATGAAGTCTTCCTTCTGCAACCGGCGCATATTCCGTGGCAGGTGTCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTTGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGata                                                                                                                                                                                   <  2:107029/251‑1 (MQ=255)
                                                                  gAAGTCGTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTa                                                                                                                                                                              <  2:142041/251‑1 (MQ=255)
                                                                          ccTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAgg                                                                                                                                                                      <  1:80858/251‑1 (MQ=255)
                                                                           ctgctgCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGG‑‑GGGGTTCATCCGA‑AA‑ATGCTGG‑ATGGGGGGGCGGA‑GAAGGGGGGAGCCCCGCCGTAATTCGGCCCTTTCCCGGTATAGTATGCGTAGGc                                                                                                                                                                     >  1:372958/1‑251 (MQ=255)
                                                                             gctgcAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGcc                                                                                                                                                                    <  1:227461/250‑1 (MQ=255)
                                                                                    ccGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATGA‑gg                                                                                                                                                          <  1:242009/251‑1 (MQ=255)
                                                                                           agatTCCGTGGCAGGTGGCAGAAGTGGCCGAAGCCTGCGTTCAGGCTGCACACTGGAGTGGCGATGGGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGAGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑AGAA‑GGCG‑TTCAc                                                                                                                                                  <  2:60768/249‑1 (MQ=255)
                                                                                             atTCCGTGGCAGGTGTCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCg                                                                                                                                                                   <  2:154269/235‑1 (MQ=255)
                                                                                               tCCGTGGCAGGTGGCAGAAGGGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAACCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCAcgcc                                                                                                                                               <  2:372958/250‑1 (MQ=255)
                                                                                                           ggcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGGGCCATCAGA‑AA‑ACGCTGGGATGGG‑TGGCGAA‑GGAAGCGGGAGCCGGGCCGTAAATAGGCCTTAGACAGGAATATTATGCGTGGGGCGG‑GTAA‑GGGG‑GTCACGCCGCACCCGGCCttt                                                                                                                                  >  1:307693/1‑251 (MQ=255)
                                                                                                               gaagtggcagaggCCGGCGTTCAGCCGGCACACTGGAGTGGCGAGGTGGAGACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCGGGGGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑c                                                                                                                             <  2:231218/251‑1 (MQ=255)
                                                                                                                        aaagCCTGCGTTCAGCCTGCACACGGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAT‑CATAGTGCCt                                                                                                                    <  1:350548/250‑1 (MQ=255)
                                                                                                                           gCCTGCGTTCAGCCTGCACCCTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGAAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CAGAGTGCCTGAt                                                                                                                 <  2:120919/251‑1 (MQ=255)
                                                                                                                           gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGAt                                                                                                                 <  1:140796/251‑1 (MQ=255)
                                                                                                                           gCCTGCGTTCAGCCTGCACACGGGAGTGGGGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAGCGGCGG‑GTTTGG‑TGGGAGCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTAGGCGTAGGCCGG‑ATAA‑TGCG‑GTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGAt                                                                                                                 <  2:214459/251‑1 (MQ=255)
                                                                                                                            ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGACAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTGG‑GGGGGTACATCAGA‑AA‑ACGATGG‑GTGGGCTGGGGGA‑GAAGGGGGAAGACGCGCAGGAATTCGGGCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGGGGTTC‑CCCCGCATCCGGCATTTggggcaaacgccttgtg                                                                                                                 >  1:266297/1‑235 (MQ=255)
                                                                                                                             cTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGGTCCATCAGA‑AAACT‑CTGG‑ATGGGCTGGCGGA‑AAAGGCGGAAGACGCGCCGTAAATCGGCCTTTGCCAGAAATAGTATTCGTAGGCCGGGTTAA‑GG‑G‑GTCACGCCGCCTCCGGCCTTTTAG‑CATAGTGCCCGATgg                                                                                                               >  1:160665/1‑250 (MQ=255)
                                                                                                                             cTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGAACATATTCTGGTGATGAGCAACGGCGG‑TTTTTG‑GGGGGTCCATCAGA‑AA‑AAGCTGG‑ATGGGGTGGCGCA‑GAAGGGGGGAGACGCGCAGTAATTCGGCCTTTGCCAGATATAGTATTCGTAGGCCGG‑ATAA‑GGGG‑TTTCCGCCGCATCCgggctttgtggctagtgcccgatgg                                                                                                                >  1:171445/1‑228 (MQ=255)
                                                                                                                                      gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGGTCCCTCAGA‑AA‑ATGCTGG‑ATGGGGTGGCGAA‑GAAGGCGGGAGCCGCGCAGGAATTTGGCCTTTGCCAGATATAATATGCGTAGGCCGG‑GTAA‑GGCG‑GTCACGCCGCATCCCGCCTTTTAG‑GATAGTGGCTGATGCGACGCTTgg                                                                                                      >  1:134124/1‑250 (MQ=255)
                                                                                                                                              acTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCtt                                                                                              <  2:145771/251‑1 (MQ=255)
                                                                                                                                                  gAGTGGCGATGTGGATCCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCa                                                                                          <  1:342929/251‑1 (MQ=255)
                                                                                                                                                  gAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTGGAGGCGTCTTATCa                                                                                          <  1:111388/251‑1 (MQ=255)
                                                                                                                                                    gTGGCGATGTGGATACGCTGGCAGAAATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGGTCCATCAGA‑AA‑CCGCTGG‑GGAGGGTGGCGGA‑GAAGGAGGAAGGCGCGCAGGAAATTGGGCTTTGCCCGAGATAGTATGCGTAGGCCGG‑GTAA‑GGCG‑GTCAACCCGCATCCGGGACTTGAGACA‑AGTGGCCGTTGCGACGCTTGATGCGTCTTATCAgg                                                                                        >  1:352561/1‑251 (MQ=255)
                                                                                                                                                                  aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGGTCCATCAGA‑AA‑CTGCTGG‑GTGGGCTGGCGAA‑GAAAGCGGGAGCCGCGCAGGAATTTGGCCTTAGCCCGATATAGTATTCGTAGGCCGG‑GTAA‑GGCG‑GTCACGCCGCATCCGGCATTTGAG‑CATAATGGCTGGTGCGACGCTTGATGGGTCTTTTCAGGCCcacaaatcgga                                                                           >  1:34600/1‑243 (MQ=255)
                                                                                                                                                                  aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCTGCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCGTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCAcc                                                                          <  1:417296/251‑1 (MQ=255)
                                                                                                                                                                  aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTAc                                                                                   <  2:34600/242‑1 (MQ=255)
                                                                                                                                                                       ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTAt                                                                      <  1:368752/250‑1 (MQ=255)
                                                                                                                                                                       ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTGG‑GGGGGTCCCTCCGA‑AA‑CTGCTGG‑AGGGGGTGGCGGA‑GAAAGCGGAAGCCGCGCAGCAAATTGGCCTTTGCCCGATATAGGATGGGGAGGCCGG‑GTAAAGGGG‑TCCA‑GCCGCATCCGGCATTTGAGGCTAAGGGCTTG‑TGCGGCGCTTGATGCGGCTTTTCAAGCCTGCAGGTGCACCGTagc                                                                     >  1:237916/1‑249 (MQ=255)
                                                                                                                                                                          aGATATGGTGGTGAAAACCGCGCAGCCTGGCGACCATATGCTGGTGATGAGCAACGGCGG‑TTGTGT‑TGGGATCCAGCAGA‑AA‑CTGCTGA‑ATGGGCTGGCGAA‑GAAGGCGAAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATTCGTGGTCCGG‑AAAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATCCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCgg                                                                  <  1:407110/251‑1 (MQ=255)
                                                                                                                                                                           gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTG‑GGGGGTCCATCCGA‑AA‑CCGCTTG‑GTGGGCTGGCGAA‑GAAGGCGGGAGCCGCGCAGTAATTTGGCCCTTGCCAGATATAGTATGCGTAAGCCGG‑GTAA‑GGCG‑TTCACCCCGCATCCGGCATTTGAG‑CCTAGTGCCTGGTGCGACGCTTGTTGCGTCTTATCAGGCCTCCAGGTGCACCGTATcccgc                                                                 >  1:321100/1‑248 (MQ=255)
                                                                                                                                                                              atGGTGGTGAAAACCGCGCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTGGG‑TGGGATCCAGCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TGCACGCCGCATCCGGCTTTGGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCATCGTATCCGGCaaa                                                              <  2:381663/251‑1 (MQ=255)
                                                                                                                                                                                ggtggtGAAAACCGCTCAGCCTGGCGACGATATTCTGGTGTTGAGCAATGGCGG‑TTGTGG‑TGGGAGCCATGAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑TAAGGCGGAAGCGCCGCAGTAAGTTGACTTGAGTGAGATAGGGTAGGCATAGGCCGG‑AGAA‑GGCG‑TTCGCGCCTGAGCCGGGATTTGAG‑GATGGTGCCTGAGGCGTCGCGTGAGTCGGCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAAcc                                                            <  1:197542/251‑1 (MQ=255)
                                                                                                                                                                                ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAAcc                                                            <  1:106316/251‑1 (MQ=255)
                                                                                                                                                                                  gggtgAAAACCGCTCAGCCGGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGGATGCGTAGGCCGG‑ATAA‑GGCG‑GTCACGCCGCATCCGGCATTTGAG‑CATAGTTCCTGAGGCGACGCGTGATGCGTCGTATGAGGCGTACAGGGGCAGCGTATCCGGCAAACCAt                                                          <  2:254723/250‑1 (MQ=255)
                                                                                                                                                                                           cctcTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTGGG‑TGGGAGCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGGAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CAGAGTGCCTGATGCGACGCTTGATGCGTCTTTTCAGGCCTACAGGTGCACCGTGTCCGGCAAGCCATACTCACATc                                                 <  2:88827/248‑1 (MQ=255)
                                                                                                                                                                                             gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑GGGGATCCATCAGA‑AA‑CTACTGG‑GCGGGCTGGCGCA‑GAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCCAACATAGTATGCGAAGGCCGG‑ATAA‑AGCG‑TTCCCGCCGCCTCCGGCATTTTGG‑AAAAGTGCCCGGTGCGCCGCTTGATGCGTCTTATCAGGCCCCCAGGGGCACCGCATCCGGCAAAACATACTCACATca                                                >  1:355096/1‑250 (MQ=255)
                                                                                                                                                                                                  gCCTGGCGACCATGGTCTGGTGCGGAGCAACGGCGG‑TTGTGG‑TGGGATCCATCAGA‑AA‑CGGCGGG‑AGTGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGGAGGCAGG‑ATAA‑GGCG‑TTCACGCCGCAGCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCGTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATTCTCACATCAACAACg                                          <  2:79407/251‑1 (MQ=255)
                                                                                                                                                                                                   ccTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGGTCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCCGTAATTCGGCCTTAGCCAGATATAGTATGGGGAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTTAG‑CATAGGGCCTGATGCGACGCTTGGTGCGTCTTATCAGGCCTAAAGGGGCACCGTATCCGGCAAAACATACTCACATCAACAACg                                          >  1:92865/1‑250 (MQ=255)
                                                                                                                                                                                                    cTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGGTCCATCAGA‑AA‑CTGCTGG‑GTGGGGTGGCGAA‑AAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑GTCACGCCGCATCCCGCATTTGAG‑CATAGTGCCTGGTGCGACGCTTTATGCGGCTTATTAGGCCTACCGGTGCACCGTATCCGGCAAACCCTAATCACATcaacaaaga                                         >  1:371221/1‑247 (MQ=255)
                                                                                                                                                                                                      ggCGACCATAGTCTGGTGAGGATCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCATTAATCCGGCCTTTGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATGTGAG‑CATAGTGTCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGaaa                                       <  2:321309/250‑1 (MQ=255)
                                                                                                                                                                                                         gACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGGTCCACCAGA‑AA‑CTGCTGG‑GTGGGCCGGCGAA‑GAAGGCGGAAGGCGCGCAGGAATTCGGCCTTAGGCAGAAATAGTATGAGTAGGGCGG‑ATTA‑GGCG‑GTCACGCCGCATCCGGCCTTTTAG‑CATAGGGCCCGATGCGGCGCTTGATGCGGCTTTTCAGGGCTACAAGGGGACCGTATCCGGCAAAACCTACTCACATCAACAACGAAAAtt                                    >  1:203870/1‑250 (MQ=255)
                                                                                                                                                                                                           ccATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGGCTTATCAGGCCTACAGGTGCACCGTAGCCGGCAAACCATACTCACATCAACAACGAAAATTAcc                                 <  1:209884/251‑1 (MQ=255)
                                                                                                                                                                                                             atatTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCt                               <  2:134124/251‑1 (MQ=255)
                                                                                                                                                                                                                            gCAACGGGGG‑TTGGGG‑TGGGAGCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAGGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATTTAGTATGCGTAGCCTGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTGGATGCGTGTGATCAGGCCTACAGGGGCACCGTTTCCGGCAAGCCATACTCACATCAACAACGAAAATTGGCCTGCGTTCGCCGCTAAc                <  2:56832/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ggcgg‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTcgcg           >  1:229306/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gg‑TTTTGG‑TGGGACCCTTCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGAGGACGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGCATGAGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGATTGATGCGACTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTCCCCTTCGTTCGCCGCTAACTCgcgc          >  2:254247/1‑249 (MQ=255)
                                                                                                                                                                                                                                          ggg‑tggGATCCATCGTA‑AT‑CTGCTGG‑AGGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACt   <  2:17366/250‑1 (MQ=255)
                                                                                                                                                                                                                                           gg‑tggGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTGAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTg  <  1:205272/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                            
CGCTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGG‑TGGGATCCATCAGA‑AA‑CTGCTGG‑ATGGGCTGGCGAA‑GAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGG‑ATAA‑GGCG‑TTCACGCCGCATCCGGCATTTGAG‑CATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTG  >  NC_000913/4456865‑4457349

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