Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 843,749 C→A 100% D578Y (GAT→TAT)  ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,7490CA97.6% 551.3 / 3.6 170D578Y (GAT→TAT) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (85/81);  total (87/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01

Only 100 of 170 total aligned reads displayed.

GAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAAT  >  NC_000913/843653‑843843
                                                                                                |                                                                                              
gAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATcc                                                                                             >  2:3454419/1‑100 (MQ=255)
  atgatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAg                                                                                           <  2:3957951/100‑1 (MQ=255)
    gatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtg                                                                                         >  1:2577990/1‑100 (MQ=255)
     atgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtgt                                                                                        >  2:3745695/1‑100 (MQ=255)
       gtgATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTc                                                                                      >  2:3415173/1‑100 (MQ=255)
       gtgATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTc                                                                                      <  1:3583183/100‑1 (MQ=255)
       gtgATACGCCCCGGTATCCTGACGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGAGtt                                                                                       >  1:4208941/1‑99 (MQ=255)
         gATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCAt                                                                                    >  2:2405733/1‑100 (MQ=255)
           tACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATCGTCACCAAATATCCAGTGTTCATGc                                                                                  >  1:1360751/1‑100 (MQ=255)
            aCGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGcc                                                                                 >  1:601868/1‑100 (MQ=255)
              gCCCCGGCATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGCTAGTCCCCACATATCCTGTGTTCATGcc                                                                                 >  2:2304539/1‑98 (MQ=255)
                 ccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGtttt                                                                            <  2:507631/100‑1 (MQ=255)
                 ccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGtttt                                                                            >  2:4335228/1‑100 (MQ=255)
                 ccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGtttt                                                                            >  1:3522540/1‑100 (MQ=255)
                 ccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGttt                                                                             <  1:2881180/99‑1 (MQ=255)
                 ccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGCGCCGGTCAACGGCCCGATAGTCACCAAGTATCCAGTGTGCATGCCGt                                                                               <  1:3758706/97‑1 (MQ=255)
                  cGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTc                                                                           <  1:4335197/100‑1 (MQ=255)
                   ggTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCa                                                                          <  1:2136928/100‑1 (MQ=255)
                    gTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCGAGTGTTCATGCCGTTTTCaa                                                                         <  1:2685822/100‑1 (MQ=255)
                    gTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCaa                                                                         <  2:930322/100‑1 (MQ=255)
                      aTCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAAc                                                                       <  2:700066/100‑1 (MQ=255)
                       tCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCCCCAAATATCCATTGATCATGCCGATTTCAAACt                                                                      >  1:423073/1‑100 (MQ=255)
                        ccTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGCCACCAAATATCCAGTGTTCATGCCGCTCTCa                                                                          >  1:3866650/1‑95 (MQ=255)
                          tGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaa                                                                   <  1:3649682/100‑1 (MQ=255)
                          tGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaa                                                                   <  2:1752844/100‑1 (MQ=255)
                          tGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaa                                                                   <  2:2012995/100‑1 (MQ=255)
                           gTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaat                                                                  >  1:87906/1‑100 (MQ=255)
                             cGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCACCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataa                                                                >  1:2435252/1‑100 (MQ=255)
                                 cGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaat                                                              >  1:1785584/1‑98 (MQ=255)
                                 cGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAc                                                            <  1:1986403/100‑1 (MQ=255)
                                 cGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAc                                                            >  2:74743/1‑100 (MQ=255)
                                  gCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAcc                                                           <  2:1407328/100‑1 (MQ=255)
                                   cccACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCCGTGTTCATGCCGTTTTCAAACTGAATAAATAccc                                                          >  1:2181314/1‑100 (MQ=255)
                                     cACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCg                                                        >  1:987292/1‑100 (MQ=255)
                                      aCGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCgg                                                       <  1:4161765/100‑1 (MQ=255)
                                      aCGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCgg                                                       <  1:70026/100‑1 (MQ=255)
                                       cGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGt                                                      <  2:2454316/100‑1 (MQ=255)
                                       cGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGt                                                      <  2:2559287/100‑1 (MQ=255)
                                       cGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGt                                                      >  2:3779846/1‑100 (MQ=255)
                                        ggaactaaTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGTTg                                                     <  1:1157174/94‑1 (MQ=255)
                                        ggAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTg                                                     <  2:3588032/100‑1 (MQ=255)
                                        ggAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTg                                                     <  2:3480435/100‑1 (MQ=255)
                                        ggAGCGAATCGACATCCGTTCCACTGTGCAGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTg                                                     <  1:2501969/100‑1 (MQ=255)
                                             gAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataat                                                >  2:31902/1‑100 (MQ=255)
                                             gAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataat                                                >  1:1870056/1‑100 (MQ=255)
                                             gAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATACTCACCAAATATCCAGCGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGata                                                  >  1:689543/1‑98 (MQ=255)
                                              aaTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataatn                                               <  2:4202190/100‑2 (MQ=255)
                                              aaTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataata                                               <  2:300330/100‑1 (MQ=255)
                                              aaTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataata                                               <  2:4376841/100‑1 (MQ=255)
                                                 cGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCACTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCt                                            <  1:2597066/100‑1 (MQ=255)
                                                  gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCtt                                           <  2:4150107/100‑1 (MQ=255)
                                                  gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCtt                                           >  2:1385078/1‑100 (MQ=255)
                                                  gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCtt                                           >  2:1597803/1‑100 (MQ=255)
                                                  gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCtt                                           >  1:784820/1‑100 (MQ=255)
                                                    cATCCGTTCCACTGTGTCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTc                                         <  2:4292420/100‑1 (MQ=255)
                                                    cATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTc                                         >  1:1665812/1‑100 (MQ=255)
                                                    cATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTc                                         <  1:3024647/100‑1 (MQ=255)
                                                       ccGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAcc                                      >  2:3914186/1‑100 (MQ=255)
                                                          ttCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGc                                   <  1:2242926/100‑1 (MQ=255)
                                                          ttCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGc                                   <  2:2245326/100‑1 (MQ=255)
                                                          tgccACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGCGTTCATGCTGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAc                                       <  1:1947971/94‑1 (MQ=255)
                                                           tCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCg                                  >  1:465526/1‑100 (MQ=255)
                                                           tCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCg                                  >  2:3601796/1‑100 (MQ=255)
                                                            ccACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGt                                 >  1:180990/1‑100 (MQ=255)
                                                             cACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGtt                                >  2:1315228/1‑100 (MQ=255)
                                                              aCTGTGCCGGTAAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCACGGTGATAATATCTTTCACCAGCGttt                               >  1:3202884/1‑100 (MQ=255)
                                                                tgtgCCGGTCAACGGCGCGTTAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTATGc                             <  1:3271998/100‑1 (MQ=255)
                                                                tgtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGc                             <  2:3157980/100‑1 (MQ=255)
                                                                 gtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATATCCCGGTGATAATATCTTTCACCAGCGTTTGCg                            <  1:1853159/100‑1 (MQ=255)
                                                                 gtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                            >  2:21827/1‑100 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCCAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                            >  2:875722/1‑98 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                          >  2:3739334/1‑100 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGa                           >  1:3631884/1‑99 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTGTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                            <  1:3360334/98‑1 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGCTTTCAAACCGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCgca                          >  1:1931434/1‑98 (MQ=255)
                                                                    ccGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAAc                         >  2:2733366/1‑100 (MQ=255)
                                                                    ccGGTCAAAGGCCAGATAGTCACCCAATATCCAGTGTTCATGCCGTTTTCAAACTGTATTAAAACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                          <  2:644458/99‑1 (MQ=255)
                                                                     cGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAAcc                        <  1:2945673/100‑1 (MQ=255)
                                                                     cGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACACCGGTGATAATATCTTTCACCAGCGttt                               >  1:2355633/1‑93 (MQ=255)
                                                                      ggTCACCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCCAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCa                       >  2:3755239/1‑100 (MQ=255)
                                                                      ggTCAACGGCCCGATAGTCACCAAATATCCGGTGTTCATGCCGCCTTGAATCTGAATAAATACCCCGGTGATGATAGCTTTCACCAGCGCTTGCGAACCa                       >  1:4005409/1‑100 (MQ=255)
                                                                         cAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAAc                    >  2:1476537/1‑100 (MQ=255)
                                                                         cAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAAc                    <  1:1187932/100‑1 (MQ=255)
                                                                          aaCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACg                   >  2:2132024/1‑100 (MQ=255)
                                                                           aCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACga                  >  2:2395486/1‑100 (MQ=255)
                                                                            cGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACgag                 >  2:1561443/1‑100 (MQ=255)
                                                                            cGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACgag                 <  1:1558598/100‑1 (MQ=255)
                                                                               cccGTTAGGCACCAAATATCCAGTGTTCATGCCGGTTTCAAACTAAGTAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACgag                 <  1:3898412/97‑1 (MQ=255)
                                                                                ccGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAg             >  2:1601739/1‑100 (MQ=255)
                                                                                  gATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcc           >  2:1149729/1‑100 (MQ=255)
                                                                                   aTAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCa          <  2:3493297/100‑1 (MQ=255)
                                                                                       tCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAcc      <  1:2351705/100‑1 (MQ=255)
                                                                                       tCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAcc      <  1:2634550/100‑1 (MQ=255)
                                                                                        cACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcagatcgg     >  1:3890213/1‑93 (MQ=255)
                                                                                        cACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCt     >  1:1592137/1‑100 (MQ=255)
                                                                                         aCCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCt     >  2:2721025/1‑99 (MQ=255)
                                                                                         aCCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTa    >  2:3318582/1‑100 (MQ=255)
                                                                                         aCCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTa    >  2:2865518/1‑100 (MQ=255)
                                                                                           cAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAAt  >  2:562434/1‑100 (MQ=255)
                                                                                            aaaTATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAACGTTTACAAACCAAACAAGATAGCCAGACCTaa   >  1:481383/1‑98 (MQ=255)
                                                                                                |                                                                                              
GAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAAT  >  NC_000913/843653‑843843

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