Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,825 G→C 15.1% P88A (CCC→GCC)  ftsX ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,8250GC15.1% 329.9 / 41.3 117P88A (CCC→GCC) ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (51/48);  new base C (11/7);  total (62/55)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 122 total aligned reads displayed.

GTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTG  >  NC_000913/3601727‑3601924
                                                                                                  |                                                                                                   
gTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGc                                                                                                    <  2:372064/100‑1 (MQ=255)
    aGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCg                                                                                                <  1:1510817/100‑1 (MQ=255)
     gCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATCCCGCGTCGCCGCCTGGTTAACGTTTTCGTACACCATATAACAGACGCTGGGCCGCGt                                                                                               >  2:2147524/1‑100 (MQ=255)
         tttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCaga                                                                                           >  1:146359/1‑100 (MQ=255)
         tttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgt                                                                                           <  1:328536/100‑4 (MQ=255)
         tttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgt                                                                                           <  1:2409205/100‑4 (MQ=255)
          ttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCagag                                                                                          <  1:1555577/100‑1 (MQ=255)
             tcgaggcaaaCAGTGATTTGCGGTGACCTATAATACTGCCTCCCGCCCTGGTCAACGTTGTTGTACACCATACAACAGTTTCTAGGCAGCTTCag                                                                                            <  2:1689586/88‑1 (MQ=255)
             tGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGTTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgttata                                                                                       <  2:2167582/100‑8 (MQ=255)
              gCAGATAAACAGTGATTTGCGGTGACTGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgttatat                                                                                      <  2:2473356/100‑9 (MQ=255)
              gCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAAt                                                                                      >  2:570734/1‑100 (MQ=255)
              gCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAAt                                                                                      >  1:2131328/1‑100 (MQ=255)
              gCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAAt                                                                                      >  1:848562/1‑100 (MQ=255)
               cAGATAAACAGTGATTTGCGGTGACGGATAATACTGCCTCGCCGCCTGGGTAACGTTTTTGTACACCATATAACAGACGCTGGCCAGGGTCAGAGAAATg                                                                                     >  2:2407025/1‑100 (MQ=255)
                aGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgttatatgg                                                                                    <  1:1282879/100‑11 (MQ=38)
                aGATAAACAGTGAGTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgttatatg                                                                                     <  2:2488453/99‑10 (MQ=38)
                 gATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGc                                                                                   <  1:2099074/100‑1 (MQ=255)
                 gATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGc                                                                                   >  2:1949068/1‑100 (MQ=255)
                  aTAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCg                                                                                  <  1:383927/100‑1 (MQ=255)
                     aaCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATa                                                                               >  1:1954615/1‑100 (MQ=255)
                      aCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATaa                                                                              >  2:896938/1‑100 (MQ=255)
                       cAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAc                                                                             >  1:918106/1‑100 (MQ=255)
                       cAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAc                                                                             >  2:859870/1‑100 (MQ=255)
                        aGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAcc                                                                            <  1:431450/100‑1 (MQ=255)
                         gTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGATAAATGGCGATAACCa                                                                           <  2:76304/100‑1 (MQ=255)
                           gATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATc                                                                         <  2:2186982/100‑1 (MQ=255)
                           gATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATc                                                                         >  1:1334616/1‑100 (MQ=255)
                              ttGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCAc                                                                       <  1:2241692/99‑1 (MQ=255)
                               tGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCg                                                                     <  1:1671831/100‑1 (MQ=255)
                                gCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCGGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGt                                                                    <  2:1619893/100‑1 (MQ=255)
                                 cGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGtt                                                                   <  2:802786/100‑1 (MQ=255)
                                 cGGTGACGGATAATACTGCGTCACCGACTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGtt                                                                   >  2:706792/1‑100 (MQ=255)
                                     gACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTaaaa                                                               >  2:2252726/1‑100 (MQ=255)
                                      aCGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTaaaaa                                                              >  1:2213948/1‑100 (MQ=255)
                                        ggATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACg                                                            <  1:946425/100‑1 (MQ=255)
                                        ggATAATACTGCGTCACCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACg                                                            >  1:2144812/1‑100 (MQ=255)
                                          ataataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTg                                                          >  1:2322371/1‑100 (MQ=255)
                                            aataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGc                                                        <  2:1830/100‑1 (MQ=255)
                                             ataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCg                                                       <  1:1913749/100‑1 (MQ=255)
                                             ataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCCACAGAGAAATGGCGATGACCATCACGGTTAAACACCTgga                                                        >  1:2394112/1‑98 (MQ=255)
                                              taCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGa                                                      >  1:803147/1‑100 (MQ=255)
                                              taCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGa                                                      >  2:2272744/1‑100 (MQ=255)
                                               aCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTACAAACGTGGCGac                                                     >  1:1491614/1‑99 (MQ=255)
                                               aCTGCGTCGCCGCCTGGTTAACGTTTCTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAGTGGCGATAACCATCACCGTTAAAAACGTGGCGaa                                                     >  2:904386/1‑100 (MQ=255)
                                                cTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAAc                                                    <  1:1927248/100‑1 (MQ=255)
                                                cTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAAc                                                    >  2:2052335/1‑100 (MQ=255)
                                                  gCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACgg                                                  >  2:1092488/1‑100 (MQ=255)
                                                    gTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGtt                                                >  2:1401717/1‑100 (MQ=255)
                                                     tCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGttt                                               <  1:1911917/100‑1 (MQ=255)
                                                     tCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGttt                                               <  1:1912430/100‑1 (MQ=255)
                                                      cgccgcCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTg                                              <  1:859878/100‑1 (MQ=255)
                                                      cgccgcCTGGTTAACGTTTTTGTACACCATATAACAGACGCCGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGttt                                               <  1:2271038/99‑1 (MQ=255)
                                                       gccgccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGc                                             <  1:1355516/100‑1 (MQ=255)
                                                       gccgccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGc                                             <  2:2326710/100‑1 (MQ=255)
                                                         cgcgTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGCGAACTTGCGATAAACATCACAGTTAAAAACGTGGCGAACCGTTTGCtt                                           >  2:2176817/1‑100 (MQ=255)
                                                         cgccTGGTTAACGTTTTTGTACACCATATAACAGACGCGGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGc                                             <  2:669914/98‑1 (MQ=255)
                                                          gccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCttt                                          <  1:848847/100‑1 (MQ=255)
                                                          gccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCttt                                          <  2:1951993/100‑1 (MQ=255)
                                                           ccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCtttt                                         <  1:607109/100‑1 (MQ=255)
                                                            cTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTc                                        <  1:1180121/100‑1 (MQ=255)
                                                            cTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTc                                        <  2:379080/100‑1 (MQ=255)
                                                              ggTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAg                                      <  2:1214871/100‑1 (MQ=255)
                                                               gTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGa                                     >  2:2328306/1‑100 (MQ=255)
                                                                ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGAt                                    <  1:867553/100‑1 (MQ=255)
                                                                 tAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATc                                   >  1:1757895/1‑100 (MQ=255)
                                                                 tAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATc                                   <  1:2257482/100‑1 (MQ=255)
                                                                  aaCGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATcc                                  <  1:2154180/100‑1 (MQ=255)
                                                                    catttttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTg                                >  1:1730476/3‑100 (MQ=255)
                                                                    cGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTg                                >  1:1904354/1‑100 (MQ=255)
                                                                     gTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCTATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGc                               <  1:721281/100‑1 (MQ=255)
                                                                     gTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGc                               <  1:1337454/100‑1 (MQ=255)
                                                                      tttttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCa                              >  1:2276765/1‑100 (MQ=255)
                                                                            tACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATgcgc                        >  2:751802/1‑100 (MQ=255)
                                                                            tACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATgcgc                        <  2:548724/100‑1 (MQ=255)
                                                                              cacCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCg                      >  2:752476/1‑100 (MQ=255)
                                                                              cacCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCg                      >  1:988597/1‑100 (MQ=255)
                                                                                ccATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTg                    <  2:1895678/100‑1 (MQ=255)
                                                                                ccATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTg                    <  2:1730464/100‑1 (MQ=255)
                                                                                ccATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTg                    >  1:2494466/1‑100 (MQ=255)
                                                                                 cATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTgg                   <  1:1445185/100‑1 (MQ=255)
                                                                                 cATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTgg                   <  1:1826221/100‑1 (MQ=255)
                                                                                 cATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTgg                   <  1:725521/100‑1 (MQ=255)
                                                                                  atatAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGa                  >  1:1027104/1‑100 (MQ=255)
                                                                                    ataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAg                <  2:2155355/100‑1 (MQ=255)
                                                                                     taACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAgg               <  2:1183841/100‑1 (MQ=255)
                                                                                      ctctgACGCTGCCCAGCGTCTGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc              >  1:436129/5‑100 (MQ=255)
                                                                                      ctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc              >  2:2002285/5‑100 (MQ=255)
                                                                                      ctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc              >  1:67296/5‑100 (MQ=255)
                                                                                       tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa             >  1:2012904/4‑100 (MQ=255)
                                                                                       aCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa             >  1:7057/1‑100 (MQ=255)
                                                                                        ctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCAt            >  1:2330470/3‑100 (MQ=255)
                                                                                         tgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATa           >  1:2124802/2‑100 (MQ=255)
                                                                                         aGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATa           >  2:1089832/1‑100 (MQ=255)
                                                                                          gACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAg          >  1:127489/1‑100 (MQ=255)
                                                                                              cTGGGCAGCGTCAGAGAAATGTCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAgcgc       <  2:1832280/99‑1 (MQ=255)
                                                                                              cTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCATGCAATGCGCAGAGCACGGCATaaaaac      >  1:1108592/1‑95 (MQ=255)
                                                                                               tGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAc     >  2:286171/1‑100 (MQ=255)
                                                                                                gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc    <  2:2529996/100‑1 (MQ=255)
                                                                                                gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc    <  1:1156174/100‑1 (MQ=255)
                                                                                                gcgCAGCGTCAGATAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc    >  1:1725740/3‑100 (MQ=255)
                                                                                                  |                                                                                                   
GTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTG  >  NC_000913/3601727‑3601924

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: