Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,224,239 A→T 21.1% intergenic (+332/+40) ycgI → / ← minE pseudogene/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,224,2390AT21.1% 283.0 / 71.1 133intergenic (+332/+40)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (47/57);  new base T (20/8);  total (68/65)
Fisher's exact test for biased strand distribution p-value = 1.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 147 total aligned reads displayed.

TGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCC  >  NC_000913/1224140‑1224338
                                                                                                   |                                                                                                   
  tAAAGGGGGAGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAgaga                                                                                                     <  2:2892378/98‑1 (MQ=255)
  tAAAGGGGGAGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATa                                                                                                   <  2:1507827/100‑1 (MQ=255)
  tAAAGGGGGAGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATa                                                                                                   <  2:2090839/100‑1 (MQ=255)
  tAAAGGGGGAGGAGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTCCCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATa                                                                                                   >  1:1295486/1‑100 (MQ=255)
   aaaGGGGGAGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATaa                                                                                                  >  1:672566/1‑100 (MQ=255)
    aaGGGGGAGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAc                                                                                                 <  1:1531423/100‑1 (MQ=255)
    aaGGGGGAGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAc                                                                                                 >  1:1549604/1‑100 (MQ=255)
    aaGGGGGAGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAc                                                                                                 >  1:2399887/1‑100 (MQ=255)
    aaGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAc                                                                                                 >  2:140583/1‑100 (MQ=255)
        gggAGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAgagagatcggaa                                                                                             >  1:2248664/1‑92 (MQ=255)
          taGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                                                           <  1:846754/99‑1 (MQ=255)
          gAGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATATTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                                                           <  1:4125757/100‑1 (MQ=255)
          gAGGTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                                                           <  1:686628/100‑1 (MQ=255)
            ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                         >  2:1666204/1‑100 (MQ=255)
            ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                         <  2:4161035/100‑1 (MQ=255)
             gtggagagtgCAAAGAGAGGTATCATGATGGATACTCGCCGCTTTACCACAACGTGGATGATAGAGATGACAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                        <  2:4021217/91‑1 (MQ=255)
             gTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                        <  1:301626/100‑1 (MQ=255)
             gTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                        <  1:741876/100‑1 (MQ=255)
             gTGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                                                        <  1:1948747/100‑1 (MQ=255)
             gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                        <  1:1213028/100‑1 (MQ=255)
              tGGGAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAACCTCAAGGCAGAGTTATCTCTGCCTTGa                                                                                       <  1:3144858/100‑1 (MQ=39)
                ggaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagt                                                                                     <  1:111659/98‑3 (MQ=39)
                ggaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                                                     <  1:1777598/98‑1 (MQ=255)
                ggaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                                                     <  1:1340772/98‑1 (MQ=255)
                 gaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGa                                                                                    <  2:705677/99‑1 (MQ=255)
                    tagaagtGAGGAGTATCAGTAAGAAGAATAGGCGCTTGACGACAACGTGGATAAGAAAAAAGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                      <  2:1306356/89‑1 (MQ=255)
                    tAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaa                                                                                 >  2:2558445/1‑100 (MQ=255)
                     aGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaaa                                                                                >  2:310484/1‑100 (MQ=255)
                      gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                               >  1:3086650/1‑100 (MQ=255)
                       cAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATg                                                                              <  1:3205933/100‑1 (MQ=255)
                        aaTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                                                             >  1:2440358/1‑100 (MQ=255)
                        aaTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                                                             >  2:462816/1‑100 (MQ=255)
                             ggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCttt                                                                        >  1:1972509/1‑100 (MQ=255)
                                   tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGc                                                                  <  1:4132679/100‑1 (MQ=255)
                                   tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGc                                                                  <  2:1658158/100‑1 (MQ=255)
                                     aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                                >  1:2087622/1‑100 (MQ=255)
                                     aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                                <  1:3980508/100‑1 (MQ=255)
                                      gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                               >  1:3297134/1‑100 (MQ=255)
                                      gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                               >  1:3136025/1‑100 (MQ=255)
                                       cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                                                              >  2:1939910/1‑100 (MQ=255)
                                         aGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCtt                                                            >  2:1498128/1‑100 (MQ=255)
                                         aGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCtt                                                            >  1:3969024/1‑100 (MQ=255)
                                          gAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTa                                                           <  1:3222420/100‑1 (MQ=255)
                                           aaTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                                                          <  1:4116317/100‑1 (MQ=255)
                                            aTACTCACCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAtt                                                         >  1:3740699/1‑100 (MQ=255)
                                              aCTCGACGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                                                       <  2:3980460/100‑1 (MQ=255)
                                               cTCGCCGCTTTACCACACCGTGCATGCGAGGGCTGCAACACTCAAGCCCCAGATAACTCTGCCTTGAAGATAAATCCGCTTTTACAGCGGGCTTCTTTc                                                       >  2:3454504/1‑99 (MQ=255)
                                               cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                                                       >  2:785753/1‑99 (MQ=255)
                                               cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                      >  1:2325624/1‑100 (MQ=255)
                                               cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                      >  1:3888486/1‑100 (MQ=255)
                                               cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGATTTTACAGCGGGCTTATTTCa                                                      >  2:182685/1‑100 (MQ=255)
                                                tCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                     <  1:632728/100‑1 (MQ=255)
                                                 cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                     <  2:3334903/99‑1 (MQ=255)
                                                 cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                    <  1:2844675/100‑1 (MQ=255)
                                                 cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                    <  1:3630052/100‑1 (MQ=255)
                                                    cgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGctct                                                 <  2:1737860/100‑1 (MQ=255)
                                                      cttcatcacaACGTCGATGAGAGGGATGAAATACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                     <  2:2482491/88‑1 (MQ=255)
                                                      cTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                               <  2:3132947/100‑1 (MQ=255)
                                                         tACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGc                                            <  2:2261435/100‑1 (MQ=255)
                                                         tACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGc                                            <  1:1693690/100‑1 (MQ=255)
                                                          aCCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCt                                           >  2:3779393/1‑100 (MQ=255)
                                                           ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                          >  2:3416400/1‑100 (MQ=255)
                                                                aCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt                                     >  1:632033/1‑100 (MQ=255)
                                                                    ggATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAgg                                 <  2:1253664/100‑1 (MQ=255)
                                                                      aTGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTc                               <  1:2555464/100‑1 (MQ=255)
                                                                       tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCa                              >  2:986917/1‑100 (MQ=255)
                                                                            gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                         <  1:960344/100‑1 (MQ=255)
                                                                            gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGCCTTATTTCAGCTCTTCTTCTTCCGGTAAGGTCACGtt                          >  2:2259439/1‑99 (MQ=255)
                                                                             ggATGAAAAACTCAAGGCAGAGATAACTCTGCCTCGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCa                        <  1:3525283/100‑1 (MQ=255)
                                                                              gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAg                       >  2:3830392/1‑100 (MQ=255)
                                                                               aTGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGc                      >  1:426564/1‑100 (MQ=255)
                                                                                tGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCt                     >  1:1811947/1‑100 (MQ=255)
                                                                                 gAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTc                    <  1:915626/100‑1 (MQ=255)
                                                                                  aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCa                   >  2:1258702/1‑100 (MQ=255)
                                                                                    aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAg                 <  1:2445141/100‑1 (MQ=255)
                                                                                    aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAg                 >  2:2753991/1‑100 (MQ=255)
                                                                                     tcttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                >  1:3638749/4‑100 (MQ=39)
                                                                                      cttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                >  2:2146444/3‑99 (MQ=255)
                                                                                       ttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAAt              >  2:3767931/2‑100 (MQ=255)
                                                                                       ttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAAt              >  1:2117580/2‑100 (MQ=255)
                                                                                       cTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAAt              >  1:3920034/1‑100 (MQ=255)
                                                                                        tCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa             >  1:2224584/1‑100 (MQ=255)
                                                                                         ccaGGCGGAGCTAACTCTGCCTTGAAGAGAAATGCGCTTTTACAGAGTGCTTATGTAGGCTCTTCTGCTGCTGCTAAGGGCACGTCCAGCTCACGa                >  2:3860200/3‑96 (MQ=255)
                                                                                         cAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAg            >  2:3950453/1‑100 (MQ=255)
                                                                                         cAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAg            >  1:857429/1‑100 (MQ=255)
                                                                                         cAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAg            <  1:2245042/100‑1 (MQ=255)
                                                                                          aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  2:3021464/1‑100 (MQ=255)
                                                                                          aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  1:1469557/1‑100 (MQ=255)
                                                                                          aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  2:394998/1‑100 (MQ=255)
                                                                                          aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  1:3294747/1‑100 (MQ=255)
                                                                                          aaGGCAGAGATAATTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  2:2490934/1‑100 (MQ=255)
                                                                                          aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           <  1:3698512/100‑1 (MQ=255)
                                                                                          aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           <  1:720648/100‑1 (MQ=255)
                                                                                            ggCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaaa         <  1:3586335/100‑1 (MQ=255)
                                                                                                 agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCg    <  1:3474897/100‑1 (MQ=255)
                                                                                                 agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCg    >  2:721792/1‑100 (MQ=255)
                                                                                                  gaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGc   >  2:3591291/1‑100 (MQ=255)
                                                                                                  gaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGc   >  2:3264322/1‑100 (MQ=255)
                                                                                                  gaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGc   >  1:751497/1‑100 (MQ=255)
                                                                                                   aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGcc  <  1:1084154/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCC  >  NC_000913/1224140‑1224338

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: