Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,946,036 A→C 18.9% intergenic (‑60/+79) ruvA ← / ← yobI component of RuvABC resolvasome, regulatory subunit/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,946,0360AC18.9% 322.7 / 73.6 139intergenic (‑60/+79)ruvA/yobIcomponent of RuvABC resolvasome, regulatory subunit/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (65/46);  new base C (14/12);  total (81/58)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 152 total aligned reads displayed.

GTTGTTTTTCAATGATGATGCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAGCAAATATT  >  NC_000913/1945937‑1946134
                                                                                                   |                                                                                                  
gttgttTTTCAATGATGATGCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGgataga                                                                                                    <  2:16992/100‑1 (MQ=255)
  tgttTTTCAATGATGATGCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGgatagata                                                                                                  >  2:682519/1‑100 (MQ=255)
    tttttCAATGATGATGCCTCTGAGTCTGCCTATCACATAAAGCTACTGCGTAATTAATCAAAGATAAAAAAAAAAAAAAAAAAAATGCTGGATAGatat                                                                                                 <  2:2422378/99‑1 (MQ=255)
      tttCAATGATGATGCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCa                                                                                              >  2:2971012/1‑100 (MQ=255)
       ttCAATGATGATGCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAg                                                                                             >  2:805818/1‑100 (MQ=255)
        tCAATGATGATGCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGc                                                                                            <  2:3618640/100‑1 (MQ=255)
          aaTGATGATGCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATATCTATCCAGcat                                                                                          <  1:1396951/100‑3 (MQ=255)
          aaTGATGATGCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGa                                                                                          >  1:3074606/1‑100 (MQ=255)
           atgatgatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAGAAAATGCTGGATAGATATCCAGCGaa                                                                                         >  1:1577643/1‑100 (MQ=255)
            tgatgatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAg                                                                                        >  2:3359901/1‑100 (MQ=255)
            tgatgatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAg                                                                                        >  2:3358703/1‑100 (MQ=255)
            tgatgatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAg                                                                                        <  1:2207952/100‑1 (MQ=255)
            tgatgatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAg                                                                                        <  2:1614747/100‑1 (MQ=255)
             gatgatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAgg                                                                                       <  1:2740968/100‑1 (MQ=255)
              atgatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATATCTATCCAGcattttt                                                                                      <  2:720928/100‑7 (MQ=255)
              atgatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGa                                                                                      >  1:782683/1‑100 (MQ=255)
                gatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATATCTATCCAGcattttttt                                                                                    <  2:3588564/100‑9 (MQ=38)
                gatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATATCTATCCAGcattttttt                                                                                    <  1:238075/100‑9 (MQ=38)
                gatgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATATCTATCCAGcattttttt                                                                                    <  1:3103132/100‑9 (MQ=38)
                 atgCCTCTGAGTCTACCTATCACATGACGCTACTGCGTAATGAAGCAAAGATAATGCTGTATGATAAAAAAATGCTGGATATCTATCCAGcattttttt                                                                                    <  2:811786/99‑9 (MQ=37)
                  tgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATATCTATCCAGcatttttttat                                                                                  <  1:1031263/100‑11 (MQ=38)
                  tgCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATATCTATCCAGcatttttttat                                                                                  <  1:3143614/100‑11 (MQ=38)
                   gCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATgaag                                                                                 <  1:3246131/100‑1 (MQ=255)
                    cttccgatcttGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATgaaga                                                                                <  1:3451201/90‑1 (MQ=255)
                    ccTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATgaaga                                                                                >  2:430001/1‑100 (MQ=255)
                    ccTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATgaaga                                                                                >  1:4132803/1‑100 (MQ=255)
                    ccTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATgaaga                                                                                <  2:866076/100‑1 (MQ=255)
                    ccTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATgaaga                                                                                <  1:3321454/100‑1 (MQ=255)
                     ctctGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAACGATAATGCTGTATGATCAAAAAATGCTGGATAGATATCCAGCGAAGGATgaagaa                                                                               >  1:3269169/1‑100 (MQ=255)
                     ctctGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGCTgaagaa                                                                               >  2:2743212/1‑100 (MQ=255)
                       ctGAGTCTGCCTATCACATGACGCGCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGTTAGATATCCAGCGAAGGATGAAGaaaa                                                                             <  2:480772/100‑1 (MQ=255)
                        tGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAAc                                                                            >  2:300867/1‑100 (MQ=255)
                            tCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGc                                                                        <  2:959331/100‑1 (MQ=255)
                             cTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGTTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCg                                                                       >  1:1022830/1‑100 (MQ=255)
                             cTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATATATATCCAGCGAAGGATGAAGAAAACTTGCg                                                                       <  1:1380329/100‑1 (MQ=255)
                             cTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCg                                                                       <  2:115935/100‑1 (MQ=255)
                              tGCCTATGACATGACGCTCCTGCGGAATGAATCAAACATAATGCTGTAGGAGAAAAAAATGCTGGAGAGGTAGCCAGCGAAGGATGAAGAGAACTTGCGa                                                                      >  2:2442180/1‑100 (MQ=255)
                              tGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATgaagaagatcggaag                                                                      >  2:3451185/1‑91 (MQ=255)
                               gCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAg                                                                     >  1:3142144/1‑100 (MQ=255)
                                 cTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGgt                                                                   <  1:2607676/100‑1 (MQ=255)
                                 cTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGgt                                                                   >  1:424475/1‑100 (MQ=255)
                                  tATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGgtg                                                                  >  1:121530/1‑100 (MQ=255)
                                  tATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGgtg                                                                  <  1:3050640/100‑1 (MQ=255)
                                  tATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGgtg                                                                  >  1:2074761/1‑100 (MQ=255)
                                  tATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGgtg                                                                  <  1:478609/100‑1 (MQ=255)
                                   aTCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGgtgt                                                                 >  1:338161/1‑100 (MQ=255)
                                   aTCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGgtgt                                                                 >  1:1310677/1‑100 (MQ=255)
                                   aTCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGgtgt                                                                 >  2:1836932/1‑100 (MQ=255)
                                     cacaTGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGtct                                                               >  1:1091904/1‑100 (MQ=255)
                                      acaTGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGtct                                                               >  2:312881/1‑99 (MQ=255)
                                       caTGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCg                                                             <  2:2255700/100‑1 (MQ=255)
                                         tGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGAt                                                           >  1:5917/1‑100 (MQ=255)
                                         tGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGAt                                                           <  1:3634157/100‑1 (MQ=255)
                                              cTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACt                                                      >  2:1435392/1‑100 (MQ=255)
                                              cTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACt                                                      <  1:3419240/100‑1 (MQ=255)
                                               tCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTg                                                     <  1:3694245/100‑1 (MQ=255)
                                                ccTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGa                                                    >  1:2071930/1‑100 (MQ=255)
                                                 cTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGaa                                                   >  2:789840/1‑100 (MQ=255)
                                                  tGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGaaa                                                  >  1:364075/1‑100 (MQ=255)
                                                   gCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGaaaa                                                 >  2:60730/1‑100 (MQ=255)
                                                      taataaATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAAGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATg                                              <  2:1233192/100‑1 (MQ=255)
                                                       aaTGAAATAAAGATAATGCTGTAGGAGAAAAAAATGCTGGATAGATAACCAGCGAAGAATGAAGAAAACTTGCGAGGTGTCAAAATAAAAAGAAAAATg                                              <  2:1443824/99‑1 (MQ=255)
                                                          gAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGcac                                          <  1:1956048/100‑1 (MQ=255)
                                                          gAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGcac                                          >  1:3820413/1‑100 (MQ=255)
                                                             tCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGt                                       >  1:1270984/1‑100 (MQ=255)
                                                                aaGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGtata                                    >  2:1358256/1‑100 (MQ=255)
                                                                    tAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAAttt                                <  2:3879668/100‑1 (MQ=255)
                                                                          tGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTa                          >  2:1942784/1‑100 (MQ=255)
                                                                          tGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTa                          >  2:3281455/1‑100 (MQ=255)
                                                                          tGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTa                          <  1:1514726/100‑1 (MQ=255)
                                                                           gTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTAc                         >  2:2344857/1‑100 (MQ=255)
                                                                           gTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTAc                         >  2:2485623/1‑100 (MQ=255)
                                                                             aTGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAg                       <  1:2386446/100‑1 (MQ=255)
                                                                             aTGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAg                       <  1:2807310/100‑1 (MQ=255)
                                                                             aTGATAAAAAAATGCTGGATAGATATCCAGCAAAAGATGAAGAAAACTTGCCAGGTGTCTAAATAAAAAAAAAAATGGCACAGTATAATTTATTCTACa                        <  2:3764982/99‑1 (MQ=255)
                                                                              tGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAgg                      >  2:3439526/1‑100 (MQ=255)
                                                                                aTAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGtt                    >  1:2293764/1‑100 (MQ=255)
                                                                                aTAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGtt                    >  1:1142355/1‑100 (MQ=255)
                                                                                 tAAAAAAATGGTGGATAGCTGGGCAGCGGAGGTTGAAGACAACTTGCGAGGTGTCTCGATGAACTGAAAAGTGACACAGTATAATTTATTAGACAGGAta                   >  2:419015/1‑100 (MQ=255)
                                                                                  aaaaaaaTGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTtat                  >  2:3088608/1‑100 (MQ=255)
                                                                                   atccttcgCTGGATATCTATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTtata                 >  1:728875/8‑100 (MQ=255)
                                                                                   aaaaaaTGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGCACTGAAACATGGCACAGTATAATTTATTCTACAGGTtata                 >  1:2548928/1‑100 (MQ=255)
                                                                                      cttcgCTGGATATCTATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTg              >  1:2628547/5‑100 (MQ=255)
                                                                                       ttcgCTGAATATTTATCCAGCGAAAGATGAAGAAATCTTATAATTTGTATCGATTAATTAAAAAATGGCACAATATAATTTATTCTACAAATTATATTg              >  2:1644164/4‑99 (MQ=255)
                                                                                        tcgCTGGATATCTATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTg              >  2:2511332/3‑98 (MQ=255)
                                                                                         tGCTGGATAGATATCCAGCGAAGGATGAAGAAAACCTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGaa           <  1:381931/100‑1 (MQ=255)
                                                                                         cgCTGGATATCTATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGaa           >  1:65918/2‑100 (MQ=255)
                                                                                         cgCTGGATATCTATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGaa           >  1:473592/2‑100 (MQ=255)
                                                                                         cgCTGGATATCTATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGCCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGaa           >  2:3092416/2‑100 (MQ=255)
                                                                                          gCTGGATATCTATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAg          >  2:1259467/1‑100 (MQ=255)
                                                                                          gCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAg          <  1:3682973/100‑1 (MQ=255)
                                                                                           cTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAGc         >  1:3048020/1‑100 (MQ=255)
                                                                                           cTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAGc         >  1:1774513/1‑100 (MQ=255)
                                                                                             ggATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAGCaa       >  1:1368717/1‑100 (MQ=255)
                                                                                             ggATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAGCaa       >  1:971954/1‑100 (MQ=255)
                                                                                             ggATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAGCaa       >  2:2457966/1‑100 (MQ=255)
                                                                                                tagataTCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAGCAAata    >  1:840916/1‑100 (MQ=255)
                                                                                                nagataTCCAGCGAAGGATGAGGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACGGTATAATTTATTCTACAGGTGATATTGGAAGCAAata    >  2:1010712/2‑100 (MQ=255)
                                                                                                  gataTCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAGCAAATAtt  >  1:339903/1‑100 (MQ=255)
                                                                                                  gataTCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAGCAAATAtt  >  1:339835/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GTTGTTTTTCAATGATGATGCCTCTGAGTCTGCCTATCACATGACGCTCCTGCGTAATGAATCAAAGATAATGCTGTATGATAAAAAAATGCTGGATAGATATCCAGCGAAGGATGAAGAAAACTTGCGAGGTGTCTCGATGAACTGAAAAATGGCACAGTATAATTTATTCTACAGGTTATATTGGAAGCAAATATT  >  NC_000913/1945937‑1946134

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