Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,824 G→C 17.1% P88R (CCC→CGC) ‡ ftsX ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,8240GC17.1% 504.6 / 93.4 201P88R (CCC→CGC) ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (102/63);  new base C (16/18);  total (118/83)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 214 total aligned reads displayed.

TCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCT  >  NC_000913/3601725‑3601923
                                                                                                   |                                                                                                   
tCGTCCAGCGTTTTTGGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTgg                                                                                                     <  2:4332793/100‑1 (MQ=255)
 cGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGCTGAACGTTTTTGTACACCATATAACAGACCCTggg                                                                                                    >  1:2375328/1‑100 (MQ=255)
  gTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGc                                                                                                   <  1:3409611/100‑1 (MQ=255)
  gTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGc                                                                                                   >  2:3581010/1‑100 (MQ=255)
    ccAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAg                                                                                                 >  1:3978120/1‑100 (MQ=255)
      aGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCg                                                                                               <  2:4147294/100‑1 (MQ=255)
      aGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCg                                                                                               <  2:3172494/100‑1 (MQ=255)
       gCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGt                                                                                              <  1:1272872/100‑1 (MQ=255)
          ttttttTCAGATAAACAGTGAATTGCGGTGACGGATACTACTGCGTCGCCGCCTCGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTct                                                                                            <  1:3553836/99‑2 (MQ=255)
          ttttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCa                                                                                            >  1:2844706/1‑99 (MQ=255)
           tttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgt                                                                                          <  1:2624967/100‑4 (MQ=255)
            ttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCagag                                                                                         <  1:3811727/100‑1 (MQ=255)
            ttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCagag                                                                                         <  2:1582612/100‑1 (MQ=255)
             tttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCagag                                                                                         <  1:1303509/99‑1 (MQ=255)
             tttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgtta                                                                                        <  1:972467/100‑6 (MQ=255)
              ttGCAGATAAACAGTGATTTGCGGTGACGGTTAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGaa                                                                                       >  2:2486430/1‑100 (MQ=255)
               tGCAGATAATCAGTGATTTGCGGTGACGGATAATACTGCCTCGCCGCGTGGTCAACGTTTTTGTACACCATATAACAGACGCTGCGCAGCGTctgttata                                                                                      <  1:1540904/100‑8 (MQ=255)
               tGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgttata                                                                                      <  2:2813070/100‑8 (MQ=255)
               tGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCAGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGaaa                                                                                      >  2:2159704/1‑100 (MQ=255)
                gCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAAt                                                                                     >  2:2544010/1‑100 (MQ=255)
                  aGATAAGCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACCGACGCTGCCCAGCGTctgttatatgg                                                                                   <  2:1874716/100‑11 (MQ=38)
                  aGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATgg                                                                                   >  2:2494767/1‑100 (MQ=255)
                     tAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGa                                                                                >  2:3332905/1‑100 (MQ=255)
                      aaaCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGAt                                                                               <  2:1035421/100‑1 (MQ=255)
                       aaCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATa                                                                              >  1:1115781/1‑100 (MQ=255)
                        aCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATaa                                                                             >  2:1411808/1‑100 (MQ=255)
                         cAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAc                                                                            <  2:555288/100‑1 (MQ=255)
                         cAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAc                                                                            >  2:3974275/1‑100 (MQ=255)
                          aGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAcc                                                                           >  2:1666393/1‑100 (MQ=255)
                           gTGATTGGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATa                                                                              <  1:4110558/96‑1 (MQ=255)
                             gATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATc                                                                        <  1:1013895/100‑1 (MQ=255)
                              aTTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCa                                                                       >  2:2722259/1‑100 (MQ=255)
                               tttGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCAc                                                                      <  1:3252818/100‑1 (MQ=255)
                               tttGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCAc                                                                      >  1:1728804/1‑100 (MQ=255)
                                 tGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCg                                                                    >  2:1476428/1‑100 (MQ=255)
                                 tGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGACAGCGTCAGAGAAATGGCGATAACCATCACCg                                                                    >  1:1880153/1‑100 (MQ=255)
                                  gCGGTGACGGATCATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGt                                                                   >  2:2895377/1‑100 (MQ=255)
                                       gACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTaaaa                                                              <  2:4108826/100‑1 (MQ=255)
                                         cGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAAc                                                            <  2:1389559/100‑1 (MQ=255)
                                           gATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCCCCGTTAAAAACGt                                                          >  1:1634639/1‑100 (MQ=255)
                                           gATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGt                                                          <  1:3722756/100‑1 (MQ=255)
                                            ataataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTg                                                         <  1:2829277/100‑1 (MQ=255)
                                               ataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCg                                                      >  2:3358731/1‑100 (MQ=255)
                                                taCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGa                                                     >  2:2121619/1‑100 (MQ=255)
                                                  cTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAAc                                                   >  1:3511385/1‑100 (MQ=255)
                                                    gCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACgg                                                 >  1:1293297/1‑100 (MQ=255)
                                                     cGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGt                                                >  2:3879374/1‑100 (MQ=255)
                                                         gccgccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGc                                            <  2:3388644/100‑1 (MQ=255)
                                                         gccgccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGc                                            >  1:2324166/1‑100 (MQ=255)
                                                            gccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCttt                                         >  1:602760/1‑100 (MQ=255)
                                                            gccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCttt                                         >  2:3768012/1‑100 (MQ=255)
                                                            gccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCttt                                         >  1:1617412/1‑100 (MQ=255)
                                                              cTGGTTAAAGTTTTTGTACACCATATAACAGACGCTGAGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTc                                       <  1:361157/100‑1 (MQ=255)
                                                                ggTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAg                                     >  2:2948984/1‑100 (MQ=255)
                                                                  ttAGCGTTTTTGTACACCATATAACAGACGCTGGGAATCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGGTTTTCAGAt                                   <  1:896226/100‑1 (MQ=255)
                                                                  ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGa                                    >  2:1953632/1‑99 (MQ=255)
                                                                     aCGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCt                                >  1:3867431/1‑100 (MQ=255)
                                                                     aCGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACCGTTTGCTTTTCAGATCCt                                >  1:4324246/1‑100 (MQ=255)
                                                                      cGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTg                               >  2:3615760/1‑100 (MQ=255)
                                                                      cGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTg                               >  2:2651288/1‑100 (MQ=255)
                                                                       gTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGTTTTGCTTTTCAGATCCTGc                              <  2:3243667/100‑1 (MQ=255)
                                                                       gTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGc                              <  2:2385401/100‑1 (MQ=255)
                                                                        tttttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCa                             >  1:3698281/1‑100 (MQ=255)
                                                                        tttttGTAAACCATAAAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCa                             <  1:2257141/100‑1 (MQ=255)
                                                                              tACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATgcgc                       >  1:473474/1‑100 (MQ=255)
                                                                              tACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATgcgc                       >  2:789901/1‑100 (MQ=255)
                                                                               acacCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGcc                      >  2:1744190/1‑100 (MQ=255)
                                                                                cacCATATACCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTGTGCTTTTCAGATCCTGCAATGCGCCg                     <  2:2681854/100‑1 (MQ=255)
                                                                                 acCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGt                    <  2:809162/100‑1 (MQ=255)
                                                                                 acCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGt                    >  1:1682368/1‑100 (MQ=255)
                                                                                  ccATATAAAAGACGCTGGGCAGCGTCAGAGAACTGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCGTTTCAGAGCCTGCAATGCGCCg                     >  1:3038871/1‑98 (MQ=255)
                                                                                   cATATGACAGACGCTGGGCAGCGTCCGAGAAATGGCGATAACCATCACCGGTAAAAACGCGGAGAACGGTTTGCGTTGCAGATCATGCAATGCGCAGTgg                  >  1:2699236/1‑100 (MQ=255)
                                                                                   cATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAGCCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTg                   <  1:2908217/99‑1 (MQ=255)
                                                                                    atatAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGa                 <  1:403635/100‑1 (MQ=255)
                                                                                       tctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAgg              >  2:2621503/6‑100 (MQ=255)
                                                                                       taACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAgg              >  1:34428/1‑100 (MQ=255)
                                                                                       taACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAgg              >  1:4211677/1‑100 (MQ=255)
                                                                                        ctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc             >  2:1795438/5‑100 (MQ=255)
                                                                                        ctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc             >  1:3224092/5‑100 (MQ=255)
                                                                                        ctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc             >  2:1416313/5‑100 (MQ=255)
                                                                                        aaCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGATTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc             >  2:700003/1‑100 (MQ=255)
                                                                                         tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa            >  1:4352163/4‑100 (MQ=255)
                                                                                         tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa            >  2:523338/4‑100 (MQ=255)
                                                                                         tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa            >  2:710838/4‑100 (MQ=255)
                                                                                         tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa            >  2:112738/4‑100 (MQ=255)
                                                                                         tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa            >  2:1324277/4‑100 (MQ=255)
                                                                                         tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCa            >  2:3499966/4‑100 (MQ=255)
                                                                                           tgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATa          <  2:1857276/99‑1 (MQ=255)
                                                                                            gACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAg         >  1:2490291/1‑100 (MQ=255)
                                                                                              cGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTACAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGACGGCATAgcg       >  1:155338/1‑100 (MQ=255)
                                                                                              cGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAgcg       >  1:1823232/1‑100 (MQ=255)
                                                                                              cGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAgcg       >  1:77973/1‑100 (MQ=255)
                                                                                               gCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGTCGAACGGTTTGCTTTTCAGCTCCTGCAATGAGCCGTCGAGGCCATa          >  1:3357277/1‑96 (MQ=255)
                                                                                                cTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCa     >  1:258258/1‑100 (MQ=255)
                                                                                                 tGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAc    <  1:3113650/100‑1 (MQ=255)
                                                                                                  gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc   <  2:3104824/100‑1 (MQ=255)
                                                                                                  gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc   <  2:2218902/100‑1 (MQ=255)
                                                                                                  gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc   <  1:2121611/100‑1 (MQ=255)
                                                                                                  gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc   <  1:3498073/100‑1 (MQ=255)
                                                                                                   ggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCt  >  1:2069026/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCT  >  NC_000913/3601725‑3601923

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: