Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,506,939 T→C 20.3% intergenic (+51/‑66) petC → / → sspA ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,506,9390TC20.3% 125.7 / 26.3 64intergenic (+51/‑66)petC/sspAubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A
Reads supporting (aligned to +/- strand):  ref base T (23/28);  new base C (6/7);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

AGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGCCGAGAAG  >  NC_002947/1506793‑1507088
                                                                                                                                                  |                                                                                                                                                     
aGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGc                                                                                                                                                    <  2:1709452/150‑1 (MQ=255)
   cAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg                                                                                                                                                 >  2:798683/1‑150 (MQ=255)
   cAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg                                                                                                                                                 <  1:1345097/150‑1 (MQ=255)
   cAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg                                                                                                                                                 <  1:1146565/150‑1 (MQ=255)
     gcgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGc                                                                                                                                               >  1:286784/1‑150 (MQ=255)
       gcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt                                                                                                                                                      >  1:66171/1‑141 (MQ=255)
       gcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt                                                                                                                                                      <  2:66171/141‑1 (MQ=255)
               aCCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGgc                                                                                                                                     >  1:576409/1‑150 (MQ=255)
                   aTGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGt                                                                                                                                 <  2:576409/150‑1 (MQ=255)
                      tgtTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCCGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGttca                                                                                                                              <  1:1106497/150‑3 (MQ=38)
                              tATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTGATTTaaa                                                                                                                      >  1:1438035/1‑150 (MQ=255)
                                  tGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGTCGCGTTTTCATTTAAAGTTg                                                                                                                  >  2:248532/1‑150 (MQ=255)
                                   ggCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGc                                                                                                                 <  2:493263/150‑1 (MQ=255)
                                    gCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGca                                                                                                                >  2:1357374/1‑150 (MQ=255)
                                       ttcttcttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCGCCCCGGGCGCGTTTTCATTTAAAGTTGCACaa                                                                                                             <  1:1356256/150‑1 (MQ=255)
                                        tcttcttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAg                                                                                                            <  1:1357374/150‑1 (MQ=255)
                                                  aTTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGc                                                                                                                                                    <  2:738029/100‑1 (MQ=38)
                                                   ttCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCga                                                                                                 <  1:1703058/150‑1 (MQ=255)
                                                   ttCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCga                                                                                                 <  1:146625/150‑1 (MQ=255)
                                                      gCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgagga                                                                                              >  1:747664/1‑150 (MQ=255)
                                                          aCTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCg                                                                                          <  2:1291229/150‑1 (MQ=255)
                                                           cTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTaaa                                                                                                                      <  2:1696870/121‑1 (MQ=255)
                                                           cTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTaaa                                                                                                                      >  1:1696870/1‑121 (MQ=255)
                                                            ttGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGCCACCCCg                                                                                                                                        >  1:1116196/1‑102 (MQ=25)
                                                            ttGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGCCACCCCg                                                                                                                                        <  2:1116196/102‑1 (MQ=25)
                                                            ttGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCGGTGGGAGGAGCCGCt                                                                                        >  1:63653/1‑150 (MQ=255)
                                                               cTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACa                                                                                     <  2:1311245/150‑1 (MQ=255)
                                                                tCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGcc                                                                                                                                                   >  2:686273/1‑87 (MQ=255)
                                                                tCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGcc                                                                                                                                                   <  1:686273/87‑1 (MQ=255)
                                                                tCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGc                                                                                                           >  2:1350430/1‑127 (MQ=255)
                                                                tCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGc                                                                                                           <  1:1350430/127‑1 (MQ=255)
                                                                tCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACAt                                                                                    >  1:1050563/1‑150 (MQ=255)
                                                                  aaGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGGGAGGAGCCGCTACATgg                                                                                  <  2:1050563/150‑1 (MQ=255)
                                                                  aaGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATgg                                                                                  <  1:1597267/150‑1 (MQ=255)
                                                                       tGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGACCGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGca                                                                             <  1:426877/150‑1 (MQ=255)
                                                                        gAATACTGGAAGGACGTGCACTGATGACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATAGCGAGGGCGCCCTGTGCGCGTTTTCATTTAAAGCTGCACAATCATCCTTGCGAGTAGCCGCTACATGGGCAcaa                                                                            <  2:520502/150‑1 (MQ=255)
                                                                         aaTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaac                                                                           <  2:747664/150‑1 (MQ=255)
                                                                             ggggAAGGGCGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGCACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCGCCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGCGCCGCTACATGGGCGccaccaac                                                                       <  2:286784/148‑1 (MQ=255)
                                                                              tGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACa                                                                      >  2:1608155/1‑150 (MQ=255)
                                                                                       gTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCt                                                             >  2:1111082/1‑150 (MQ=255)
                                                                                         gCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCTCCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGc                                                           <  2:356878/150‑1 (MQ=255)
                                                                                                 cACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc                                                                                                                                               >  1:584731/1‑58 (MQ=14)
                                                                                                 cACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc                                                                                                                                               <  2:584731/58‑1 (MQ=14)
                                                                                                  aCGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGtt                                                                                                                                >  1:1395749/1‑72 (MQ=14)
                                                                                                  aCGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGtt                                                                                                                                <  2:1395749/72‑1 (MQ=14)
                                                                                                     cAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCa                                                                                                                            >  2:1554610/1‑73 (MQ=255)
                                                                                                     cAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCa                                                                                                                            <  1:1554610/73‑1 (MQ=255)
                                                                                                      agtagtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCg                                              >  2:1720492/1‑150 (MQ=255)
                                                                                                          gtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGa                                          <  1:1111082/150‑1 (MQ=255)
                                                                                                                 gttgttGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTCTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTAt                                   <  1:798683/150‑1 (MQ=255)
                                                                                                                           gcgcCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATgg                                                                                  <  1:1240661/93‑1 (MQ=255)
                                                                                                                           gcgcCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATgg                                                                                  >  2:1240661/1‑93 (MQ=255)
                                                                                                                              cccGGGGTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATccc                                               <  2:728938/125‑1 (MQ=37)
                                                                                                                              cccGGGGTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATccc                                               >  1:728938/1‑125 (MQ=38)
                                                                                                                              cccGGGGGGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTc                      >  2:553822/1‑150 (MQ=255)
                                                                                                                                    gTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTgg                >  1:1469848/1‑150 (MQ=39)
                                                                                                                                    gTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTAtt                                                       >  1:1334312/1‑111 (MQ=255)
                                                                                                                                    gTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTAtc                                                       <  2:1334312/111‑2 (MQ=255)
                                                                                                                                      ggCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATc                            <  1:1080517/136‑1 (MQ=39)
                                                                                                                                      ggCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATc                            >  2:1080517/1‑136 (MQ=39)
                                                                                                                                       gCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTc             >  2:16154/1‑150 (MQ=255)
                                                                                                                                         cTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTg           >  2:205149/1‑150 (MQ=255)
                                                                                                                                              aGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGCCga      >  2:1183603/1‑150 (MQ=255)
                                                                                                                                                tATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATc                            >  1:139551/1‑126 (MQ=255)
                                                                                                                                                tATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATc                            <  2:139551/126‑1 (MQ=255)
                                                                                                                                                 aTTGCCAGGCCTCCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGCCGAGaa   <  1:248532/150‑1 (MQ=255)
                                                                                                                                                  ttGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGCCGAGAAg  <  1:1183603/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGCCGAGAAG  >  NC_002947/1506793‑1507088

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