Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,825 G→C 16.5% 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,8250GC16.5% 520.1 / 90.3 201P88A (CCC→GCC) ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (101/66);  new base C (16/17);  total (118/83)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

Only 100 of 212 total aligned reads displayed.

CGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTG  >  NC_000913/3601726‑3601924
                                                                                                   |                                                                                                   
     aGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCg                                                                                                <  2:4147294/100‑1 (MQ=255)
     aGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCg                                                                                                <  1:805903/100‑1 (MQ=255)
     aGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCg                                                                                                <  2:3172494/100‑1 (MQ=255)
      gCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGt                                                                                               <  1:1272872/100‑1 (MQ=255)
        gTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCCGCGTc                                                                                              >  1:993769/1‑99 (MQ=255)
         ttttttTCAGATAAACAGTGAATTGCGGTGACGGATACTACTGCGTCGCCGCCTCGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTct                                                                                             <  1:3553836/99‑2 (MQ=255)
          tttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCaga                                                                                           <  1:2321065/100‑1 (MQ=255)
          tttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgt                                                                                           <  2:1401614/100‑4 (MQ=255)
          tttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgt                                                                                           <  1:2624967/100‑4 (MQ=255)
          tttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgt                                                                                           <  2:4028528/100‑4 (MQ=255)
           ttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACTGACGCTGCCCAGCGTctgtt                                                                                          <  1:2384578/100‑5 (MQ=255)
           ttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGAACACCATATAACAGACGCTGCCCAGCGTctgtt                                                                                          <  1:671039/100‑5 (MQ=255)
            tttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCagag                                                                                          <  1:1303509/99‑1 (MQ=255)
             ttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGaa                                                                                        >  1:3312153/1‑100 (MQ=255)
              tGCAGATAATCAGTGATTTGCGGTGACGGATAATACTGCCTCGCCGCGTGGTCAACGTTTTTGTACACCATATAACAGACGCTGCGCAGCGTctgttata                                                                                       <  1:1540904/100‑8 (MQ=255)
               gCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAAt                                                                                      <  1:2863849/100‑1 (MQ=255)
               gCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAAt                                                                                      >  1:3388626/1‑100 (MQ=255)
                 aGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACTGACGCTGCCCAGCGTctgttatatgg                                                                                    <  2:1874600/100‑11 (MQ=37)
                 aGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATgg                                                                                    >  2:2494323/1‑100 (MQ=255)
                  gATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGc                                                                                   >  2:1095498/1‑100 (MQ=255)
                  gATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGc                                                                                   <  2:1883301/100‑1 (MQ=255)
                  gATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGc                                                                                   >  1:4396442/1‑100 (MQ=255)
                    tAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACCGACGCTGGGCAGAGTCAGAGAAATGGCg                                                                                  >  1:4180752/1‑99 (MQ=255)
                    tAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGa                                                                                 >  2:3332905/1‑100 (MQ=255)
                      aaCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATATCAGACGCTGGGCAGCGTCAGAGAAATGGCGAt                                                                                <  1:4328175/99‑1 (MQ=255)
                      aaCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATa                                                                               >  1:1115781/1‑100 (MQ=255)
                       aCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATaa                                                                              >  2:2572982/1‑100 (MQ=255)
                        cAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAc                                                                             <  2:1351947/100‑1 (MQ=255)
                        cAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAc                                                                             <  2:555288/100‑1 (MQ=255)
                        cAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAc                                                                             >  1:914463/1‑100 (MQ=255)
                          gTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCa                                                                           >  1:2174110/1‑100 (MQ=255)
                          gTGATTGGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATa                                                                               <  1:4110558/96‑1 (MQ=255)
                            gATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATc                                                                         <  1:1013895/100‑1 (MQ=255)
                              tttGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCAc                                                                       <  1:3252818/100‑1 (MQ=255)
                              tttGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCAc                                                                       >  1:1728804/1‑100 (MQ=255)
                               ttGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCAcc                                                                      <  1:587256/100‑1 (MQ=255)
                                tGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCg                                                                     >  2:2931349/1‑100 (MQ=255)
                                tGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGACAGCGTCAGAGAAATGGCGATAACCATCACCg                                                                     >  1:1880153/1‑100 (MQ=255)
                                  cGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGtt                                                                   <  1:3440093/100‑1 (MQ=255)
                                   ggTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTa                                                                  <  2:244497/100‑1 (MQ=255)
                                         ggATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACg                                                            >  2:271594/1‑100 (MQ=255)
                                            taataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTgg                                                         <  2:1833000/100‑1 (MQ=255)
                                             aataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCGTATAACAGACGCTGGGCAGCGTGAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGc                                                        <  1:2734167/100‑1 (MQ=255)
                                              ataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCg                                                       <  2:2543503/100‑1 (MQ=255)
                                              ataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCg                                                       >  1:2836705/1‑100 (MQ=255)
                                               taCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATGTAGCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCg                                                       <  1:3701962/99‑1 (MQ=255)
                                               taCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGa                                                      >  2:2121619/1‑100 (MQ=255)
                                                 cTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAAc                                                    <  2:4264194/100‑1 (MQ=255)
                                                   gCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACgg                                                  >  1:347863/1‑100 (MQ=255)
                                                   gCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACgg                                                  >  1:1293297/1‑100 (MQ=255)
                                                        gccgccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGc                                             <  2:3388644/100‑1 (MQ=255)
                                                           gccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCttt                                          >  1:1617412/1‑100 (MQ=255)
                                                             cTGGTTAAAGTTTTTGTACACCATATAACAGACGCTGAGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTc                                        <  1:361157/100‑1 (MQ=255)
                                                             atGGTTAACGTTTTTGTACACCATATAACAGACGCAGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGATTGCTTTTc                                        >  2:3064851/2‑100 (MQ=255)
                                                              tGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCa                                       <  1:4303521/100‑1 (MQ=255)
                                                                 ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGAt                                    >  2:2818756/1‑100 (MQ=255)
                                                                 ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGAt                                    >  1:3679174/1‑100 (MQ=255)
                                                                  tAACGTTTTTTTACACCATATAAAAGACGATGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGTTc                                   >  1:4014128/1‑100 (MQ=255)
                                                                  tAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAAAGTGGCGAACGGTTTGCTTTTCAGATc                                   >  2:164755/1‑100 (MQ=255)
                                                                     cGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTg                                >  2:2651288/1‑100 (MQ=255)
                                                                     cGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTg                                >  2:1527975/1‑100 (MQ=255)
                                                                      gTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGTTTTGCTTTTCAGATCCTGc                               <  2:3243667/100‑1 (MQ=255)
                                                                      gTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGc                               >  1:3765922/1‑100 (MQ=255)
                                                                       tttttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCa                              <  2:1301465/100‑1 (MQ=255)
                                                                        ttttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAGAAACGTGGCGAACGGTTTGCTTTTCAGAGCCTGCaa                             >  2:1458787/1‑100 (MQ=255)
                                                                        ttttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCaa                             >  2:1998966/1‑100 (MQ=255)
                                                                          gatctcacCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATg                           <  2:170892/95‑1 (MQ=255)
                                                                             tACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATgcgc                        >  2:789901/1‑100 (MQ=255)
                                                                                acCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGt                     <  2:809162/100‑1 (MQ=255)
                                                                                acCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGt                     >  2:3860289/1‑100 (MQ=255)
                                                                                acCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGt                     >  1:1682368/1‑100 (MQ=255)
                                                                                 ccATATAAAAGACGCTGGGCAGCGTCAGAGAACTGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCGTTTCAGAGCCTGCAATGCGCCg                      >  1:3038871/1‑98 (MQ=255)
                                                                                  catttctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTgg                   >  2:3582509/10‑100 (MQ=255)
                                                                                  cATATGACAGACGCTGGGCAGCGTCCGAGAAATGGCGATAACCATCACCGGTAAAAACGCGGAGAACGGTTTGCGTTGCAGATCATGCAATGCGCAGTgg                   >  1:2699236/1‑100 (MQ=255)
                                                                                      tctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAgg               >  2:2621503/6‑100 (MQ=255)
                                                                                      taACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAgg               >  1:3814275/1‑100 (MQ=255)
                                                                                      taACAGACGCTGAGGAGCGTAAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCCTTTCAGACACTGCAATGCGCCGTGGACgg               <  1:3508455/100‑1 (MQ=255)
                                                                                       ctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc              >  2:1415248/5‑100 (MQ=255)
                                                                                       ctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc              >  1:3224092/5‑100 (MQ=255)
                                                                                       aaCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc              >  2:3956008/1‑100 (MQ=255)
                                                                                       aaCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGATTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc              >  2:700003/1‑100 (MQ=255)
                                                                                        tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa             >  2:3499966/4‑100 (MQ=255)
                                                                                        tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa             >  2:2762598/4‑100 (MQ=255)
                                                                                        tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa             >  2:523338/4‑100 (MQ=255)
                                                                                        tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa             >  2:710838/4‑100 (MQ=255)
                                                                                         ctgACGCTGCCCAGCGTCAGAGAAATGGCGATGACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGACGTGGAAGGCa             >  1:1629054/3‑99 (MQ=255)
                                                                                         cAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCAt            >  2:2202986/1‑100 (MQ=255)
                                                                                           gACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAg          >  1:2490291/1‑100 (MQ=255)
                                                                                            aCGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAgc         <  2:3529884/100‑1 (MQ=255)
                                                                                             cGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTACAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGACGGCATAgcg        >  1:155338/1‑100 (MQ=255)
                                                                                             cGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAgcg        >  1:1823232/1‑100 (MQ=255)
                                                                                               cTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCa      >  1:258258/1‑100 (MQ=255)
                                                                                                tGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAc     <  1:3113650/100‑1 (MQ=255)
                                                                                                tGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAc     >  2:177958/1‑100 (MQ=255)
                                                                                                 gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc    <  2:3765946/100‑1 (MQ=255)
                                                                                                 gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc    <  2:3104824/100‑1 (MQ=255)
                                                                                                 gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc    <  1:2121611/100‑1 (MQ=255)
                                                                                                  ggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCt   >  1:2069026/1‑100 (MQ=255)
                                                                                                  ggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCt   >  2:2384226/1‑100 (MQ=255)
                                                                                                  ggCAGCGTCAGAGAAAAGGCGAAAACCATCACCGTTAAAAACGTGGCGAACGGTTTGATTTTCAAATCCTGCAACGCGCCGCGGAAGGCATAGCGCAcc    <  1:3953451/99‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTG  >  NC_000913/3601726‑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: