Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 443,022 G→A 19.0% intergenic (+25/+29) thiI → / ← yajL tRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913443,0220GA19.0% 441.0 / 83.0 194intergenic (+25/+29)thiI/yajLtRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone
Reads supporting (aligned to +/- strand):  ref base G (83/74);  new base A (25/12);  total (108/86)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 200 total aligned reads displayed.

GGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCC  >  NC_000913/442923‑443121
                                                                                                   |                                                                                                   
ggggTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTgg                                                                                                     <  1:1297703/100‑1 (MQ=255)
 gggTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTgg                                                                                                     >  1:2172197/1‑99 (MQ=255)
 gggTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGt                                                                                                    >  1:2033344/1‑100 (MQ=255)
   gTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCa                                                                                                  <  1:3478960/100‑1 (MQ=255)
   gTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCa                                                                                                  <  2:4072384/100‑1 (MQ=255)
    tgatgaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAg                                                                                                 <  2:2198098/100‑1 (MQ=255)
      atgaGCCGTCTGCAGGCGCTCTATCTGCGCAAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGAt                                                                                               <  1:3660909/100‑1 (MQ=255)
        gaGCCGTCTGCAGGCGCTCTATCTGCGCGAGTAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGc                                                                                             <  1:2849003/100‑1 (MQ=255)
          gCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGa                                                                                           >  2:2293011/1‑100 (MQ=255)
            cGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGACg                                                                                         <  1:231465/100‑1 (MQ=255)
             gTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGACg                                                                                         <  1:908473/99‑1 (MQ=255)
               cTGCAGGCGCTCTATCTGCGCAAGCAGGGCCTTCACAATGTGAAGGTATAACGCCCGTAATTTGTGGTGTTTACGTCGCATCTAACCAGATACGACGt                                                                                        <  1:874741/98‑1 (MQ=255)
                tacAGGCGCTCTAACTGCGCGAGCAGGGCTTAAACAATGTGACGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTg                                                                                     >  2:2969042/3‑100 (MQ=255)
                tGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTg                                                                                     >  1:404451/1‑100 (MQ=255)
                 gCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGc                                                                                    >  1:3190238/1‑100 (MQ=255)
                 gCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGc                                                                                    >  1:186064/1‑100 (MQ=255)
                 gCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGCGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGc                                                                                    >  1:2961521/1‑100 (MQ=255)
                   aGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCg                                                                                  <  1:784320/100‑1 (MQ=255)
                   aGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCg                                                                                  <  2:4091040/100‑1 (MQ=255)
                   aGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCg                                                                                  <  1:100055/100‑1 (MQ=255)
                    ggCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGc                                                                                 >  2:927812/1‑100 (MQ=255)
                    ggCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGc                                                                                 >  2:3856370/1‑100 (MQ=255)
                      cgcTCTATCTGCGTGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTGTTTACGTCGCATCTGACCAGATGCGACGtaaaaaccac                                                                               <  1:459586/100‑10 (MQ=25)
                       gcTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATc                                                                              <  1:1136615/100‑1 (MQ=255)
                       gcTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCACATTTGGTCAGATACGACGTTTGCCGCATc                                                                              >  1:538116/1‑100 (MQ=255)
                         tctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCg                                                                            >  1:611521/1‑100 (MQ=255)
                            aTCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGaca                                                                         >  1:755382/1‑100 (MQ=255)
                              cTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACact                                                                       >  2:2940928/1‑100 (MQ=255)
                              cTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACact                                                                       >  2:216630/1‑100 (MQ=255)
                               tGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACacta                                                                      >  2:1521845/1‑100 (MQ=255)
                                    gAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGt                                                                 >  2:1254551/1‑100 (MQ=255)
                                      gCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaa                                                               >  2:3311235/1‑100 (MQ=255)
                                         gggCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaataa                                                            >  2:466531/1‑100 (MQ=255)
                                         gggCGTTAACAATGAGAAGGTATATCGCCCGTAATTTGTGGTTTTTCCGTCGCATGTCGGCAGATGCGAGGTTTACCGCATGCTACCCTACTGTtaata                                                             >  1:105423/1‑99 (MQ=255)
                                          ggCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTaataat                                                           <  1:3837396/100‑1 (MQ=255)
                                           gCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTaataat                                                           <  1:4162180/99‑1 (MQ=255)
                                           gCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAAtt                                                          >  1:1601760/1‑100 (MQ=255)
                                           gCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAAtt                                                          >  1:2215465/1‑100 (MQ=255)
                                            cTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCTACGTTTGCCGCATCCGACACTACTCGTAATAATta                                                         <  1:154916/100‑1 (MQ=255)
                                             tttAACAAGGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATtat                                                        <  1:1891886/100‑1 (MQ=255)
                                               tAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTaataatagatc                                                      <  2:616488/100‑6 (MQ=255)
                                               tAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATaa                                                      >  1:1630849/1‑100 (MQ=255)
                                                 ataaTGGGAAGGTATATCGCCCGGAATTTGTGGGTTTTACGGCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTATTAATTATaaa                                                     <  1:3098913/97‑1 (MQ=255)
                                                 aCAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAAt                                                    >  2:2350340/1‑100 (MQ=255)
                                                   aaTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATcc                                                  <  2:1654946/100‑1 (MQ=255)
                                                     tgtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCt                                                <  1:803669/100‑1 (MQ=255)
                                                     tgtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGAAGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCt                                                >  2:76216/1‑100 (MQ=255)
                                                      gtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTg                                               >  1:2472815/1‑100 (MQ=255)
                                                      gtgAAGGTATATCGCCCGTAATCTGTGGTTTCTACGTTGCATCTGGTTAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTg                                               <  1:2970126/100‑1 (MQ=255)
                                                           ggTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgccgcc                                          >  2:1865396/1‑100 (MQ=255)
                                                             tataTCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGCCACTACTCGTAATAATTATAAATCCCTGCCGCCAt                                        >  2:3261995/1‑100 (MQ=255)
                                                             tataTCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCAt                                        >  1:2616036/1‑100 (MQ=255)
                                                              atatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATc                                       >  1:1812847/1‑100 (MQ=255)
                                                              atatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATc                                       >  1:355971/1‑100 (MQ=255)
                                                              atatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATc                                       <  2:582920/100‑1 (MQ=255)
                                                               tatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCa                                      >  2:3760153/1‑100 (MQ=255)
                                                               tatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCa                                      <  1:1852180/100‑1 (MQ=255)
                                                                    cccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATCATAAATCCCTGCCGCCATCACCAGt                                 >  2:531370/1‑100 (MQ=255)
                                                                    cccGTAATCTGGGGTTTTTACGTCGCATCTGGTCAGATGCGAGGTGTCCCGGATCCGACACTACTCGTAATAAGTATAAATCCCTGGCGCCATCACCAg                                  <  2:4383116/99‑1 (MQ=255)
                                                                     ccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGtt                                >  1:577918/1‑100 (MQ=255)
                                                                     ccGTAATTTGTGGGTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGtt                                <  1:2058558/100‑1 (MQ=255)
                                                                          aTTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGAt                           >  1:3939509/1‑100 (MQ=255)
                                                                           tttGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATg                          <  2:393420/100‑1 (MQ=255)
                                                                           tttGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATg                          <  2:4225720/100‑1 (MQ=255)
                                                                             tgtgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                        <  2:3841614/100‑1 (MQ=255)
                                                                             tgtgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                        <  1:3311219/100‑1 (MQ=255)
                                                                             tgtgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                        <  1:424777/100‑1 (MQ=255)
                                                                             tgtgGTTTTTACGGCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                        <  2:1677445/100‑1 (MQ=255)
                                                                               tgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCAc                      <  2:185998/100‑1 (MQ=255)
                                                                                ggTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACt                     >  2:3224777/1‑100 (MQ=255)
                                                                                    caaaCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAt                 >  2:1906459/4‑100 (MQ=39)
                                                                                     ttACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATg                <  1:1683488/100‑1 (MQ=255)
                                                                                      aaCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATgg               >  2:2643837/2‑100 (MQ=255)
                                                                                      aaCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATgg               >  2:1028290/2‑100 (MQ=255)
                                                                                      aaCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATgg               >  2:2092619/2‑100 (MQ=255)
                                                                                       aCGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATggg              <  1:3239232/100‑1 (MQ=255)
                                                                                       aCGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATggg              <  2:4362233/100‑1 (MQ=255)
                                                                                       aCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCTCTTCATggg              >  1:358/1‑100 (MQ=39)
                                                                                       aCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATggg              >  1:2913049/1‑100 (MQ=255)
                                                                                       aCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATggg              >  1:1376546/1‑100 (MQ=255)
                                                                                       aCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATggg              >  2:3214581/1‑100 (MQ=255)
                                                                                        cGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGc             >  2:710378/1‑100 (MQ=255)
                                                                                        cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGc             >  2:3660401/1‑100 (MQ=255)
                                                                                        cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGc             >  2:1554751/1‑100 (MQ=255)
                                                                                        cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGc             >  1:4373467/1‑100 (MQ=255)
                                                                                         gTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt            >  2:3452616/1‑100 (MQ=255)
                                                                                         gTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt            >  2:233319/1‑100 (MQ=255)
                                                                                         gTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt            >  2:2754717/1‑100 (MQ=255)
                                                                                         gTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt            >  2:875944/1‑100 (MQ=255)
                                                                                         gTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt            >  1:1663223/1‑100 (MQ=255)
                                                                                            gcagatGGTCAGATGTGCCGTATGCCGCAGCCGACACTGCTAGTAATAATTATAAATCCCTGCCGCCATCACCAGGTGGGATGCTACTTCGTGGGCt            <  1:1481850/92‑1 (MQ=255)
                                                                                            gCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCtttt         <  2:2319961/100‑1 (MQ=255)
                                                                                            gCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCtttt         >  2:2638275/1‑100 (MQ=255)
                                                                                             cATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCttttt        <  2:3061370/100‑1 (MQ=255)
                                                                                                 tGGTCAGATGCGACGTTTGCCGCATCGGACACTAGTCGTAATAAGTATAAAGCCCTGCTGCCATCACCAGTTGTGATGCAACTTAATGGGCTTTTTCa      >  1:3783097/1‑98 (MQ=255)
                                                                                                 tGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACg    <  2:2503021/100‑1 (MQ=255)
                                                                                                   gTCAGATGCGACGTTTGCCGCATCGGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATCCCACTTCATGGGCTTGTTCACGcc  >  1:3641921/1‑100 (MQ=255)
                                                                                                   gTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGcc  <  2:4152412/100‑1 (MQ=255)
                                                                                                   gTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGcc  <  1:676788/100‑1 (MQ=255)
                                                                                                   gTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGcc  >  2:1124792/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCC  >  NC_000913/442923‑443121

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: