Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG23.5% 129.8 / 13.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 (9/43);  new base G (16/0);  total (25/43)
Fisher's exact test for biased strand distribution p-value = 1.39e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCCGCGCTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAAT  >  NC_000913/4456860‑4457355
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                   
aGCCGCGCTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTTGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCTGTTGTGGTTGGATCCAt                                                                                                                                                                                                                                                               <  2:305786/251‑1 (MQ=255)
 gCCGCGCTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATc                                                                                                                                                                                                                                                              <  2:441076/251‑1 (MQ=255)
      gcTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTTGTGGGATACATccgca                                                                                                                                                                                                                                                          >  1:360885/1‑246 (MQ=255)
      gcTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGACTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGGGATGAGCCACGGCGGTTTTGGGGGGaaaccagcgaa                                                                                                                                                                                                                                                           >  1:206869/1‑240 (MQ=255)
      gcTCGAAGACCATGAAGATGGTGATCTGCAAAGACGATCTGGCACCTGCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAGCCGGAGTATAGTCCGTGGCAGGGGGCAGAAGTGGCAGAAGCCTGGGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATAGGGTGTTGGAAAGCGCTCAGACTGGCGTACAGATTCGGGTGATGAGCAATGTCGGGTGTTGTGGAATCCATCAGaaa                                                                                                                                                                                                                                                         <  2:506590/251‑1 (MQ=255)
           aaTACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATccaaaaatgct                                                                                                                                                                                                                                                    >  1:4540/1‑241 (MQ=255)
            aTACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGAGCCATCCGAAActgctg                                                                                                                                                                                                                                                  >  1:132063/1‑251 (MQ=255)
                cATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATAGTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCATGCGTTCAGCCTGCCCACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAGACCGCTCAGCCTGGCGACCATATTCTGGTGATGATCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATg                                                                                                                                                                                                                                              <  1:36644/251‑1 (MQ=255)
                 aTGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAatgcttggtgg                                                                                                                                                                                                                                             >  1:412142/1‑245 (MQ=255)
                 aTGAAAATGGGGATCTGCAAAGACGATCTGGCACCGTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTAAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATgg                                                                                                                                                                                                                                             <  1:150533/251‑1 (MQ=255)
                       aTGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGc                                                                                                                                                                                                                                       <  1:19815/251‑1 (MQ=255)
                                tGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTTCCTCACAAACTGGTTGGTGGGCTggggagaaaagg                                                                                                                                                                                                                              >  2:372248/1‑241 (MQ=255)
                                 gCAAAGACGATCTGGCTCCTTCATAAGGTCGGGCCGATGAAGTCTGCCTGCTGCAACCGGCGCATATTCCGTGGCAGTTGGCAGAAGTTGCAGAAGCATGCGTTGAGCCGGCACACTGGAGTGGTGATGTGGATACGCTGGCTGATATGGGGTGGAAAACCGCTCAGCCTGGCGACCAGATTCTGGTGGTTAGCAACGGCGGTTTTGGTGGGATCCTGCATAAACTGCTGGTTGGGGTGGGTAAGAAGGCg                                                                                                                                                                                                                            <  2:244790/251‑1 (MQ=255)
                                 gCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTTTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGaagaag                                                                                                                                                                                                                               <  2:482529/248‑1 (MQ=255)
                                 gCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCg                                                                                                                                                                                                                            <  2:48234/251‑1 (MQ=255)
                                 gCAAAGACGATCTGGCACCTTCAGTAGTTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGGGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCGCAGCCGGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTGGTTTGGATCCATCATAAACGGCTGGAAGGGCTGGCTAAGAAGGc                                                                                                                                                                                                                             <  2:238491/250‑1 (MQ=255)
                                  cAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGAGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATCCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGACCATCCGAAAATGCTGGATGGGcggggggagaaagcgg                                                                                                                                                                                                                            >  1:481232/1‑236 (MQ=255)
                                            cTGGCACCTTCATTAGGTCGTGCCGCGGAAGTCTTCCTGCTGCACCCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGACGCCTGCGTTCAGCCTGCACACTGGAGTGGCGAAGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGGCCATATTCTGGTGATGAGCAACGACGGGTTTGGTGGGGTCACACAGAAACTGCTGGCTGGGGTggaggcagaagcgggaacggcgcg                                                                                                                                                                                                                  >  2:359463/1‑229 (MQ=255)
                                              ggctcctgcatgagggcGTGCCGATGAAGTCTGCCTGCTGCAACCGGGGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTTCACACTGGATTGGCGATTTGGATACGCGGGCAGAGATGGTGGTGAAAACCGCTCAGCCTGGCGACCAGATGCTGGTGATGAGCAACGGCGGCTTTGGTGGGATCCATCAGAAACGGCTGGATGGGTTGGCGAAGGAGGCGGAAGcgtcgcag                                                                                                                                                                                                               <  2:225450/234‑8 (MQ=255)
                                              ggCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCa                                                                                                                                                                                                                <  1:215196/249‑1 (MQ=255)
                                               ccACCTTCGTTAGGTCGTGCCTATGAGGTCTGCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGGAGGAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGGTCAGCCTGGCGACCATATGCTGGTGATGAGCAACGGCGGTTTTTGTGTGATCCATCAGAAACTGCTGGATGTGCTGGCGAAGAAAGCGGAAGCCGCTCAGTa                                                                                                                                                                                                            <  2:66911/250‑1 (MQ=255)
                                                      caatAGGGCGTGCCGATGAAGTCGTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATGCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCGGGATGGGCTGTCGAAGAAGGCGGAAGGCGCGCAGTAATTCGGc                                                                                                                                                                                                     <  1:438021/248‑1 (MQ=255)
                                                           ggTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGCGCAACGGCGGTTTTTGGGGGATCCCTCATAAACTGCTGGGTGGGGTGGCGAAGAAAGCGGGAGCCGCGCCGTAATTTGGcccta                                                                                                                                                                                                 >  1:238661/1‑247 (MQ=255)
                                                                   cGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTTCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTCGCCCGata                                                                                                                                                                                         >  1:274556/1‑250 (MQ=255)
                                                                    gATGAAGTCTGCCTGCTGCAACCGGCGCATATTCCGTGGAAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGAGACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGAGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCGTAGCCAGAtata                                                                                                                                                                                       <  1:417499/251‑1 (MQ=255)
                                                                         aGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGCTATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCCTTCGACACTGCCGGGTGGGCTGGCGGAGGAGGCGGGAGGCCGGCAGGAAATTCGCCTTAGGCCGATATAGTATg                                                                                                                                                                                 >  2:323947/1‑251 (MQ=255)
                                                                            cctCCTGCTGCAACCGGCGCATAGTCCGTGGCAGGTGGCAGAAGGGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTTGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCg                                                                                                                                                                               <  2:350238/248‑1 (MQ=255)
                                                                                    tgcAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGGCCATCCGAAAAAGCTGGATGGGGTGGCGGAAAAAGGGGAAGCACCGGCGAAATTCGGCCTTAGCCAGGTATAGTATGCGGAGGGCggg                                                                                                                                                                      >  1:225073/1‑250 (MQ=255)
                                                                                         ccGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAgg                                                                                                                                                                 <  1:410702/251‑1 (MQ=255)
                                                                                           ggCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTGGGGGGGTCCATCAGGAAATGCTGGATGGGCTGGCGAAGAAGGCGGAAAACGCGCAGTAATTCGGGCTTATCCCGAAATAATAAGCGTAGGCCGGGAAAGGCg                                                                                                                                                               >  1:277306/1‑251 (MQ=255)
                                                                                               cATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGACCATCAGAAACTGCTGGGTGGGCGGGGGGAGGAGGCGGAAGGCGCGCAGGAATTTGGCCTTAGCCAGATATAGTATTCGGAGGCCGGGTAAGGCggtc                                                                                                                                                            >  1:514024/1‑247 (MQ=255)
                                                                                                    tCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                                       <  1:90701/250‑1 (MQ=255)
                                                                                                       gtggcaggttgtagAAGTGGCAGAAGCTTGCGTTCAGCCTTCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                                   <  2:269321/251‑1 (MQ=255)
                                                                                                         ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGTCCCTCAAAAAACGCTGGGTTGGGTGGGGGAAAAGGCGGGAGCCGCGCCGGAATTAGGCCTTCTCCCGATATAGTTTGCGTAGGCCGGGGAAGGCGTTCAcgccccctcc                                                                                                                                                 >  2:227411/1‑245 (MQ=255)
                                                                                                             ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCGGATATAGTATGCGTAGGCCGGATAAGGCTTTCACGCCGCATCCGGCa                                                                                                                                            <  1:438927/251‑1 (MQ=255)
                                                                                                                 gaagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTg                                                                                                                                        <  1:347511/249‑1 (MQ=255)
                                                                                                                                      gTTCAGCCTGCACACTGGAGGGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACg                                                                                                                  <  2:162924/251‑1 (MQ=255)
                                                                                                                                            ccTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAAAAAATGCTCGGTGGGCTGGCGAAGAAAGCGGAAGCCGCCCCGTAATTTAGCCTTAGACAGATTTAGTTTGCGTAGGCCGGGTAAGGCGTTTACGCCGCATCCCGCATTTGAGCATAGAGCCTGATGCGGAGCTTGAt                                                                                                           >  1:174889/1‑251 (MQ=255)
                                                                                                                                              tGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGGGAAACCAGCTCAGCCTGGCGAACATATTCAGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTTAGAAACTAATGGGAGGGCTGGGGAAGAAGGGGGAAGACGCGCCGGAAATTGGCCTTAGCCAGATATAGACTGCGGAGGCCGGATAAGGGGGTCACGCCGCATCCGGCGTTTGCGCCTAGTGACTGATGAGACGCTTGATGc                                                                                                         >  1:342534/1‑251 (MQ=255)
                                                                                                                                                 acacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCt                                                                                                                           <  1:349900/231‑1 (MQ=255)
                                                                                                                                                 acacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGAGCCATTAGAAACTGCCGGAAGGGCTGGGGGAGAAAGGGGAAGGCGCGCCGGAATTCGGCCTTTGCCAGATATAGGATTCGTAGGGCGGGTAAGGCGTGCAAGCCGCATCCGGCATTTGAGCATAGTGCCtagatcggaagaggggcgtgg                                                                                                       >  2:349900/1‑231 (MQ=255)
                                                                                                                                                        ggtgtcGATGTGTATACGCTGGCAGATATGGTGTTGAAAACCGCGCAGCCTGGCGACCATATTCTGGTGATTAGCATTGGCGGTTGTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATATTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAg                                                                                               <  2:468190/246‑1 (MQ=255)
                                                                                                                                                            gCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTAGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCt                                                                                           <  1:94975/251‑1 (MQ=255)
                                                                                                                                                                 gtgGATACGCTGGCAGATATGGTGGGGAAAACCGCACAGCCTGGCGACCATATTCTGGTGATGAGCATCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAGGGCGTGCACGCCGCATCCGGCATGTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAg                                                                                       <  2:430664/250‑1 (MQ=255)
                                                                                                                                                                  tgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCCTATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCAAAAAATTCTGGGTGGGCTGGGGAAGAAGGCGGAAGGCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGGGGAGGCCGGATAAAGGGGTCAAGCCGCCTCCGGCCTTTTGACATAGTGCCTGATGCGTCGCtttgttagccttattagtcccaccggg                                                                                      >  2:126732/1‑226 (MQ=255)
                                                                                                                                                                              cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAACAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCg                                                                         <  2:155070/251‑1 (MQ=255)
                                                                                                                                                                                     ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGAAGAAGGCGGAAGCCGCGCCGTAATTCGGCCTTTGCCAGATATAGGATGCGGAGGCCGGGTAAGGCGGTCACGCCGCATCCGGCATTTTAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATTAGGCCTACAGGGGCACCCTATCCGGCAAAcc                                                                  >  1:167538/1‑251 (MQ=255)
                                                                                                                                                                                      gttgggAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCa                                                                 <  1:163464/246‑1 (MQ=255)
                                                                                                                                                                                       tggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTTGGGGGGCCCCCCAAACACAGCTGGGTTGGGTGGGGCAGAAGGAGGGAGGCACGCAGGAATTCCACCTTAGCCAGATATAGTAATGGTAAGGCGGAGAAGGAGGTCAAGCCCCCTCCGGCATTTGTGACAAGTGGCAGATGCGACACGTGTATCGGCTTATCCGGCCCACAGGTGCACCGTATCCGGCCAACCa                                                                 >  2:194870/1‑250 (MQ=255)
                                                                                                                                                                                              aaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAATGCTGGATGGGCTGGCGAAAAAAGCGGGAGCCGCGCCGTAATTCGGCCTTAGCCAGAAATAGTAAGCGGAGGCCGGATAAGGCGTTCACCCCGCATCCGGCATTTGTGCCTAAAGCCTGGAGCGACGCTTGATGGGTCTTATTAGGCCTACAAGAGCACCGCAACCGGCAAACCCTACttcc                                                          >  1:37931/1‑248 (MQ=255)
                                                                                                                                                                                                   gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATcaa                                                     <  1:404097/251‑1 (MQ=255)
                                                                                                                                                                                                      cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATcaacaa                                                  <  1:435279/251‑1 (MQ=255)
                                                                                                                                                                                                      cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATcaaca                                                   <  1:477655/250‑1 (MQ=255)
                                                                                                                                                                                                      cAGCCTGGCGACCATATTCTGGTGAGGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCAGACTCACATcaacaa                                                  <  2:379843/251‑1 (MQ=255)
                                                                                                                                                                                                        gccgggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACg                                                <  1:88241/247‑1 (MQ=255)
                                                                                                                                                                                                         ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTTCATCAGAAACAGATGGCTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAAGTTGGCCTTAGCAAGGTGTTGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCATTTGAGGCGAATGCCTGATGGGACGCTTGGTGCGGCCTATTTGGCCTACAGGGGCACCCTATCCCGCAAACAATACCCCCTTCAACAACg                                                >  1:378680/1‑250 (MQ=255)
                                                                                                                                                                                                            ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGaa                                              <  1:234317/249‑1 (MQ=255)
                                                                                                                                                                                                                     atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCtt                                    <  1:429274/250‑1 (MQ=255)
                                                                                                                                                                                                                     atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTCATGGGTGGGGTGGCGAAGAAGGCGGGAACCGCGCAGTAATTTGGCCTTCGCCAGATATAGTATGCGTAGGGCGGGTAAGGCGGTCACGCCCCATCCGGCATTTTAGCCTATCGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCGCCGCATCCGGCAACCCATACTCACCTCAACAAAGaaactttcccctc                                   >  1:261844/1‑241 (MQ=255)
                                                                                                                                                                                                                     atTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCAACAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTTGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTCCAGGGGGACCGTATCCGGGCAACCCTACTCACATCAAAAACGAAAATTAACCTTc                                   >  1:60849/1‑251 (MQ=255)
                                                                                                                                                                                                                       tCTGGTGATGAGCAATGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCAGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCg                                  <  1:468806/250‑1 (MQ=255)
                                                                                                                                                                                                                           gTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTtctc                             <  1:483724/251‑1 (MQ=255)
                                                                                                                                                                                                                             gttgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACAACAACAACGAAAAGTACCCTTCGTTCTCCg                           <  1:432706/249‑1 (MQ=255)
                                                                                                                                                                                                                             gatgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCg                           <  1:323947/251‑1 (MQ=255)
                                                                                                                                                                                                                              atgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGc                          <  1:497564/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ggtttggttggtaTCCATCGAAATCTGCTGGCTGGGCTGGCTAAGAAGGCGAAAGCCGCGCAGTAATTCGTCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGTGTTGACGCCGCATCCGGCATTTGCGCATAGTGCCTGATGCGACGCTTGATTCGTCTTACCAGGCCTACAGGCGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAg              <  2:265137/239‑1 (MQ=255)
                                                                                                                                                                                                                                              tttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGGGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGAt            <  2:384239/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTg        >  2:263909/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    ggggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGa      <  1:238986/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       gATCCATGAGAAACTGCTGGATGGGCTGGCGAAGAAGGCTGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATAGAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGaaaa   <  2:277062/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        aTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAAt  >  1:73459/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                   
AGCCGCGCTCGAATACCATGAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAAT  >  NC_000913/4456860‑4457355

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: