Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,853,769 C→G 6.3% T150S (ACC→AGC)  mazG → modified nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,853,7690CG 6.3% 177.3 / 2.2 59T150S (ACC→AGC) mazGmodified nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base G (2/2);  total (31/28)
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

CCGACTGGCGAGCTGTACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACG  >  NC_002947/1853624‑1853909
                                                                                                                                                  |                                                                                                                                            
ccGACTGGCGAGCTGTACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCAGGTCCTGCTTGATGACGTGCCGGCAACCTTg                                                                                                                                          <  1:748905/150‑1 (MQ=255)
 cGACTGGCGAGCTGTACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGAGAGCGTCGTGCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGc                                                                                                                                         >  2:748905/1‑150 (MQ=255)
           tcTGTACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATc                                                                                                                               >  2:979115/2‑150 (MQ=255)
              gTACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATcgcg                                                                                                                            >  1:565294/1‑150 (MQ=255)
               tACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCtt                                                                                                                                           <  2:171533/134‑1 (MQ=255)
               tACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCtt                                                                                                                                           >  1:171533/1‑134 (MQ=255)
                  gcgcCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCa                                                                                                                        >  2:1316955/1‑150 (MQ=255)
                  gcgcCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCCGCGATAGCGCCGGCa                                                                                                                        >  1:205217/1‑150 (MQ=255)
                    gcCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGc                                                                                                                      >  1:724960/1‑150 (MQ=255)
                    gcCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCCGCGATAGCGCCGGCaag                                                                                                                      <  2:205217/150‑3 (MQ=255)
                              acGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCaa                                                                                                            >  1:451070/1‑150 (MQ=255)
                                gCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCaaaa                                                                                                          >  2:1081237/1‑150 (MQ=255)
                                      ccTGAGCGAAGCCCAGGTCAAGTCGCGATTTGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGc                                                                                                    <  1:714626/150‑1 (MQ=255)
                                       cTGAGCGAGGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTTCAAAAACGTGCg                                                                                                   <  2:1196174/150‑1 (MQ=255)
                                        tGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCgg                                                                                                  >  2:348268/1‑150 (MQ=255)
                                                      gTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAAGGTTGCCGgcac                                                                                                                                   >  2:93128/1‑101 (MQ=255)
                                                      gTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAAGGTTGCCGgcac                                                                                                                                   <  1:93128/103‑3 (MQ=255)
                                                      gTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTc                                                                                    >  2:1398861/1‑150 (MQ=255)
                                                      gTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGAGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTc                                                                                    >  2:820147/1‑150 (MQ=255)
                                                       tCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGc                                                                                                 >  2:221492/1‑136 (MQ=255)
                                                       tCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGc                                                                                                 <  1:221492/136‑1 (MQ=255)
                                                         aaGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGAc                                                                                 <  2:518077/150‑1 (MQ=255)
                                                           gTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTg                                                                               <  1:158065/150‑1 (MQ=255)
                                                            tCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTgg                                                                              >  1:1221294/1‑150 (MQ=255)
                                                                    ggAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGAAACAGCAGCTGTCTTTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCggcggc                                                                      >  1:604423/1‑150 (MQ=255)
                                                                         aaaTCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGc                                                                 <  2:682681/150‑1 (MQ=255)
                                                                         aaaGCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTCCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGc                                                                 <  1:98818/150‑1 (MQ=255)
                                                                          aaTCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGcc                                                                >  1:163316/1‑150 (MQ=255)
                                                                            tCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCgg                                                              >  2:780938/1‑150 (MQ=255)
                                                                                ggCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCt                                                          <  2:927278/150‑1 (MQ=255)
                                                                                 gCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATTGCGGGCAGCCAAACTGCAAAAACGTGGTGCCACTGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTg                                                         <  1:1316955/150‑1 (MQ=255)
                                                                                 gCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTg                                                         <  1:348268/150‑1 (MQ=255)
                                                                                  ccGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTgg                                                        >  2:880416/1‑150 (MQ=255)
                                                                                    gAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGAc                                                      >  1:995772/1‑150 (MQ=255)
                                                                                       gagCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAg                                                   <  2:163316/150‑1 (MQ=255)
                                                                                          cgcgcgGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTc                                                <  1:979115/150‑1 (MQ=255)
                                                                                          cgcgcgGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTc                                                <  1:880416/150‑1 (MQ=255)
                                                                                                 gagaAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGaa                                          <  2:565294/150‑1 (MQ=255)
                                                                                                 gagaAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCAGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGaa                                          >  1:532516/1‑150 (MQ=255)
                                                                                                         ccagccCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGa                                  >  2:92801/1‑150 (MQ=255)
                                                                                                          cagccCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAAAATTGCCACCGATATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGAt                                 >  1:248879/1‑150 (MQ=255)
                                                                                                              ccGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAg                             <  2:604294/150‑1 (MQ=255)
                                                                                                              ccGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAg                             >  1:604294/1‑150 (MQ=255)
                                                                                                                gAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGtg                           >  2:22468/1‑150 (MQ=255)
                                                                                                                     gCTGTCGCTGCTTGATGACGTGCCGGCAACCTTCCCTGCTGCATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCa                      <  1:784707/150‑1 (MQ=255)
                                                                                                                      cTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAg                     >  1:135889/1‑150 (MQ=255)
                                                                                                                       tGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAgg                    <  2:532516/150‑1 (MQ=255)
                                                                                                                          cGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGCCGTTGCCGGTGCTGGACACGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCa                 <  2:65564/150‑1 (MQ=255)
                                                                                                                            ctgctTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATg               <  1:820147/150‑1 (MQ=255)
                                                                                                                            ctgctTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATg               <  2:724960/150‑1 (MQ=255)
                                                                                                                             tgctTTATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAAACAAACCGCCAAAACGCGCTGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATgg              >  2:380559/1‑150 (MQ=255)
                                                                                                                              gctTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGc                                                           >  2:140045/1‑104 (MQ=255)
                                                                                                                              gctTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGc                                                           <  1:140045/104‑1 (MQ=255)
                                                                                                                              gctTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGc             <  2:452477/150‑1 (MQ=255)
                                                                                                                                ttGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCg           >  2:400783/1‑150 (MQ=255)
                                                                                                                                   atgaCCTGCCGGGAAGCTTGCCGGCGCTATTGCGGGCAGCGAAACTGCAAAACCGTGCGGCCGCGGTCGGTTTCGACTGGCCGGCGGCGTGGCCGGTGCTGGACAAGGTCCGCGAATAGCTGGATGAAGTGCTGCAGGCCATGGCCGACg        <  2:604423/150‑1 (MQ=255)
                                                                                                                                     gaCGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGt      <  2:135889/150‑1 (MQ=255)
                                                                                                                                        gTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGAc   >  1:1414323/1‑150 (MQ=255)
                                                                                                                                         cgCCGGCAAGGTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGCAGAACTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACg  >  2:1402635/2‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CCGACTGGCGAGCTGTACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACG  >  NC_002947/1853624‑1853909

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