Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,901 A→C 28.5% R661R (CGA→CGC rsxC → SoxR [2Fe‑2S] reducing system protein RsxC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,9010AC28.5% 68.4 / 19.0 63R661R (CGA→CGCrsxCSoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base A (28/17);  new base C (16/2);  total (44/19)
Fisher's exact test for biased strand distribution p-value = 6.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGA  >  NC_000913/1708763‑1709047
                                                                                                                                          |                                                                                                                                                  
ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAGGCTCGGCAGCACGGAGGTCCGCGCAAAGCGGCCGGGGAAGGCGCGATTGCCCGCGGCACAgcgca                                                                                                                                         >  1:153759/1‑149 (MQ=2)
        cgcgCAAACCCGCCCTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCATAACAACATTTCGATCCGCTCAAATCCGGCGTCGAATCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                                                                                                                 <  2:1228620/150‑1 (MQ=17)
             aaaCCCGCCGTCGAAGCAGCTATTGCCCGGGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAg                                                                                                                            <  1:1760614/150‑1 (MQ=255)
                gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCCGCAGTAGCAGCGGCTATTGCCCGAGCGAAGGCCAGCAACCGGGAACAGCag                                                                                                                         >  2:1525804/1‑149 (MQ=11)
                        cGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCCACGCCCGCAAACGGGACCAGCACCCGGCTaa                                                                                                                 >  2:919432/1‑150 (MQ=17)
                             ccGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGACCAGCACCCGGCTAATGCgg                                                                                                            >  1:328675/1‑150 (MQ=255)
                                 tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCGAACGGGAGCAGCCACCGGCTAATGCGGGGcc                                                                                                        >  1:68091/1‑150 (MQ=255)
                                     gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTATACGCCGCTATTGCCACAGCCAAACCGCGCACACCGGACCAGAAACCCGCTAATGCCGAGACagaa                                                                                                    >  1:1120080/1‑150 (MQ=14)
                                          tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCCGTAGACGCGGCTATTGCCCGCGCGAAAGCCCGCAAACCCGCCCAGCAGCAGGATAATGCGGAGCCAGAAGAGCa                                                                                               >  2:367006/1‑150 (MQ=14)
                                           gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTAGCCGCCGCTATTGCCCGCGCAAAAGCGCGCAAACGGGAACAGCACCAGGATAATGCGGAGCAAGAAGACCAg                                                                                              >  2:1790340/1‑150 (MQ=14)
                                             cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGCCGCCGCTATTGCCCGCGCCAAAGCCAGCAAACGCGAACAGCAACAGGCTGATACGGAGCCAGAAGAACACGt                                                                                            >  1:1589260/1‑150 (MQ=14)
                                                   tcgaaaGCTTTAACAGCAACAGGCTAATGCGGATCAAGAACTACATGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc                                                                                      <  1:1927985/146‑1 (MQ=14)
                                                   tcGCAAGCTTTAACAGCAACAGTCTAATGCGTATCAAGAACAACATGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc                                                                                      <  1:1937414/149‑1 (MQ=14)
                                                          cTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGAGCAAACCCGAACAGCACCAGGGTGAAGCCGAGCCAGAAGAACAGGTGGATACGCGaaaa                                                                               >  1:2239114/1‑150 (MQ=11)
                                                               aCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCCAAAGCCCGCAAACGGGACCAGCAACCGGCTAAAGCGGAGCCAGAAGAACACGTAGATACGCACAAagcagc                                                                          >  1:1668055/1‑150 (MQ=14)
                                                                cAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCCAAAGCCCGCAAACGCGAACAGCACCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGca                                                                         >  2:1942519/1‑149 (MQ=17)
                                                                  gCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCGGCTATTGCCCGCGCGAAAGCCAGCAAACGGGAACAGCAACAGGCTAATGAGGAGCCAGAAGAACAGGTTGATCAGCGCAAAGCGGCAGt                                                                       >  2:1890714/1‑150 (MQ=14)
                                                                   cAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCAGCAAACGGGAACAGCAACAGGCTAATGAGGAGCCAGAAGAACAGGTTCGTACGCGCAACGCGGCAGtt                                                                      >  2:1587928/1‑149 (MQ=14)
                                                                     aCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                                                                                                                >  2:1873305/1‑90 (MQ=255)
                                                                     aCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                                                                                                                <  1:1873305/90‑1 (MQ=255)
                                                                      cAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTCCCAGAGCCAAAGAGAGCAAACGGGAACAGAAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATAAGCGCAAAGCTGCCGTCGac                                                                   >  2:593129/1‑149 (MQ=17)
                                                                       aGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGtagcac                                                                  >  1:1981059/1‑145 (MQ=255)
                                                                        ggCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCGGCTATTGCCCGCGCAAAAGCCCGCAAACGCGAACAGCAACAGGCTAATGCGGAGCAATAAGAACAGTTAGATCCGCGCAAAGCAGCAGttgccgc                                                                 >  2:2288682/1‑144 (MQ=11)
                                                                          cTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTAGAAGCCGCTATTGCCCGCGCAAAAGCCAAAAAACACGAACAGCAACACGCTAAAGCGGAGACAGAAGAACAAGTTGATACACGAAAAGCTGACGTTGCAGccg                                                               >  2:1696469/1‑148 (MQ=9)
                                                                           cactGCCGAGCCAGAACAACAGGTCGATCCGCGCAACGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGc                                                              <  1:593129/147‑1 (MQ=255)
                                                                              tGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTAt                                                           >  1:1697181/1‑150 (MQ=255)
                                                                               gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGCGCCAAAGCCCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGACCAGGTGGATCCGCGCAAAGCGGCCGGCGCCGCGGCTagt                                                          >  1:299918/1‑148 (MQ=14)
                                                                               gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAACCGGGAACAGCAACCGGCTAATGCGGGGCCAGAAGAACAGGTTGATCCGCGCATAGCGGCCGTGGCAGCGGCTAtt                                                          >  1:2373825/1‑150 (MQ=14)
                                                                               gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCCCGCAAACGGGAACAGCACCCGGGTAATGCGGAGCCAGAAGAACCGGTTGATCCGCGCAAAGCGGCAGTGGAAGCGGCTAtt                                                          >  1:1989277/1‑150 (MQ=17)
                                                                                ccGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTAGATACGAGCAAAGCTGCCGTCTAAACGGCTATTg                                                         >  2:1767809/1‑150 (MQ=17)
                                                                                 ggaaccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                                                                                                   >  1:1352380/5‑107 (MQ=255)
                                                                                 ggaaccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                                                                                                   <  2:1352380/103‑1 (MQ=255)
                                                                                  gAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACACGTTGATCCGCGCAAAGCTGCAGTAGCCGCGGCTAttccc                                                       >  1:1664811/1‑147 (MQ=17)
                                                                                    gCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCAcg                                                     >  1:23000/1‑150 (MQ=255)
                                                                                       aGAACAACACGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCCGCTAATGCCGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCTGACGTCGACGCGCCTATTGCACgcgc                                                  >  1:638107/1‑150 (MQ=17)
                                                                                         aacaacAGGTCTATCCTCGCAAAGCCGCCGTCGAATCCGATATTTCCCGAGCCAAAGCGCGCAAACGTTAACTTCAACCTGCTAATGTGGATCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCa                                                <  2:1205165/150‑1 (MQ=21)
                                                                                            aacaGCTCGATCCGGGCCATGCCTCGGTCGAAGCCGCTATCTCCCGAGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAg                                             <  1:922303/150‑1 (MQ=21)
                                                                                             acaGGTCTATCCGCCCAACCCCCCCGCAGAAGCCCCTACTGCGCGAGCCCATTCACGCCTGCGGTCACCCGAACCGGCTAATGCGGACCCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGc                                            <  1:287666/150‑1 (MQ=9)
                                                                                             acaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                                                                                  <  1:185231/96‑1 (MQ=255)
                                                                                             acaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                                                                                  >  2:185231/1‑96 (MQ=255)
                                                                                             acaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGCACGCGCCAAAGc                                            >  2:35442/1‑150 (MQ=255)
                                                                                             acaGGTCGATCCCCGCACAGCCTCCTGCGACCCCGCTATTGCCCGCGCCCAAGCGCGCAACCGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGc                                            <  1:2282650/150‑1 (MQ=14)
                                                                                               aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAg                                                                                                          <  2:1003247/86‑1 (MQ=255)
                                                                                               aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAg                                                                                                          >  1:1003247/1‑86 (MQ=255)
                                                                                                ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                                                                                                                >  2:734185/1‑63 (MQ=255)
                                                                                                ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                                                                                                                <  1:734185/63‑1 (MQ=255)
                                                                                                ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGACCAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCAGCGCAAAGCTGCCGTAGCAGCGGCTATTGCACGCGCCAAAGccag                                         >  1:983407/1‑147 (MQ=17)
                                                                                                cgTCGATCCGCGCAAAGCCGCCGTCTAAGACGCTATTGCCCGAGCCAAAGCTCGCAATCTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACg                                         <  2:68091/149‑1 (MQ=14)
                                                                                                     aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGATGCCCTGTCAGCGGCTATTCCACGGGGCGAAGCACGCaacc                                    >  1:244322/1‑148 (MQ=17)
                                                                                                           gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGCCGCGGCTATTGCCCGCGCCCACGCCCGCAAggaggacc                              >  1:2456770/1‑143 (MQ=17)
                                                                                                                gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCaa                         >  2:350939/1‑150 (MQ=255)
                                                                                                                       tCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCg                                                                        >  1:590301/1‑96 (MQ=255)
                                                                                                                       tCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCg                                                                        <  2:590301/96‑1 (MQ=255)
                                                                                                                           aGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCACGCGCCAACGCCCGCAAGCTGGCAGAGCAACAGGCTACTGc              >  1:2477703/1‑150 (MQ=25)
                                                                                                                           aGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCAAGCAAGCAGGAACAGCAACAGGCTAATGc              >  1:1401561/1‑150 (MQ=255)
                                                                                                                           aGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCCCGCGCAAAAGCCCAAAAACACGAACAGCAAAAGGATaaaaa              >  2:1053626/1‑147 (MQ=18)
                                                                                                                             ccGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTCGAACAGCAACAGGCTAATGCgg            >  1:557230/1‑150 (MQ=255)
                                                                                                                              cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGt           >  2:926815/1‑150 (MQ=255)
                                                                                                                                   ttGCCCGATCCAAAGCGCGCAAACTGTAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCag      <  1:1767809/150‑1 (MQ=255)
                                                                                                                                     gCCCGCGCCAAAGCGCGCAATCTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATc                                                                                       <  1:2177913/67‑1 (MQ=14)
                                                                                                                                     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCa                                       >  1:1827850/1‑115 (MQ=255)
                                                                                                                                     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCa                                       <  2:1827850/115‑1 (MQ=255)
                                                                                                                                       ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGAACGCAACAAAGAACTCAAGCTGGAACAGCAAAAGGCTAATGCGGAACCagaaga  >  2:1789992/1‑150 (MQ=17)
                                                                                                                                          |                                                                                                                                                  
GGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGA  >  NC_000913/1708763‑1709047

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