Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,140,400 G→A 8.3% R405W (CGG→TGG)  asmA ← putative assembly protein AsmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,140,4000GA 8.3% 180.1 / 2.9 60R405W (CGG→TGG) asmAputative assembly protein AsmA
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (3/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAATTC  >  NC_000913/2140256‑2140548
                                                                                                                                                |                                                                                                                                                    
gtgCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCt                                                                                                                                                 <  2:490923/150‑1 (MQ=255)
gtgCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGACt                                                                                                                                                 <  2:1297909/150‑1 (MQ=255)
   ccTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGa                                                                                                                                              <  2:754908/150‑1 (MQ=255)
      gTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGt                                                                                                                                           >  1:672060/1‑150 (MQ=255)
        ccTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCACGGCTGGaaaacg                                                                                                                                         <  2:549338/150‑5 (MQ=255)
         cTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTa                                                                                                                                        >  2:1255059/1‑150 (MQ=255)
          ttGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCACGGCTGGaaaacgtt                                                                                                                                       <  1:163768/150‑7 (MQ=255)
           tGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATc                                                                                                                                      >  2:1209956/1‑150 (MQ=255)
             ccAGTTGTGGCGGAATGCGTCGGCATCTATGTCAACACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCcg                                                                                                                                    >  2:186540/1‑150 (MQ=255)
               aGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCcgcg                                                                                                                                  <  1:723084/150‑1 (MQ=255)
                     ggCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTa                                                                                                                            >  1:1219046/1‑150 (MQ=255)
                       cGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGTCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAtt                                                                                                                          >  1:1385900/1‑150 (MQ=255)
                       cGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAtt                                                                                                                          >  2:1075064/1‑150 (MQ=255)
                          aaTGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGAt                                                                                                                       <  2:1124675/150‑1 (MQ=255)
                             gCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTc                                                                                                                    >  2:249912/1‑150 (MQ=255)
                               gTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTcgc                                                                                                                  >  1:872059/1‑150 (MQ=255)
                                           gTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGcc                                                                                                      <  2:1217331/150‑1 (MQ=255)
                                             cAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCg                                                                                                    <  2:872059/150‑1 (MQ=255)
                                                   cGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTg                                                                                              <  2:1219046/150‑1 (MQ=255)
                                                        aaaTCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAAcc                                                                                         <  2:311616/150‑1 (MQ=255)
                                                          aTCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTg                                                                                       <  1:1209956/150‑1 (MQ=255)
                                                             aCCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTcc                                                                                    <  1:295718/150‑1 (MQ=255)
                                                                agccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCAcc                                                                                 <  1:1255059/150‑1 (MQ=255)
                                                                 gccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCg                                                                                >  2:468056/1‑150 (MQ=255)
                                                                     gTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTaa                                                                            >  2:1097342/1‑150 (MQ=255)
                                                                       gACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGTATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGt                                                                          >  1:1418048/1‑150 (MQ=255)
                                                                         cATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCa                                                                                   >  2:1182275/1‑139 (MQ=255)
                                                                         cATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCa                                                                                   <  1:1182275/139‑1 (MQ=255)
                                                                            ttttCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGcc                                                                     <  1:418221/150‑1 (MQ=255)
                                                                             tttCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCATCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCTTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGcgc                                                                    <  2:720686/150‑3 (MQ=255)
                                                                                 cGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCa                                                                >  2:1065701/1‑150 (MQ=255)
                                                                                  ggTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAg                                                               >  2:1009288/1‑150 (MQ=255)
                                                                                  ggTCAACGAAATCGGATAGTTAAACGCCTTCAGAATAGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAg                                                               >  1:314444/1‑150 (MQ=255)
                                                                                      aaCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCTGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGa                                                           >  2:53814/1‑150 (MQ=255)
                                                                                       aCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAg                                                          >  2:1950/1‑150 (MQ=255)
                                                                                          aaaTCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTaa                                                       <  1:1075064/150‑1 (MQ=255)
                                                                                            aTCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAAtt                                                     >  1:516959/1‑150 (MQ=255)
                                                                                               ggATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTcc                                                  >  1:44534/1‑150 (MQ=255)
                                                                                                 aTAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAg                                                <  1:1950/150‑1 (MQ=255)
                                                                                                  tAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGc                                               >  1:119342/1‑150 (MQ=255)
                                                                                                      tAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTa                                                                             <  1:698140/116‑1 (MQ=255)
                                                                                                      tAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTa                                                                             >  2:698140/1‑116 (MQ=255)
                                                                                                            ccTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACt                                     <  1:1097342/150‑1 (MQ=255)
                                                                                                             cTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACtt                                    >  2:838404/1‑150 (MQ=255)
                                                                                                                cAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgtt                                 <  1:1065701/150‑1 (MQ=255)
                                                                                                                   aaTGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGt                              <  1:468056/150‑1 (MQ=255)
                                                                                                                       gTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCAtt                          >  2:1035309/1‑150 (MQ=255)
                                                                                                                         aCCAATCTCAACGTTTTCCAGCCGTGGATGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTg                        <  2:44534/150‑1 (MQ=255)
                                                                                                                              tctcAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTgg                   <  2:672060/150‑1 (MQ=255)
                                                                                                                                tcAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCa                 >  1:371242/1‑150 (MQ=255)
                                                                                                                                 cAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCaa                >  1:182473/1‑150 (MQ=255)
                                                                                                                                 cAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCa                 >  2:647370/1‑149 (MQ=255)
                                                                                                                                 cAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCa                 <  1:647370/149‑1 (MQ=255)
                                                                                                                                     gTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATc            <  2:516959/150‑1 (MQ=255)
                                                                                                                                     gTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATc            <  1:1109778/150‑1 (MQ=255)
                                                                                                                                     aacttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATc            >  1:405326/4‑150 (MQ=255)
                                                                                                                                      acttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCt           >  1:1238959/3‑150 (MQ=255)
                                                                                                                                       cttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTg          >  1:818527/2‑150 (MQ=255)
                                                                                                                                            aGCCGTGGCTGGAAGTTTATCCCCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTaaaa     >  2:725867/1‑150 (MQ=255)
                                                                                                                                               cGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAattc  <  2:119342/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAATTC  >  NC_000913/2140256‑2140548

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