Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 94,278 G→A A371A (GCG→GCA)  murE → UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091394,2780GA96.2% 153.1 / ‑2.3 53A371A (GCG→GCA) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (24/26);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

TTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCGGG‑AATGTTAGATGCCGGACATGCCAAAGTGATGGAAGGCCGTGCTGAAGCGGTGACTTGCGCCG  >  NC_000913/94033‑94519
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
tttAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGcgcg                                                                                                                                                                                                                                                >  1:364362/1‑251 (MQ=255)
 ttAGCTCAAGTGGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGcgcg                                                                                                                                                                                                                                                <  1:234821/250‑1 (MQ=255)
                   atggaaaaaTTGAAAGCCACTTGATGGGCGCTTTTCATGTCTGCACCCCGCTGCTCGCGCGGGCTACACGGTTGGCACTCGGCTATCCCCTGGCTGATCTGCTGAAACCCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCTGGCAAg                                                                                                                                                                                                                             <  2:115936/245‑1 (MQ=255)
                    tGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAgg                                                                                                                                                                                                                            >  1:184954/1‑250 (MQ=255)
                          aaTTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTgtggt                                                                                                                                                                                                                      <  1:356544/251‑1 (MQ=255)
                                 aGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGGGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGTGTCTGCACTGTGCGGGAAAGCTgtggagggtgtt                                                                                                                                                                                                               >  2:172275/1‑243 (MQ=255)
                                  gCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGAGACACTGTTGGCACTCGGCTCTCCACTGGCTGATCTGCTGAAAACCGCCGCGCGACTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAACCCGACA‑GTGGTGGTGGATTCCGCGCATACGCCGGATGCACTGGAAAACGCCTTCCAGGCAGCGAGGCAGCACTGTGCGGGCAGGCtgaggggtggattt                                                                                                                                                                                                              >  2:266966/1‑239 (MQ=255)
                                       cTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGGGTGTCTTTGGCtg                                                                                                                                                                                                         >  1:71330/1‑251 (MQ=255)
                                                  ttttAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATc                                                                                                                                                                                              >  1:193991/1‑251 (MQ=255)
                                                    ttAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGCTCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGACTTACAGGCAGCGCGTCTGCACTGGGCGGGCTCGCTGTGGTGTGTCTTTGGCTGGGGTGGCGGTcgg                                                                                                                                                                                            >  2:434309/1‑250 (MQ=255)
                                                      aactgcAGCAACCTGCTGCTCGCGCTGTCGACACTGTTGGCACTCGGCTATCCACTGGCGTATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCTCCAGGCAAACCGACG‑GTGGTGGTTGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATcgcg                                                                                                                                                                                           <  1:428747/245‑1 (MQ=255)
                                                        cgccAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCTCTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATa                                                                                                                                                                                        <  1:51286/248‑1 (MQ=255)
                                                                cccTGGTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACGGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTTCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAg                                                                                                                                                                                <  1:273529/250‑1 (MQ=255)
                                                                aCCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCACGATAAAGGTa                                                                                                                                                                                  >  1:280257/1‑249 (MQ=255)
                                                                  ctgctgCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCg                                                                                                                                                                              >  1:315444/1‑251 (MQ=255)
                                                                  cggctgctCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCg                                                                                                                                                                              <  1:103370/249‑1 (MQ=255)
                                                                    gctgctCGCGCTGGCGACACTGTTGGCACTCGGCTATCCAATGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GAGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCTGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATCAAGGTAAGCGTc                                                                                                                                                                            >  1:265341/1‑251 (MQ=255)
                                                                    gctgcgCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTc                                                                                                                                                                            <  1:465797/251‑1 (MQ=255)
                                                                        ctCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGGGGTGGCTTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCCAGCTGTGGTGTGTCTTTGGCTGGGATGGCGATCGCGCTAAAGGTAGGCGTCCACt                                                                                                                                                                        >  2:365857/1‑251 (MQ=255)
                                                                             gcTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATgg                                                                                                                                                                   <  1:362502/251‑1 (MQ=255)
                                                                             gcTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATgg                                                                                                                                                                   <  1:219106/251‑1 (MQ=255)
                                                                                ggCGACTCTGTTGGCACTCGGCTAGCCACTGGCTGAGCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATggg                                                                                                                                                                  <  2:71330/249‑1 (MQ=255)
                                                                                     cacTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACACCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCACGATGTGGTGTGTCTTTGGCGGGGGGGACGGTCGCGATGACGGGAAGCGTCCACTGAgggggccaagt                                                                                                                                                           >  2:118280/1‑244 (MQ=255)
                                                                                        tGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCACAATTGcc                                                                                                                                                        >  1:215117/1‑251 (MQ=255)
                                                                                                       aTCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGAc                                                                                                                                         >  1:400080/1‑251 (MQ=255)
                                                                                                         ccACGGGCTGATCGGCTGAAAACCGCCGCGCTTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACg                                                                                                                                        <  2:65567/250‑1 (MQ=255)
                                                                                                              ggCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCgg                                                                                                                                  >  1:176164/1‑251 (MQ=255)
                                                                                                                  gttctgcaagaaaCCGACGCGCGGCGGCTACCGGTTTGCGGGCGGATGGAA‑TGTGCCGCGGCTCCAGGCAACCCGACGCGTG‑TGGTGGATTACGCGCATACGCCGGATGCTCTGGAAAAAGCCTTACAGGCAGCGCGGCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGAGCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCggtgg                                                                                                                               <  2:55303/240‑1 (MQ=255)
                                                                                                                       gctgAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTga                                                                                                                            <  1:463156/248‑1 (MQ=255)
                                                                                                                          gAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGGGGTgacggacg                                                                                                                      >  1:302517/1‑251 (MQ=255)
                                                                                                                            aaaCCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGAt                                                                                                                    >  1:424355/1‑251 (MQ=255)
                                                                                                                                   cacgcgTCTGCAAACGGTGTGCGGACGCATGGAAGTTG‑TCGCTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACGGTTCGGGCTAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTTCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCcgc                                                                                                             <  2:323068/249‑1 (MQ=255)
                                                                                                                                     cgcgTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGGGGCGGACGCTAAACCGCGt                                                                                                           >  1:181505/1‑251 (MQ=255)
                                                                                                                                      gcgTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGGGGCGGACGGTAAcccccgg                                                                                                           >  1:38572/1‑246 (MQ=255)
                                                                                                                                          cTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGCTAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCggtgggggcagcagcaaacccgcgaacc                                                                                                        >  2:295675/1‑226 (MQ=255)
                                                                                                                                            gCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCg                                                                                                       <  1:1605/248‑1 (MQ=255)
                                                                                                                                                             gTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCatcatc                                                                                  <  1:302303/251‑1 (MQ=255)
                                                                                                                                                                agggaattggtCACTGCGCCAGGCAACCCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCTCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTGTGGCGGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTAACGAAGAACCGCG‑TGCCATCATCaa                                                                                <  2:365636/241‑1 (MQ=255)
                                                                                                                                                                               gcgcCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCgg                                                                  <  1:295675/250‑1 (MQ=255)
                                                                                                                                                                               gcgcCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCg                                                                   <  1:266221/249‑1 (MQ=255)
                                                                                                                                                                                 gcCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGGGGCGGACGATAAACCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCgggg                                                                >  1:183880/1‑249 (MQ=255)
                                                                                                                                                                                   cAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGGCGGGCGATAAACCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCgggga                                                               >  1:235547/1‑247 (MQ=255)
                                                                                                                                                                                        aaaCCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCGGG‑AATGTTAGa                                                       <  1:326137/251‑1 (MQ=255)
                                                                                                                                                                                                     gtggtggATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGAGGGGGCGGGAGATAACCCGCGAACCGGAAAAACGCGGTGCCA‑CATCACCGATATTTTGGCGGGGAATGtagaaggacggaacgg                                            >  2:38398/1‑236 (MQ=255)
                                                                                                                                                                                                          ggATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTTCAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGGTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATGCTGGCGGG‑AATGTTAGATGCCGGACATGCCaaa                                       <  2:461754/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggcATTCGCCGGATGCACTGGAAAAAGCGTTACAGGCAGCGCGTCTGCACTGTGGGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCTATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TCCCATCATCAACGATATTCTGGCGGG‑AATGTTAGATGCCGGACATGCCAAAGTGATGGaa                              <  1:266966/250‑1 (MQ=255)
                                                                                                                                                                                                                  cgcATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCGGG‑AATGTTAGATGCCGGACATGCCAAAGTGATGGaa                              <  1:71368/251‑1 (MQ=255)
                                                                                                                                                                                                                     aTACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGGCGGGCGGTAACCCCCGGACCGAAGAAACGCG‑GGCCCTCATCAACGGAATTCTGGCGGG‑GGTGGTAGATGTCGGAccaaccacaaggaggggaggg                           >  2:130528/1‑231 (MQ=255)
                                                                                                                                                                                                                        cGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGGTATTCTGGCGGG‑GATGTTTGATGACGGGCATGGCCAAGTGATGGAAGGCCgg                        >  1:256721/1‑250 (MQ=255)
                                                                                                                                                                                                                        cGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCGGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATGGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCGGG‑AATGTTAGATGCCGGACATGCCAAAGTGATGGAAGGCCGt                        <  2:310374/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCGGG‑AATGTTAGATGCCGGACATGCCAAAGTGATGGAAGGCCGTGCTGAAGCGGTGa           <  2:300429/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaaGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCGGG‑AATGTTAGATGCCGGACATGCCAAAGTGATGGAAGGCCGTGCTGAAGCGGTGACtt        <  1:130528/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCGGG‑GATGTTAGATTCCGGACATGCCCAAGAGATGGGAGGCCGTGCTGAAGGGGTGGCTTTCGCCg  >  1:415538/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
TTTAGCTCAAGTTGGGGCGATGGCGAAATTGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGT‑TCACTGCGCCAGGCAAACCGACG‑GTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCG‑TGCCATCATCAACGATATTCTGGCGGG‑AATGTTAGATGCCGGACATGCCAAAGTGATGGAAGGCCGTGCTGAAGCGGTGACTTGCGCCG  >  NC_000913/94033‑94519

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