Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,323 A→T 15.2% intergenic (‑150/+17) ygfI ← / ← yggE putative DNA‑binding transcriptional regulator/oxidative stress defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,067,3230AT15.2% 467.9 / 68.6 185intergenic (‑150/+17)ygfI/yggEputative DNA‑binding transcriptional regulator/oxidative stress defense protein
Reads supporting (aligned to +/- strand):  ref base A (70/86);  new base T (12/16);  total (82/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

Only 100 of 190 total aligned reads displayed.

TGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAA  >  NC_000913/3067224‑3067422
                                                                                                   |                                                                                                   
tgatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATa                                                                                                     >  1:2090164/1‑100 (MQ=255)
tgatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATa                                                                                                     <  2:2472817/100‑1 (MQ=255)
 gatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATAc                                                                                                    >  1:1016317/1‑100 (MQ=255)
 gatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAAAAAAAAAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATAc                                                                                                    <  2:1458028/100‑1 (MQ=255)
   tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGt                                                                                                  >  1:3749329/1‑100 (MQ=255)
   tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGt                                                                                                  >  1:272970/1‑100 (MQ=255)
      aaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                                               <  1:970120/100‑1 (MQ=255)
      aaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                                               <  1:1273683/100‑1 (MQ=255)
       aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCggg                                                                                              >  2:4006022/1‑100 (MQ=255)
       aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                             <  1:2295004/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCTGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                              <  1:600657/99‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                              <  1:2634680/99‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                             <  1:3374452/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGc                                                                                             <  2:1153452/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGACGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggcc                                                                                             <  1:2794696/100‑3 (MQ=255)
         cgggAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATAGCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                             <  2:1050555/96‑1 (MQ=255)
         cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                             <  2:3311174/99‑1 (MQ=255)
         cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCt                                                                                            <  1:3605425/100‑1 (MQ=255)
         cgcgAGAAAATATAGGTGATTACGTCATGGTTTTTCAGAAAATTACAGAAACAGGAGGCAATATCGGTTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                              <  2:4126756/98‑1 (MQ=255)
           cgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTAc                                                                                          >  2:2402084/1‑100 (MQ=255)
             aGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCTaatg                                                                                        <  1:3230627/100‑4 (MQ=39)
              tcAACGATAGGTGATCACGTCATCGTCTTACAGAAATTTACAGAAAAAGGAGGAAACAACGGGTAAAGGCATTACCCCGACGAATAAGTCGCGCTACaa                                                                                        >  2:3666469/2‑99 (MQ=255)
              gCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACaaa                                                                                       <  1:2571254/100‑1 (MQ=255)
               cAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCAGTAGCCCGACGTATTCGTCGGGCTaatgcc                                                                                      <  1:3418036/100‑6 (MQ=38)
                naaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAata                                                                                     >  2:882432/2‑100 (MQ=255)
                aaaGTTATGTGATTACGTCATGGTTTTTCAGAAAATTACAGAAAAAGGAGGCAATATCGGCTAAAGGCATTAGCCCGACGAATACGTCGGGCTAc                                                                                          <  2:3413273/95‑1 (MQ=255)
                   gATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAAtatta                                                                                  >  1:570754/1‑100 (MQ=255)
                    aTAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAAtattat                                                                                 <  1:4009655/100‑1 (MQ=255)
                     tAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATattatt                                                                                <  1:1032121/100‑1 (MQ=255)
                      aGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtg                                                                               <  2:1497023/100‑1 (MQ=255)
                      aGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtg                                                                               <  1:3709789/100‑1 (MQ=255)
                      aGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtg                                                                               <  1:3506120/100‑1 (MQ=255)
                       ggTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATattatt                                                                                <  2:3160210/98‑1 (MQ=255)
                         tGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                            <  2:961074/100‑1 (MQ=255)
                           aTTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgctg                                                                          <  2:3985449/100‑1 (MQ=255)
                            ttACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgctgc                                                                         >  1:2766386/1‑100 (MQ=255)
                             tACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCa                                                                        >  2:1269590/1‑100 (MQ=255)
                              aCGTGATGGTTTTACAGAAAATTAAAGAAAAAGGAGGCAATATCGGGTAAATGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAg                                                                       <  2:3771016/100‑1 (MQ=255)
                               cGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAgg                                                                      <  2:4106707/100‑1 (MQ=255)
                               cGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATACCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAgg                                                                      >  2:2094266/1‑100 (MQ=255)
                                gTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgt                                                                     >  2:3071054/1‑100 (MQ=255)
                                   aTGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtt                                                                  <  1:3038/100‑1 (MQ=255)
                                    tGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGttt                                                                 >  1:1282072/1‑100 (MQ=255)
                                    tGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGttt                                                                 >  2:2843801/1‑100 (MQ=255)
                                       ttttACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAg                                                              >  1:3971766/1‑100 (MQ=255)
                                          tacagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCgg                                                           <  2:273304/100‑1 (MQ=255)
                                            cagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgt                                                         >  1:1421492/1‑100 (MQ=255)
                                              gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttg                                                       >  2:639824/1‑100 (MQ=255)
                                               aaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgt                                                      <  2:362477/100‑1 (MQ=255)
                                                aaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                                     <  1:2889705/100‑1 (MQ=255)
                                                 aattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                    >  1:1503249/1‑100 (MQ=255)
                                                  attacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGa                                                   <  2:466110/100‑1 (MQ=255)
                                                   ttacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGAt                                                  <  1:2457777/100‑1 (MQ=255)
                                                     ncagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcc                                                >  2:1812381/2‑100 (MQ=255)
                                                      cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                               >  1:3287918/1‑100 (MQ=255)
                                                      cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                               >  1:3202784/1‑100 (MQ=255)
                                                          aaaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAgg                                           >  1:846454/1‑100 (MQ=255)
                                                          aaaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAgg                                           >  1:3253658/1‑100 (MQ=255)
                                                            aaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGtt                                         <  1:882438/100‑1 (MQ=255)
                                                            aaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGtt                                         <  1:286905/100‑1 (MQ=255)
                                                              ggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGATGTTGATCCACAGGTTCt                                       >  2:3840215/1‑100 (MQ=255)
                                                                 ggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAAc                                    >  1:2416042/1‑100 (MQ=255)
                                                                   cAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTg                                  >  2:1678449/1‑100 (MQ=255)
                                                                     atatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGa                                <  2:2991125/100‑1 (MQ=255)
                                                                     atatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGa                                <  2:1330244/100‑1 (MQ=255)
                                                                     atatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGa                                <  2:2680398/100‑1 (MQ=255)
                                                                       atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                              >  1:1724438/1‑100 (MQ=255)
                                                                       atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                              <  1:1218701/100‑1 (MQ=255)
                                                                         cGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAc                            <  1:1559366/100‑1 (MQ=255)
                                                                          gggTAAAGGCATTAGGCCGACGAATACGTCGGGCTACGAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAcc                           <  1:4074694/100‑1 (MQ=255)
                                                                          gggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTGACTGGAAGAcc                           <  1:196550/100‑1 (MQ=255)
                                                                          gggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAcc                           <  1:379032/100‑1 (MQ=255)
                                                                          gggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGGGTTTTAGCGGGGTATTGATCCACAGTTTCTAACTGGAAGTcc                           >  2:2045350/1‑100 (MQ=255)
                                                                           ggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACca                          <  2:12871/100‑1 (MQ=255)
                                                                           ggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACca                          <  2:3077536/100‑1 (MQ=255)
                                                                           ggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGAGGATTGTTGATCCACAGGTTCTAACTGGAAGAcaa                          >  1:3117630/1‑98 (MQ=255)
                                                                              aaaGGCATTAGCCCTACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACAt                       <  2:1094981/100‑1 (MQ=255)
                                                                                    aTTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCt                 >  1:1241226/1‑100 (MQ=255)
                                                                                    aTTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCt                 <  2:1971121/100‑1 (MQ=255)
                                                                                       aGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGat              >  1:833111/1‑100 (MQ=255)
                                                                                       aGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACTACATCGACCTGat              <  1:2807796/100‑1 (MQ=255)
                                                                                       aGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGat              >  2:100552/1‑100 (MQ=255)
                                                                                       aGCCCGACGAATACGTCGGGCTACAAATATTATTCTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGat              <  1:1322434/100‑1 (MQ=255)
                                                                                        gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc             >  1:493128/1‑100 (MQ=255)
                                                                                        gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc             >  1:3173148/1‑100 (MQ=255)
                                                                                        gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc             >  1:1688459/1‑100 (MQ=255)
                                                                                        gCCCGACGTATTAGTCGGGCCACAAATATTATTGTGCTGCAGGTGTTTTAGGGAGTTGTTGTTCGCCAGGGTTTAACTGGAAGACGACATCGACTTg                >  2:2532145/1‑97 (MQ=255)
                                                                                        gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc             <  1:1984089/100‑1 (MQ=255)
                                                                                         cccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatca            <  1:2245545/100‑1 (MQ=255)
                                                                                         cccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatca            <  1:2410606/100‑1 (MQ=255)
                                                                                         cccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatca            <  2:123524/100‑1 (MQ=255)
                                                                                         cccGACGAATACGTCGGGCTACAAATATTATTGGGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatca            <  2:1598070/100‑1 (MQ=255)
                                                                                             acgaATACGCCGGGCTACAAAAATTATAGTGCTGCCGATGTTTCAGCGGCAAGTTGAACTACAGCCTCTAACTGGAAGACCGAATCGCCCCGttcatcca        >  2:1349931/1‑98 (MQ=255)
                                                                                              cgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCAACAGATTCTAACTGGAAGACCACATCGACCTGATCATCaaa       >  2:3771874/1‑100 (MQ=255)
                                                                                               gaATACGTCGGGCGACAAATATTATTGTGGTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAAc      <  2:2873192/100‑1 (MQ=255)
                                                                                                aaTACGTCGTGCTACAAACATTATTGTGCTGCAGGTATTTTAGCGCGTTGTTGATCCACAGGTTCTAACTGGAATACCACATCGACCTGATCATCAAACt     <  2:1877639/100‑1 (MQ=255)
                                                                                                aaTACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACt     <  2:2063261/100‑1 (MQ=255)
                                                                                                 aTACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTg    <  2:4043892/100‑1 (MQ=255)
                                                                                                  tACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGa   >  1:1841822/1‑100 (MQ=255)
                                                                                                   aCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACAACATTGACCTGATCATCAAACTGaa  >  2:797836/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAA  >  NC_000913/3067224‑3067422

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: