Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,122,984 C→T R119H (CGT→CAT)  ftsN ← essential cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,122,9840CT100.0% 200.5 / NA 63R119H (CGT→CAT) ftsNessential cell division protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/37);  total (26/37)

CTGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGT  >  NC_000913_3_bme_pgi/4122845‑4123129
                                                                                                                                           |                                                                                                                                                 
cTGTTCCCCTCATTGCTGCACCTGACGTTGGCCCTGTAGCGTTTGCTGACGCTATTCCGCCGTCTGTTCATTCCACGCAACTTCAACCACCTGCATTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTg                                                                                                                                          >  2:232765/1‑149 (MQ=255)
  gTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGgt                                                                                                                                        <  1:952604/149‑1 (MQ=255)
  gTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGgt                                                                                                                                        <  2:197924/149‑1 (MQ=255)
       gCTAATTTCTGCGCCGGAGGTTGGGGCTGTAGCGTTTGCTCACCCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCGTTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAg                                                                                                                                   <  2:216909/149‑1 (MQ=255)
       gCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAg                                                                                                                                   <  1:129974/149‑1 (MQ=255)
               ctgcGCCTGACGTTGGCGCTGTAGCGTTTGCTGAGGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCCGTTGCTcct                                                                                                                           >  2:64589/1‑148 (MQ=255)
                 gcgcCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGCACTTCAACCAGCTGCGTTGTCCGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGc                                                                                                                         >  1:237230/1‑149 (MQ=255)
                   gcCTGACGTTGGCGCTGTAGCGGTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGt                                                                                                                       <  1:490080/149‑1 (MQ=255)
                      tGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGtttt                                                                                                                    >  2:125667/1‑149 (MQ=255)
                      tGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTCGTTCTGGTGTCAGTTGCTCCGGCGtttt                                                                                                                    >  2:64550/1‑149 (MQ=255)
                        aCGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttca                                                                                                                  >  1:767509/1‑149 (MQ=255)
                          gTTTGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCCGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTTGTGTCAGTTGCTCCGGCGTTttcact                                                                                                                <  1:771793/149‑1 (MQ=255)
                           ttGGGGGTGTAGCGTTTGCTGACGTTGTTCGGGGGTCTGTTCATTTCACGGAACTTCAACGAGCGGCGTTGGCTGTTCTCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCACGGGCTCCGGCGATttctact                                                                                                               <  2:809973/149‑5 (MQ=255)
                            tGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttcacttc                                                                                                              <  1:285122/149‑1 (MQ=255)
                              gcgcTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTTATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttcacttcac                                                                                                            <  1:151790/149‑1 (MQ=255)
                              gcgcTGTAGCGGTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTACTGGCGCATATCAGCCTGCATTTGTTCCAGCTGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttcacttcac                                                                                                            <  2:5599/149‑1 (MQ=255)
                                 cTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTcaccac                                                                                                         >  2:465731/1‑149 (MQ=255)
                                   gTAGCGGTTGCGGACGCTGGTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGGTGCTGGCGCAGATCAGCCTGCATTTGTTCAAGAAGCCGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCg                                                                                                       <  2:197965/149‑1 (MQ=255)
                                        gTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGa                                                                                                  >  1:407944/1‑149 (MQ=255)
                                               gACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCt                                                                                           <  2:947113/149‑1 (MQ=255)
                                               gACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCt                                                                                           <  1:172435/149‑1 (MQ=255)
                                                 cGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCtgt                                                                                         <  1:316463/149‑1 (MQ=255)
                                                  gCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTGTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCtgtg                                                                                        >  1:533751/1‑149 (MQ=255)
                                                   gtGTTCCGGCGTCTGTTCATTCCCCGGAACTTCAACCAGCTGCGTTGGCTGCTAGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTgg                                                                                       <  2:766903/148‑1 (MQ=255)
                                                   cTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGACTCCATTTGTGCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAGGTTTCTGTgg                                                                                       >  2:416432/1‑149 (MQ=255)
                                                        ccGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGcgca                                                                                  >  2:219466/1‑149 (MQ=255)
                                                        ccGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGcgca                                                                                  >  1:842854/1‑149 (MQ=255)
                                                                tGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTcc                                                                          <  2:237230/149‑1 (MQ=255)
                                                                 gTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGGTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTccc                                                                         >  2:812278/1‑149 (MQ=255)
                                                                     aTTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCt                                                                     <  1:373071/149‑1 (MQ=255)
                                                                      ttCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTg                                                                    <  1:235760/149‑1 (MQ=255)
                                                                        ccACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTCGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGc                                                                  <  1:64550/149‑1 (MQ=255)
                                                                          actgAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGa                                                                <  1:424512/146‑1 (MQ=255)
                                                                                    aaCCAGCTGCGTTGGCTGCTGGCGAATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGct                                                      >  1:654433/1‑149 (MQ=255)
                                                                                     aCCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGctc                                                     >  1:198792/1‑149 (MQ=255)
                                                                                      ccAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGctct                                                    <  1:315279/149‑1 (MQ=255)
                                                                                      ccAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGctct                                                    >  1:303083/1‑149 (MQ=255)
                                                                                            gCGTTGGCTGCTGGCGCATATCGGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATg                                              <  1:982888/149‑1 (MQ=255)
                                                                                              gTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTTTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTa                                            >  1:450180/1‑149 (MQ=255)
                                                                                                 ggCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGAGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTTgcg                                         >  2:419806/1‑149 (MQ=255)
                                                                                                   ctgctgTCGCATATCCGGCTGCATTTGTAAAAGAATCTGGTGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGACGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGcc                                       <  2:381677/149‑1 (MQ=255)
                                                                                                    tgctgGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCa                                      >  2:649680/1‑149 (MQ=255)
                                                                                                      ctgGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGc                                    <  1:381247/149‑1 (MQ=255)
                                                                                                      ctgGCGCATATCAGCCGGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGc                                    <  1:823802/149‑1 (MQ=255)
                                                                                                         gcgcATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGtt                                 >  1:124704/1‑149 (MQ=255)
                                                                                                           gcATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGttct                               >  1:455977/1‑149 (MQ=255)
                                                                                                            cATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGttctt                              <  2:842854/149‑1 (MQ=255)
                                                                                                            cATATCAGCCGGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCTGCGttctg                              <  2:10742/149‑2 (MQ=255)
                                                                                                              tatCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTtcttct                            <  2:366424/149‑1 (MQ=255)
                                                                                                               atCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTg                           <  1:465731/149‑1 (MQ=255)
                                                                                                               atCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTg                           >  2:726299/1‑149 (MQ=255)
                                                                                                                tCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTgg                          <  1:939858/149‑1 (MQ=255)
                                                                                                                 cAGCCTGCATTTGTTCAGGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGt                         <  2:654433/149‑1 (MQ=255)
                                                                                                                     cTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtg                     <  1:705018/149‑1 (MQ=255)
                                                                                                                      tGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtgg                    >  1:367541/1‑149 (MQ=255)
                                                                                                                       gCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggt                   <  1:649680/149‑1 (MQ=255)
                                                                                                                        taTTTGTTGAAGAAGCGGATGTTGTTCTGGTGTCAGGTGCTCCGGCGTTTTCACGTCACCACCGGCAGAAGGGTCTGTGGGCGCACGCGCTCCCGGCGGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggtg                  <  2:474721/148‑1 (MQ=255)
                                                                                                                         aTTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggtgg                 >  2:331945/1‑149 (MQ=255)
                                                                                                                         aTTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggtgg                 <  1:219466/149‑1 (MQ=255)
                                                                                                                         aTTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggtgg                 <  2:997037/149‑1 (MQ=255)
                                                                                                                              ttCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTTGTTTTGGTGGTAGTc            <  2:373349/149‑1 (MQ=255)
                                                                                                                              ttCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGAACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTGTTCTGGTTTTGGTGGTAGTc            >  1:742012/1‑149 (MQ=255)
                                                                                                                                        tGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGt  >  1:727965/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
CTGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGT  >  NC_000913_3_bme_pgi/4122845‑4123129

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