Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,853,787 G→C 6.2% R156P (CGG→CCG)  mazG → modified nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,853,7870GC 6.2% 194.1 / 2.8 63R156P (CGG→CCG) mazGmodified nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base C (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACGCCGATGCGCTGGAAGATGAAGTCG  >  NC_002947/1853639‑1853933
                                                                                                                                                    |                                                                                                                                                  
tACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCgg                                                                                                                                                   <  1:88836/150‑1 (MQ=255)
 aCGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATTTCGTGCCGGCAACCTTGCCGGCGCTATCGCggg                                                                                                                                                  >  2:102652/1‑150 (MQ=255)
     gcCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCCGCGATAGCGCCGGCaag                                                                                                                                              <  2:868897/150‑3 (MQ=255)
     gcCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCCGCGATAGCGCCGGCaag                                                                                                                                              >  1:868897/1‑148 (MQ=255)
     gcCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCCGCGATAGCGCCGGCaag                                                                                                                                              <  1:323756/150‑3 (MQ=255)
                   ccAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAAc                                                                                                                                <  1:768724/150‑1 (MQ=255)
                    cAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACg                                                                                                                               >  2:955888/1‑150 (MQ=255)
                        cTGAGCGAAGCCAAGGTCACGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCg                                                                                                                           >  2:1149608/1‑150 (MQ=255)
                         tGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCgg                                                                                                                          >  1:1010833/1‑150 (MQ=255)
                              gAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCCGCCAAACAGCACATACGTGCGcccgcc                                                                                                                     >  2:713253/1‑147 (MQ=255)
                                  cccAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCAcggtcg                                                                                                                 <  1:631561/150‑1 (MQ=255)
                                   ccAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGAGGACGTGCCGGCAACCTTGCGGGCGCAATCGCGAGCAGCCAAACTGCAAAAACGTGCGGCCAggggcgg                                                                                                                >  1:414125/1‑150 (MQ=255)
                                    cAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCAcggtcggt                                                                                                               >  2:470101/1‑150 (MQ=255)
                                      ggTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGttt                                                                                                             <  1:713253/150‑1 (MQ=255)
                                        tCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAACTGTCGCTGCTGGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCg                                                                                                           >  2:988546/1‑150 (MQ=255)
                                           aGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAGCAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACt                                                                                                        >  1:1055870/1‑150 (MQ=255)
                                           aGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACt                                                                                                        >  2:1276262/1‑150 (MQ=255)
                                             tCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTgg                                                                                                      >  2:116208/1‑150 (MQ=255)
                                               gcgcTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGcc                                                                                                    <  2:129636/150‑1 (MQ=255)
                                                 gcTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCcgg                                                                                                  >  1:919186/1‑150 (MQ=255)
                                                   tGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGa                                                                                                          <  2:1198876/140‑1 (MQ=255)
                                                   tGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGa                                                                                                          >  1:1198876/1‑140 (MQ=255)
                                                      gaagaaATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGgcggcg                                                                                             >  2:984831/1‑150 (MQ=255)
                                                       aagaaATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGt                                                                                            >  1:34733/1‑150 (MQ=255)
                                                        agaaATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCggcggc                                                                                              <  1:545607/147‑1 (MQ=255)
                                                        agaaATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCggcggc                                                                                              >  2:545607/1‑147 (MQ=255)
                                                        agaaATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGtt                                                                                           <  1:470101/150‑1 (MQ=255)
                                                         gaaATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCg                                                                                                                           >  2:405097/1‑117 (MQ=255)
                                                         gaaATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCg                                                                                                                           <  1:405097/117‑1 (MQ=255)
                                                         gaaATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTg                                                                                          >  2:511958/1‑150 (MQ=255)
                                                          aaaTCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGc                                                                                         <  2:309193/150‑1 (MQ=255)
                                                           aaTCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGcc                                                                                        >  1:572055/1‑150 (MQ=255)
                                                              cAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGt                                                                                     >  1:89729/1‑150 (MQ=255)
                                                                aGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGc                                                                                   <  1:116208/150‑1 (MQ=255)
                                                                 ggCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCt                                                                                  >  1:1213088/1‑150 (MQ=255)
                                                                    cGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGACCGGCGGCGTTGCCGGTGCTGGa                                                                               >  1:322604/1‑150 (MQ=255)
                                                                          gcgcgcGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGt                                                                         <  1:955888/150‑1 (MQ=255)
                                                                                 gagaAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGaa                                                                  <  1:511958/150‑1 (MQ=255)
                                                                                     agagCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGc                                                              <  2:21432/150‑1 (MQ=255)
                                                                                      gagCCAGCCCGAACAGCTGTCGCTGCTTGATGAACTGGAGTCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCt                                                             <  1:66784/150‑1 (MQ=255)
                                                                                        gccagccCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCACCCTTGCCGACGCCATCCCGAAAAGCCACACCGCAAATTCGTGCGTCCATGGTCGTTTTCGACTGGCCGGCGGGGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTgg                                                           <  2:467745/150‑1 (MQ=255)
                                                                                          cagccCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCTAAGAGCTGGAt                                                         <  2:89729/150‑1 (MQ=255)
                                                                                          cagccCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGAt                                                         <  2:1194580/150‑1 (MQ=255)
                                                                                             cccGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCt                                                             <  2:737268/143‑1 (MQ=255)
                                                                                             cccGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCt                                                             >  1:737268/1‑143 (MQ=255)
                                                                                             cccGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGaa                                                      >  1:894146/1‑150 (MQ=255)
                                                                                                         tCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGc                                                                                   <  1:782182/109‑1 (MQ=255)
                                                                                                         tCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGcc                                          >  2:361178/1‑150 (MQ=255)
                                                                                                         tCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGACGGTGc                                                                                   >  2:782182/1‑109 (MQ=255)
                                                                                                            ctgctTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCCGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGATGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATg                                       >  1:1223469/1‑150 (MQ=255)
                                                                                                               ctTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGcc                                    <  2:894146/150‑1 (MQ=255)
                                                                                                                  gatgaCGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGAc                                 >  2:1131535/1‑150 (MQ=255)
                                                                                                                        gTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGAc                           >  2:322717/1‑150 (MQ=255)
                                                                                                                        gTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGAc                           >  2:1228332/1‑150 (MQ=255)
                                                                                                                              gCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACGCCGAt                     <  2:1010833/150‑1 (MQ=255)
                                                                                                                              gCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACGCCGAt                     <  1:1276262/150‑1 (MQ=255)
                                                                                                                                        cGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACGCCGATGCGCTGGAAg           <  1:102652/150‑1 (MQ=255)
                                                                                                                                          gcgcTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAg                                                                           >  2:708011/1‑84 (MQ=255)
                                                                                                                                          gcgcTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAg                                                                           <  1:708011/84‑1 (MQ=255)
                                                                                                                                          gcgcTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCt                                                             <  2:311639/98‑1 (MQ=255)
                                                                                                                                          gcgcTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCt                                                             >  1:311639/1‑98 (MQ=255)
                                                                                                                                           cgcTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACGCCGATGCGCTGGAAGATg        <  2:1055870/150‑1 (MQ=255)
                                                                                                                                           cgcTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACGCCGATGCGCTGGAAGATg        <  1:984831/150‑1 (MQ=255)
                                                                                                                                               aTCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCa                                                                                                                       <  2:603205/35‑1 (MQ=255)
                                                                                                                                               aTCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCa                                                                                                                       >  1:603205/1‑35 (MQ=255)
                                                                                                                                                 cgcgGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACGCCGATGCGCTGGAAGATGAAGTCg  <  1:1131535/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TACGCGCCGGTGGACACGCCCAGCCTGAGCGAAGCCCAGGTCAAGTCGCGCTGGGAAGAAATCAAGGCCGAAGAGCGCGCGGAGAAGAGCCAGCCCGAACAGCTGTCGCTGCTTGATGACGTGCCGGCAACCTTGCCGGCGCTATCGCGGGCAGCCAAACTGCAAAAACGTGCGGCCACGGTCGGTTTCGACTGGCCGGCGGCGTTGCCGGTGCTGGACAAGGTCCGCGAAGAGCTGGATGAAGTGCTGCAGGCCATGGCCGACGGTGACGCCGATGCGCTGGAAGATGAAGTCG  >  NC_002947/1853639‑1853933

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