Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,662,095 C→T A13T (GCA→ACA)  iscR ← DNA‑binding transcriptional dual regulator IscR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,662,0950CT98.5% 986.4 / 1.2 272A13T (GCA→ACA) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (117/151);  total (118/154)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 272 total aligned reads displayed.

AGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGG  >  NC_000913/2661998‑2662188
                                                                                                 |                                                                                             
aGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTgg                                                                                             <  1:1591557/100‑1 (MQ=255)
    gttCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCAc                                                                                         >  1:3796566/1‑100 (MQ=255)
    gttCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCAc                                                                                         >  1:1792253/1‑100 (MQ=255)
       ccAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGc                                                                                      <  1:4237386/100‑1 (MQ=255)
        cAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCa                                                                                     >  1:2021030/1‑100 (MQ=255)
          gATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATa                                                                                   >  1:275430/1‑100 (MQ=255)
          gATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATa                                                                                   >  1:2333692/1‑100 (MQ=255)
           aTAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTgg                                                                                             <  2:2235855/89‑1 (MQ=255)
           aTAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTgg                                                                                             <  2:5872574/89‑1 (MQ=255)
           aTAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAg                                                                                  <  1:446330/100‑1 (MQ=255)
            tAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGTAACGTCAAGCATTGTGGt                                                                                            >  2:4398738/1‑89 (MQ=255)
            tAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAgc                                                                                 <  1:1356386/100‑1 (MQ=255)
              aGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCa                                                                                          >  2:1310020/1‑89 (MQ=255)
               gAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCAc                                                                                         >  2:3137704/1‑89 (MQ=255)
               gAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGcc                                                                              >  1:4128865/1‑100 (MQ=255)
                  aGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGc                                                                                      >  2:4884543/1‑89 (MQ=255)
                   gggAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCa                                                                                     >  2:2419799/1‑89 (MQ=255)
                   gggAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCttt                                                                          >  1:5245702/1‑100 (MQ=255)
                   gggAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCttt                                                                          <  1:5372425/100‑1 (MQ=255)
                   gggAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCttt                                                                          >  1:2398331/1‑100 (MQ=255)
                    ggAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCAt                                                                                    >  2:2796712/1‑89 (MQ=255)
                    ggAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTa                                                                         <  1:5689489/100‑1 (MQ=255)
                    ggAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTa                                                                         >  1:166896/1‑100 (MQ=255)
                      aaaTTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAg                                                                                  <  2:5288615/89‑1 (MQ=255)
                      aaaTTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAgn                                                                       <  1:5627140/100‑2 (MQ=255)
                      aaaTTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGa                                                                       <  1:2932075/100‑1 (MQ=255)
                       aaTTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGAt                                                                      <  1:3024061/100‑1 (MQ=255)
                         ttCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCCGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAgagc                                                                               >  2:6149002/1‑89 (MQ=255)
                          tCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTc                                                                   <  1:4533962/100‑1 (MQ=255)
                            ccTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAg                                                                 >  1:6293111/1‑100 (MQ=255)
                               gACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTa                                                                         <  2:3002293/89‑1 (MQ=255)
                               gACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTa                                                                         >  2:2375041/1‑89 (MQ=255)
                               gACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTa                                                                         <  2:3761112/89‑1 (MQ=255)
                                aCGTTCGGAAATATCATCAAACGGTACCGGGCCAGGTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAg                                                                        <  2:2611423/89‑1 (MQ=255)
                                 cGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGa                                                                       <  2:6059493/89‑1 (MQ=255)
                                 cGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCa                                                            <  1:2518722/100‑1 (MQ=255)
                                 cGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCa                                                            <  1:5538650/100‑1 (MQ=255)
                                  gTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGAt                                                                      >  2:4992878/1‑89 (MQ=255)
                                  gTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGAt                                                                      >  2:2154771/1‑89 (MQ=255)
                                   ttCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATg                                                          <  1:395917/100‑1 (MQ=255)
                                   ttCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATg                                                          <  1:5320176/100‑1 (MQ=255)
                                         aaTATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTAc                                                    <  1:895628/100‑1 (MQ=255)
                                         aaTATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTAc                                                    <  1:2255851/100‑1 (MQ=255)
                                           tatCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACtt                                                  <  1:672402/100‑1 (MQ=255)
                                           tatCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACtt                                                  >  1:4658957/1‑100 (MQ=255)
                                           tatCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGAGCCCTTTAGATGTCAGTCTCATGTCTTACtt                                                  >  1:382227/1‑100 (MQ=255)
                                            atCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCa                                                            >  2:3391621/1‑89 (MQ=255)
                                              cAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATg                                                          >  2:5625257/1‑89 (MQ=255)
                                              cAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATg                                                          <  2:1757622/89‑1 (MQ=255)
                                              cAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATg                                                          <  2:3840829/89‑1 (MQ=255)
                                              cAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCAc                                               >  1:1041151/1‑100 (MQ=255)
                                              cAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCAc                                               >  1:5082835/1‑100 (MQ=255)
                                               aGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGt                                                         >  2:5210646/1‑89 (MQ=255)
                                                 ccAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCt                                                       <  2:2222064/89‑1 (MQ=255)
                                                 ccAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTc                                            <  1:6204390/100‑1 (MQ=255)
                                                  cAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCtt                                                      >  2:2026281/1‑89 (MQ=255)
                                                  cAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCtt                                                      >  2:1488666/1‑89 (MQ=255)
                                                  cAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCtt                                                      <  2:1601604/89‑1 (MQ=255)
                                                   aaCGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAg                                                                 >  2:4851763/1‑77 (MQ=255)
                                                   aaCGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCaa                                          <  1:1950009/100‑1 (MQ=255)
                                                    aCGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCaaa                                         <  1:3693958/100‑1 (MQ=255)
                                                    aCGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTAc                                                    <  2:2850335/89‑1 (MQ=255)
                                                     cGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAAc                                        >  1:2121219/1‑100 (MQ=255)
                                                      ggTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACtt                                                  >  2:4702286/1‑89 (MQ=255)
                                                       gTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTc                                      >  1:3535359/1‑100 (MQ=255)
                                                        tACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCa                                                >  2:2578863/1‑89 (MQ=255)
                                                            gggCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTc                                            >  2:5458694/1‑89 (MQ=255)
                                                            gggCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTc                                            >  2:6369505/1‑89 (MQ=255)
                                                            gggCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGcccc                                 <  1:6316397/100‑1 (MQ=255)
                                                            gggCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGcccc                                 <  1:4013565/100‑1 (MQ=255)
                                                            gggCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCAcccc                                 >  1:1188643/1‑100 (MQ=255)
                                                                   cTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCgg                          <  1:4912088/100‑1 (MQ=255)
                                                                   cTTCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCgg                          <  1:3342586/100‑1 (MQ=255)
                                                                     tCAGAGTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGcc                                   <  2:1353227/89‑1 (MQ=255)
                                                                       agagTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGcccc                                 <  2:3695962/89‑1 (MQ=255)
                                                                       agagTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGcccc                                 >  2:5950128/1‑89 (MQ=255)
                                                                       agagTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGtt                      <  1:3142244/100‑1 (MQ=255)
                                                                         agTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTg                               <  2:4658957/89‑1 (MQ=255)
                                                                          gTTGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGttttt                   >  1:1718966/1‑100 (MQ=255)
                                                                           ttGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGtttttt                  <  1:4815438/100‑1 (MQ=255)
                                                                            tGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAGCTCGCCCCTGccc                            <  2:6253826/89‑1 (MQ=255)
                                                                            tGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTa                 >  1:2201409/1‑100 (MQ=255)
                                                                            tGAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTa                 >  1:6082730/1‑100 (MQ=255)
                                                                             gAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCg                           >  2:5172457/1‑89 (MQ=255)
                                                                             gAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTAt                <  1:823614/100‑1 (MQ=255)
                                                                             gAGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTAt                <  1:1310020/100‑1 (MQ=255)
                                                                              aGCGCAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCgg                          <  2:2286101/89‑1 (MQ=255)
                                                                                 gcAACGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTa            <  1:991127/100‑1 (MQ=255)
                                                                                     cGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGttttt                   <  2:4193771/89‑1 (MQ=255)
                                                                                     cGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGttttt                   >  2:378220/1‑89 (MQ=255)
                                                                                     cGTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAg        >  1:1156416/1‑100 (MQ=255)
                                                                                      gTCAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGtttttt                  >  2:235980/1‑89 (MQ=255)
                                                                                      gTCAAGCATTGTGGTCACGCCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGtttttt                  <  2:444765/89‑1 (MQ=255)
                                                                                        cAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTAt                <  2:2540059/89‑1 (MQ=255)
                                                                                        cAAGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTgg     >  1:1277301/1‑100 (MQ=255)
                                                                                          aGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTg              <  2:3938786/89‑1 (MQ=255)
                                                                                          aGCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTgggg   >  1:5461367/1‑100 (MQ=255)
                                                                                           gCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGt             >  2:5752517/1‑89 (MQ=255)
                                                                                           gCATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCACGGGGTTTTTTATTGTAAAAGTggggg  <  1:2575430/100‑1 (MQ=255)
                                                                                            cATTGTGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTa            <  2:166896/89‑1 (MQ=255)
                                                                                                 |                                                                                             
AGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGG  >  NC_000913/2661998‑2662188

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