Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,628,638 C→A 15.4% A96E (GCG→GAG)  rspR → transcriptional repressor for rspAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,628,6380CA15.4% 304.1 / 37.4 108A96E (GCG→GAG) rspRtranscriptional repressor for rspAB
Reads supporting (aligned to +/- strand):  ref base C (46/45);  new base A (6/11);  total (52/56)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 110 total aligned reads displayed.

CCTGATTCAAATTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCGCAAGCAAC  >  NC_000913/1628540‑1628736
                                                                                                  |                                                                                                  
ccTGATTCAAATTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGcgcg                                                                                                   <  2:522000/100‑1 (MQ=255)
 ntGATTCAAATTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCgg                                                                                                  >  1:1063582/2‑100 (MQ=255)
 gtGATTCAAATTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCgg                                                                                                  >  1:1470695/2‑100 (MQ=255)
  tGATTCAAATTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGt                                                                                                 <  2:2631871/100‑1 (MQ=255)
      tCAAATTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGgcg                                                                                             >  1:767142/1‑100 (MQ=255)
      tCAAATTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGgcg                                                                                             <  2:1561308/100‑1 (MQ=255)
       cAAATTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGgcgc                                                                                            >  2:3018669/1‑100 (MQ=255)
        aaaTTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGcgcg                                                                                           <  2:2153394/100‑1 (MQ=255)
          aTTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTc                                                                                         <  2:1930779/100‑1 (MQ=255)
           ttCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGAGGTGGCGCGTCg                                                                                        >  1:1738686/1‑100 (MQ=255)
              gTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGAGGTGGCGCGTCGggc                                                                                     <  1:344202/100‑1 (MQ=255)
               tCCGCAACGTAGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGggcg                                                                                    >  1:12159/1‑100 (MQ=255)
                 cGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGggcggc                                                                                  >  2:2566922/1‑100 (MQ=255)
                 cGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGggcggc                                                                                  >  2:1608661/1‑100 (MQ=255)
                  gCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGgcggcg                                                                                 <  2:10965/100‑1 (MQ=255)
                  gCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGgcggcg                                                                                 <  2:60281/100‑1 (MQ=255)
                  gCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGAGGTGGCGCGTCGGgcggcg                                                                                 <  1:147986/100‑1 (MQ=255)
                   cAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGa                                                                                <  1:169659/100‑1 (MQ=255)
                     aCGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGc                                                                              <  2:282343/100‑1 (MQ=255)
                     aCGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGAGGTGGCGCGTCGGGCGGCGAGc                                                                              <  1:3116191/100‑1 (MQ=255)
                      cGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCa                                                                             <  2:1024486/100‑1 (MQ=255)
                      cGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCa                                                                             >  1:119706/1‑100 (MQ=255)
                       gTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGAGGTGGCGCGTCGGGCGGCGAGCAt                                                                            >  1:1225568/1‑100 (MQ=255)
                        tGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATg                                                                           <  1:2457794/100‑1 (MQ=255)
                         ggCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCATCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGa                                                                          <  2:1832930/100‑1 (MQ=255)
                             gCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTAc                                                                      <  2:240300/100‑1 (MQ=255)
                              cTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTAcc                                                                     <  2:1777077/100‑1 (MQ=255)
                               tatgTCAACTAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTAcc                                                                     <  1:743673/96‑1 (MQ=255)
                               tACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCg                                                                    >  2:1539579/1‑100 (MQ=255)
                               tACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTACCg                                                                    >  2:261309/1‑100 (MQ=255)
                                aCGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGa                                                                   <  1:1130613/100‑1 (MQ=255)
                                aCGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGa                                                                   >  2:796712/1‑100 (MQ=255)
                                 cGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGaa                                                                  >  1:2916380/1‑100 (MQ=255)
                                  gTCAACAAAATTTCCATGGCCCAGGTGCGGAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGaaa                                                                 <  1:320835/100‑1 (MQ=255)
                                    caacaaAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGc                                                               >  2:885424/1‑100 (MQ=255)
                                     aacaaAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGcc                                                              <  1:2849113/100‑1 (MQ=255)
                                     aacaaAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGcc                                                              >  2:21795/1‑100 (MQ=255)
                                     aacaaAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGcc                                                              >  2:682418/1‑100 (MQ=255)
                                      acaaAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCa                                                             >  2:2722940/1‑100 (MQ=255)
                                      acaaAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCa                                                             >  2:651415/1‑100 (MQ=255)
                                      acaaAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCa                                                             >  2:2946765/1‑100 (MQ=255)
                                         aaaTTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTg                                                          <  2:3320266/100‑1 (MQ=255)
                                         aaaTTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTg                                                          >  1:1310010/1‑100 (MQ=255)
                                          aaTTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGc                                                         >  1:3304514/1‑100 (MQ=255)
                                             ttCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTAt                                                      >  1:3160944/1‑100 (MQ=255)
                                             ttCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTAt                                                      <  2:407033/100‑1 (MQ=255)
                                              tCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATc                                                     <  1:2275337/100‑1 (MQ=255)
                                              tCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATc                                                     >  2:3089790/1‑100 (MQ=255)
                                                 aTGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCTATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAAc                                                  >  1:2243228/1‑100 (MQ=255)
                                                 aTGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAAc                                                  <  2:1924765/100‑1 (MQ=255)
                                                  tGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACt                                                 >  2:764031/1‑100 (MQ=255)
                                                   ggCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTAGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTg                                                <  1:282739/100‑1 (MQ=255)
                                                    gCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTgg                                               <  1:886626/100‑1 (MQ=255)
                                                    gCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTgg                                               >  1:3380109/1‑100 (MQ=255)
                                                    gCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTgg                                               <  2:3380138/100‑1 (MQ=255)
                                                     cccAGGTGGGCAAGGGCAGTGTTATCGTTGAGGCCATGGTGTGCGCGGTGGCGCCTCGGGCTGCGAGCATGATTACCAAAGGCGATATCTATCAAcaga                                               <  1:1212474/99‑4 (MQ=255)
                                                     cccAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGa                                              >  2:1950520/1‑100 (MQ=255)
                                                       cAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAAc                                            >  1:3071098/1‑100 (MQ=255)
                                                       cAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAAc                                            >  1:2507442/1‑100 (MQ=255)
                                                           tGCGCAACGGCAGTTTTATCCGTCAGGCCATCGAGTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACaaaa                                        <  1:1281070/100‑1 (MQ=255)
                                                              gcAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCt                                     <  1:2336641/100‑1 (MQ=255)
                                                              gcAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCt                                     <  2:2550618/100‑1 (MQ=255)
                                                              gcAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCt                                     >  2:697065/1‑100 (MQ=255)
                                                               cAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCtt                                    >  2:1306292/1‑100 (MQ=255)
                                                                 aCGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTcn                                  <  1:1245981/100‑2 (MQ=255)
                                                                 aCGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCa                                  >  1:1413159/1‑100 (MQ=255)
                                                                  cGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCAc                                 <  1:923731/100‑1 (MQ=255)
                                                                  cGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCAc                                 >  2:2354482/1‑100 (MQ=255)
                                                                  cGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCAc                                 >  2:1282269/1‑100 (MQ=255)
                                                                   ggCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCAcc                                <  1:1608654/100‑1 (MQ=255)
                                                                   ggCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCAcc                                <  1:1861740/100‑1 (MQ=255)
                                                                      aGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGc                             <  2:2142786/100‑1 (MQ=255)
                                                                      aGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGc                             <  2:1048658/100‑1 (MQ=255)
                                                                       gTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCa                            <  2:2811856/100‑1 (MQ=255)
                                                                       gTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCa                            >  2:75303/1‑100 (MQ=255)
                                                                       gTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCa                            >  2:960029/1‑100 (MQ=255)
                                                                       gTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCa                            >  2:2560746/1‑100 (MQ=255)
                                                                           tATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGc                        <  2:1654141/100‑1 (MQ=255)
                                                                           tATCCGTCAGGCCATTGAGTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGc                        <  1:2432820/100‑1 (MQ=255)
                                                                               cGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGTAAGCCAGTGCTATCAACTGGAACAAAATATTCACAAGCAACGCATTg                    >  1:570048/1‑100 (MQ=255)
                                                                                gTCAGGCCATTGAGTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGc                   <  1:2631200/100‑1 (MQ=255)
                                                                                  cAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCa                 <  1:1865436/100‑1 (MQ=255)
                                                                                    ggCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCAtt               >  1:2797442/1‑100 (MQ=255)
                                                                                     gCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTg              <  2:3101295/100‑1 (MQ=255)
                                                                                      ccATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGa             >  2:2106360/1‑100 (MQ=255)
                                                                                      ccATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGa             >  1:1209977/1‑100 (MQ=255)
                                                                                      ccATTGAGTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGa             <  2:2183546/100‑1 (MQ=255)
                                                                                       cATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAg            >  1:3026547/1‑100 (MQ=255)
                                                                                       cATTGAGTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAg            <  1:2036387/100‑1 (MQ=255)
                                                                                        aTTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAgc           <  2:3357088/100‑1 (MQ=255)
                                                                                         ttGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAgcg          <  2:1447946/100‑1 (MQ=255)
                                                                                           gAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgca        >  2:3232440/1‑100 (MQ=255)
                                                                                            aGTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgcaa       <  1:238122/100‑1 (MQ=255)
                                                                                             gTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgcaag      <  1:2946739/100‑1 (MQ=255)
                                                                                             gTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgcaag      <  1:2466481/100‑1 (MQ=255)
                                                                                             gTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgcaag      <  2:1516071/100‑1 (MQ=255)
                                                                                             gTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgcaag      >  1:1756258/1‑100 (MQ=255)
                                                                                             gTGCGAGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgcaag      >  2:454709/1‑100 (MQ=255)
                                                                                              tGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgcaagc     >  2:159260/1‑100 (MQ=255)
                                                                                                 gcgGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCGCAAGCAAc  >  2:1171998/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CCTGATTCAAATTCGTCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCGCAAGCAAC  >  NC_000913/1628540‑1628736

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: