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

GGAGCCGCGCTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAAT  >  NC_000913/4456858‑4457355
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                            
ggAGCCGCGCTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATcc                                                                                                                                                                                                                                                           <  2:270754/251‑1 (MQ=255)
                 ccATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGG‑GGGATCCATCAGAAACTGCTGGA‑t                                                                                                                                                                                                                                         >  1:23601/1‑251 (MQ=255)
                 ccATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTTACTCAGAAAATGCTGGA‑t                                                                                                                                                                                                                                         >  1:62658/1‑251 (MQ=255)
                  cATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTCCCACAGAAAATGCTGGA‑Tg                                                                                                                                                                                                                                        >  1:192660/1‑251 (MQ=255)
                             ggATCTGCAAAGACGATCTGGCACCTTCATGAGGTCGGGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCGGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCGCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCgaag                                                                                                                                                                                                                             <  2:495589/251‑1 (MQ=255)
                               aTCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGGGGGGTCCCACAGGAACTGCTGGG‑TGGGCGGGCgaagg                                                                                                                                                                                                                            >  1:284761/1‑249 (MQ=255)
                                   gCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTCCATCAGAAACAGCTGGA‑Tggggggggggagaagggg                                                                                                                                                                                                                       >  1:359935/1‑236 (MQ=255)
                                     aaaGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTACGGCGGCAACCGTCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGa                                                                                                                                                                                                                     <  2:423144/251‑1 (MQ=255)
                                          cGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGTGGGGTCCCCCCGAAACTGCTTGA‑TGGGGTggggaggacagcgggagccg                                                                                                                                                                                                                >  1:327281/1‑233 (MQ=255)
                                               tGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTGGTGGGGTCCCTCAGAAACAGCTGGG‑TGGGCTGGCGGAGAAAGCGGAAGCCgcgc                                                                                                                                                                                                             >  1:252767/1‑249 (MQ=255)
                                                ggCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAGCCGGCGCATATTCCGTGGCAGGTGGAAGAAGTGGCAGAAGCCTTCGTTCAGCCTGCACACTGGAGTGGCGAGGTGGATACGCTGGCAGATAGGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTGTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCcgc                                                                                                                                                                                                               <  2:41342/246‑1 (MQ=255)
                                                   cccTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAGAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCTG‑TTTTGGTGGGATc                                                                                                                                                                                                                                                            <  2:331539/198‑1 (MQ=255)
                                                   aCCTTCATTAGGTCGTGCAGATGAAGTCTTACTGCTGCAACCGGCGCATATTCAGTGGCAGGAGGAAGCAGTGGCAGAAGACTGCGTTCAGCCTGCAAACTGGAGTGGCGATGAGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTGGGTGGGATAC‑TCAGAAACTGGCGGG‑TGGGCTGGGGGAGAAGGCGGGAGCcggcccggaa                                                                                                                                                                                                        >  2:515595/1‑242 (MQ=255)
                                                    ccTTCATTAGGTCGGGCCGGTGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGTCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAAtt                                                                                                                                                                                                      <  1:350523/251‑1 (MQ=255)
                                                           tAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGTGGGGTCCATCAGAAAACGCTGGGGTGGG‑TGGGGAAGAAAGCGGGGGCCGCGCAGGAAATCGG‑ccct                                                                                                                                                                                              >  2:517231/1‑249 (MQ=255)
                                                             ggTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAg                                                                                                                                                                                            <  1:292606/251‑1 (MQ=255)
                                                              gTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTTGATTCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGTCGGAAGCCGCGCAGTAATTCGG‑CCTTAGc                                                                                                                                                                                           <  2:284761/251‑1 (MQ=255)
                                                                               ttCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGACCCTCAGAAAATGCTGGG‑TGGGGTGGGGGAGAAGGCGGAAGCCGCGCCGGAATTCGG‑GCTTAGCCCGATATAGTATGCGt                                                                                                                                                                           >  1:79125/1‑250 (MQ=255)
                                                                                tCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTa                                                                                                                                                                          <  1:109701/250‑1 (MQ=255)
                                                                                accTGCTGCAACCGGCGCATGTTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACTGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATAGAGTATGCGt                                                                                                                                                                           <  2:29775/248‑1 (MQ=255)
                                                                                 ccTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGGGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCGGGCGAAGAAGGCGGAAGCCGCGCAGTAATGCGG‑CCTTAGCCAGATATAGTATTCGTAgg                                                                                                                                                                        <  2:229259/251‑1 (MQ=255)
                                                                                       gcAACCGGCGCATATTCCTTGGCGGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGTTATAGTATGCGTAGGCCGGAt                                                                                                                                                                  <  2:287944/251‑1 (MQ=255)
                                                                                        cAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                                                                                                                                 <  2:456462/251‑1 (MQ=255)
                                                                                         aaCCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCCGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGTGGGATGCATTAGAAACCGGTGGG‑TGGGCTGGCGAAGAAGGGGGAAGCCCCGCAGGAAATTGG‑CCCTAGCCAGATATAGAATGGGTAAGCCgggta                                                                                                                                                                 >  1:23356/1‑247 (MQ=255)
                                                                                          aCCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGG‑TGGGGTGGCGAAGAAGGCGGGAGGCGGGCAGGAATTCGG‑CCCTTGCCAGATATAGTATGCGTAGGGCGGATa                                                                                                                                                                 >  1:37460/1‑249 (MQ=255)
                                                                                            cGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGc                                                                                                                                                             <  1:368224/251‑1 (MQ=255)
                                                                                                gcATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTc                                                                                                                                                         <  1:216509/251‑1 (MQ=255)
                                                                                                  atatTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGCCGTTCAc                                                                                                                                                       <  1:134935/251‑1 (MQ=255)
                                                                                                                      gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGGTCCATCAGAAAATGCTGGG‑TGGGGTGGGGGAGAAGGCGGGAGCCGCGCCGGAAATCGGGCCTT‑GCCCGATATAGTATGCGGAGGCCGGGTAAGGCGGTCCAGCCCCATCCGGCCttttggc                                                                                                                                   >  1:452180/1‑247 (MQ=255)
                                                                                                                               gaagCCTGCTTTCAGCCTGCACACTGGGGTGGCGATGTGGATCCGCTTGCAGATATGGTGGTGAAAACCGCTTAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGTTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATTTAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCt                                                                                                                          <  1:74855/251‑1 (MQ=255)
                                                                                                                                            aGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGTTCCATCAGAAACTGCGGGG‑TGGGGTGGGGGAGAAGGGGGGAGCCGGGCAGAAAATTGG‑CCATAGCCAGATATAATCTGGGCAAGCCGGGTAAAGCGTTTACGCCGCATCCAGCAATTGGGGATAATGCCTGGTGGGACGcctg                                                                                                             >  1:425121/1‑248 (MQ=255)
                                                                                                                                                 gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGTATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAATAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCGGATATAGTATGCGTAGGCCGGATAATGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGc                                                                                                         <  1:474946/250‑1 (MQ=255)
                                                                                                                                                       tGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCCGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGGGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCCGCATTTGAGCATAGTGCCCGGTGGGACGCTTGATGCGGCTTAt                                                                                                  >  1:495679/1‑251 (MQ=255)
                                                                                                                                                          ggTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAg                                                                                               <  2:452180/250‑1 (MQ=255)
                                                                                                                                                                 aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGTGGGGGCCATCAGAAACTGCTGGA‑TGGGCTGGCGGAGAAGAGGGGAGACGCGCAGTAATTCGG‑GCCTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCATCCGGCATTTGAGCATAATGACTGATTCGCCGCTTGATGCGTCTTATCAAGcccac                                                                                         >  1:81699/1‑247 (MQ=255)
                                                                                                                                                                 aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACa                                                                                        <  1:271925/251‑1 (MQ=255)
                                                                                                                                                                 aTGTGGATACGCTGGCAGATAGGGGGGTGAAAACCGGTCAGCCTGGCGACCATATTCTGGGGATGAGAAACGGCGG‑GTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCTTGATGCGACGCTTGATGCGTCTTATCAGGCCTACa                                                                                        <  2:327281/251‑1 (MQ=255)
                                                                                                                                                                  tgtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTCCATCAGAAAACGCTGGG‑TGGGCGGGCGGAGAAGGCGGAAGACGCGCCGAAATTTGG‑CCTTAGACCGATATAGAATGAGGAGGCCGGATAAGGGGTTCACGCCGCATCCGGCATTTGTGCATAGTGCCTGCTGCGACGCCTGATGCGGCtttacaagcctaccg                                                                                       >  1:333870/1‑238 (MQ=255)
                                                                                                                                                                      gATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTCCATCCGAAACTGCTGGG‑TGGGCTGGCGAAGAAGGCGGAAGGCGCGCCGGAATTCGG‑CCTTAGCCCGATATAGTATTCGGAAGGCGGATAAGGCGGTCACGCCGCCTCCCGCATTTGAGCATAGGGCCTGATGGGACGCTTGATGCGTCTTTTCCGGCCTACAggggc                                                                                   >  1:432349/1‑248 (MQ=255)
                                                                                                                                                                            cTGGCAGATATGGTGGTGAAAACCGCTCAGCATGGCGACCATATTCTGGTGATGAGCAACGGCGG‑ATTTGGTGGGGTCCATAAGAAAATGCTGGG‑TGGGGTGGCGCAGAAGGGGGGAGCCGGGCAGAAATTTGG‑CCCTAGCCAGATATAGAATGAGTAGGCCGGGTAAGGGGGTCAAGCCGCCTCCGGCATTTTAGCATAGAGGCTGTTGCGAAGCTTGATGCGTCTTATAAGGCCTAccggggcaacgt                                                                              >  1:225034/1‑239 (MQ=255)
                                                                                                                                                                              ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATc                                                                           <  1:405110/251‑1 (MQ=255)
                                                                                                                                                                                 aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCgg                                                                        <  1:185311/251‑1 (MQ=255)
                                                                                                                                                                                    tatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAAGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGGCTTATCCGGCCTACAGGTGCACCGTATCCCGCa                                                                      >  1:154803/1‑250 (MQ=255)
                                                                                                                                                                                       ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCAGCAAAc                                                                   <  1:148067/250‑1 (MQ=255)
                                                                                                                                                                                          ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAATCTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCAt                                                                <  1:172163/250‑1 (MQ=255)
                                                                                                                                                                                                   cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATca                                                      <  1:517231/251‑1 (MQ=255)
                                                                                                                                                                                                       cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATcaaca                                                   <  1:450051/250‑1 (MQ=255)
                                                                                                                                                                                                            tGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTGGGGGGAGCCATCCGAAACTGCTGGG‑GGGGCTGGGGGAGAAAGGGGGAGCCCCGCAGGAAATCGG‑CCTTTGCCAGGAATAGTATTCGGAGGCCGGGTAAGGCGGTCACGCCCCATCCGGCATTTGTGCCTAATGCCCGTTGTGACCCTTTAAGCGTCTTTTCAGGCCTACCAGTGCACCCCATCCGGCAAACCAAACtataataccccaccaa                                              >  2:424334/1‑235 (MQ=255)
                                                                                                                                                                                                                gaccagatattGGTGATTAGCAACGGCGG‑TTTTGGTGGTGTCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCTG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTGGATGCGTCTTATCAGGCCTACAGGTGGGCCGTAGCCGGCAAACCATACTTACATCAACAACTAAAATTa                                         <  2:225034/241‑1 (MQ=255)
                                                                                                                                                                                                                 aCCATATTCTGGTGATGATCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCTTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTTCACCGTAGCCGGCAAACCATACTCACATCAACAACGAAAATTAc                                        <  2:81699/251‑1 (MQ=255)
                                                                                                                                                                                                                   cATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTAccc                                      <  1:424334/251‑1 (MQ=255)
                                                                                                                                                                                                                     tatTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAAtt                                          <  1:321218/245‑1 (MQ=255)
                                                                                                                                                                                                                      atTCTGGGGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAGCCATACTCACATCAACAACGAAAAGTACCCtt                                    <  2:23601/250‑1 (MQ=255)
                                                                                                                                                                                                                           ggTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGtt                                <  1:89418/249‑1 (MQ=255)
                                                                                                                                                                                                                                       cggcgg‑TTTTGGTGGGGTCCCTCAAAAACTGCTGGA‑TGGGCTTGCGAAGAAGGCGGGGACCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGGATGCGTAGGCCGGGTAAGGGGTTTACGCCGCATCCGGGATTTGTGCATAGTGGCTGATGCGACGCTTGCTGTGGCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCAAACTCACATCAACAAAGGAAATTACCCTTCGTTCTCCGCtcaatccc                  >  2:418360/1‑244 (MQ=255)
                                                                                                                                                                                                                                           gg‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTCTCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCgcgc                >  2:509031/1‑249 (MQ=255)
                                                                                                                                                                                                                                                tgggtggGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCACAGATa           <  1:150869/248‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggtggGATGTATCAGAAATTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCATTAATTCGG‑CCTTAGCTAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGGGCCGTATGCGACGCTTGTTGCGTCGTATCAGGCCTATAGGTGCGCCGTATCCGGCAAACCAGACGCACATCAACAACGAAAATTACCCTTCGTGCTCCGCTACCTCGCTCAGATACTg        <  2:517538/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGCTATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGaa     >  2:411128/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      ggATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGaaa    <  1:179373/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       gATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGaaa    >  2:385032/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       gATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATCGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGaaa    >  2:400896/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        atccatccGAAACTGATGGA‑CGGGCTTGCGCAGAACGCGGAAGCAGCTCAGTAAATCGG‑TCTTAACCCGACATAGTATGAGTACGCCGGATACGACGTTCACGACGCATCCGGCATTTAAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACCGGTGCACCGTCTCCGGCCAACCCTACTCCCATCAACAACGAAAATTACCCTTCGTTCCCCGCTAACTCGCGCCGATCCTggcaact  >  2:318968/1‑246 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                            
GGAGCCGCGCTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGA‑TGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGG‑CCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAAT  >  NC_000913/4456858‑4457355

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