Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 663,759 G→A R69C (CGC→TGC)  dacA ← D‑alanyl‑D‑alanine carboxypeptidase DacA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913663,7590GA97.4% 122.7 / NA 39R69C (CGC→TGC) dacAD‑alanyl‑D‑alanine carboxypeptidase DacA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (24/14);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCGGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATT  >  NC_000913/663515‑663993
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
aaaTCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGAc                                                                                                                                                                                                                                      >  1:88951/1‑251 (MQ=255)
     ggCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCATTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTg                                                                                                                                                                                                                                 >  1:134251/1‑251 (MQ=255)
       ccATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCg                                                                                                                                                                                                                               <  2:360379/251‑1 (MQ=255)
               acacAAGCAGCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATGTTCCTGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTc                                                                                                                                                                                                                       <  2:367582/251‑1 (MQ=255)
                                             aTACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAAt                                                                                                                                                                                         >  1:330370/1‑251 (MQ=255)
                                                  gcgGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTAACGGTGGCCCATGTGTCGTTGCCGATAGTGACTAAATCAGGTTCTTTAAAGTTACCGGCTTTCATTGCGTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGGTCGCGGCGGAAAGCTGCGTTCTGGTCGGCGAGCACTGTGCCGGAGTTATAgcgaatcaaga                                                                                                                                                                                    <  2:200086/251‑11 (MQ=255)
                                                    ggATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATg                                                                                                                                                                                  >  1:112928/1‑251 (MQ=255)
                                                        cAGCTGAGAAACCGGAACCTGCATGCACGGTTCGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCCTGCGTCGTCGCCGATAGTGACTAAATCAGTTTCTTTAAATTAACCTGCTTTCATTGCCTGGCCGATAACGTAACTGTCCATTATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCTGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAg                                                                                                                                                                               <  1:143929/250‑1 (MQ=255)
                                                           cTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACt                                                                                                                                                                           >  1:90544/1‑251 (MQ=255)
                                                                                  ccGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACAccc                                                                                                                                                    <  1:359476/251‑1 (MQ=255)
                                                                                        tGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATc                                                                                                                                              >  1:317408/1‑251 (MQ=255)
                                                                                             aaCATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCATCGTTCCCCATAGTGACTAAATCAGTTTCTTTAAATTTACCGGATTTCATTCCCTGCACGATAACGTAACTGGTCATCATTATAGTCAGGCTGGCAGGATCGCAGCGGACACACGCGTTTTGTTCGGCGAGCACTTTGCCAGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCAGGGATCATAGt                                                                                                                                         >  1:258105/1‑251 (MQ=255)
                                                                                                       aagaaCCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATGGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTc                                                                                                                               <  2:272887/251‑1 (MQ=255)
                                                                                                        agaaCGTTTAACCCCCGGGTTCCTGGTCGCCCATGCCTCGACGCCGATCATGACTACATCAGTTTCTTTAAATTTCCCGGATTTCATTGCCTGGCCCATAACGTAACTGCTCATCATTTTGGTCAGGCTGGCAGGATCGAAGCGCACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGGGCTATAGTCAATCAGGTTGTAGGACTCCGCATCGATCTGCGGTACCCCCGGGATCATAGTTTTGctatcc                                                                                                                               >  2:449483/1‑248 (MQ=255)
                                                                                                                tAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGAc                                                                                                                                                                                                                                      <  2:364970/139‑1 (MQ=255)
                                                                                                                tAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGAc                                                                                                                                                                                                                                      >  1:364970/1‑139 (MQ=255)
                                                                                                                tAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCg                                                                                                                      >  1:306529/1‑251 (MQ=255)
                                                                                                                       cGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCAtgtg                                                                                                               >  1:324216/1‑251 (MQ=255)
                                                                                                                            tACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCa                                                                                                                   >  2:395279/1‑242 (MQ=255)
                                                                                                                            tACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCa                                                                                                                   <  1:395279/242‑1 (MQ=255)
                                                                                                                                ggTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCCTAGTTTTGATATTCAGGTCATCGGCATGTgcaggcg                                                                                                         >  2:342256/1‑245 (MQ=255)
                                                                                                                                   ggCCCACGCGTCATTGCCGAAAGCGACTAAATAAGATTCTTTAAATTTACCGGCGTGCATTACCTGGCCGATAACGTCACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTGCTGTTCGGCGAGCACTTTGCCGGAGTTATAGGCAAGCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATTATAGTTGTGATATTCCGGACATCGTCATGTGCAGCAGAGATa                                                                                                    <  2:282275/250‑1 (MQ=255)
                                                                                                                                        aTGCGTCGTTGCCGATAGTGACTAAATAAGTTTCTTTAAATTTACGGGCTGTCATTGCCTGGCAGATAAGGTAACTGATCATAATTTTGGTCAGGCTGGCAGGATAGCAGCGGATATCTTCGTTCTGTTCGGCGAGCAATATGCCGGAGATATAGTCAATAAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATCGTTTTGATCTTCAGGGCATTGGCATGTGCAGTAGGGATaacagct                                                                                              >  2:139892/1‑246 (MQ=255)
                                                                                                                                                     gATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCg                                                                                 >  1:375873/1‑251 (MQ=255)
                                                                                                                                                        aGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCgtg                                                                               <  1:357305/250‑1 (MQ=255)
                                                                                                                                                               aaaTCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCCTCATTTTGGTCAGGCTGGCAGGAGCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTGATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCTGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAgc                                                                         <  1:283947/249‑1 (MQ=255)
                                                                                                                                                                   cAGTTTCTTTAAATTTACCGGCTTTCATGGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAgg                                                                   <  2:360465/251‑1 (MQ=255)
                                                                                                                                                                         atttAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTATGTTCGGCGAGCACTTTGCCGGAGTTATAGGCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACAACCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGGAAAAGGACGTGGAAAGAGACGGGGTGAGCGGCAGGCGCTTc                                                             >  1:265065/2‑251 (MQ=255)
                                                                                                                                                                             aaaTTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGGGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGGGAGCGCCAGGCGCTTCATgt                                                         >  1:343657/1‑250 (MQ=255)
                                                                                                                                                                              aaTTTACCGGCTTTCATTGCCTGGCCGATCACGTAACTGGTCATCATTTTGGTCAGGCTTGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGCCTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAAGATGTGCAAAGAGCCGGTGAGAGCGTCAGGCGCTCCAtta                                                         >  2:76775/1‑248 (MQ=255)
                                                                                                                                                                                 ttACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATAc                                                      <  2:274579/250‑1 (MQ=255)
                                                                                                                                                                                                          tAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGGATTCATGGTc                            >  2:14126/1‑251 (MQ=255)
                                                                                                                                                                                                             cGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGa                         <  1:11283/251‑1 (MQ=255)
                                                                                                                                                                                                                     gTCATCATTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGGGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTANGAc                 >  2:406496/1‑251 (MQ=255)
                                                                                                                                                                                                                      tcatcaTTTTGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACgcc                 >  1:354677/1‑248 (MQ=255)
                                                                                                                                                                                                                             tttGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAgg                                                                                                                            >  2:215302/1‑136 (MQ=255)
                                                                                                                                                                                                                             tttGGTCAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAgg                                                                                                                            <  1:215302/136‑1 (MQ=255)
                                                                                                                                                                                                                                   cAGGCTGGCAGGATCGCAGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATAAGAGCGGAAAAAATGGTATTCATGGTATGAACTACGACATCCGTGATGGAAt   >  1:239640/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aGGCTGGCAGGATCGCAGCGGACCTCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTACGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGACGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGGGTTCATGGTTTGAACAACGACATACGTTATGGNAtt  >  2:417448/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
AAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTGGTCAGGCTGGCAGGATCGCGGCGGACATCTGCGTTCTGTTCGGCGAGCACTTTGCCGGAGTTATAGTCAATCAGGATGTAGGACTCCGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATT  >  NC_000913/663515‑663993

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: