Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,642,597 A→G 6.0% L83P (CTG→CCG)  rodZ ← MreB assembly cytoskeletal protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,642,5970AG 6.0% 217.0 / 4.3 67L83P (CTG→CCG) rodZMreB assembly cytoskeletal protein
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base G (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

ACCGCTCAGGCCGATAACCACAAACAACACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAG  >  NC_000913/2642458‑2642736
                                                                                                                                           |                                                                                                                                           
cccGCTCAGCCCGAGACCCCCACCCACCACCCGCCACGTGAAGGTCATCTGCCAGCCGTCGCGTTATTTGCGGCGTTTACCGAGGGAAAAACTCTGCAACGGCGCAACTTTAGCAGCCCGAAGTGGAGCCTGCTATTCCa                                                                                                                                             <  2:1511082/139‑1 (MQ=255)
aCCGCTCAGGCCGATAACCACAACCAACACCAGCCAAGTTAAGGTCAACTGCCAGCCGTCGCCTTTAGTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCATTTCCa                                                                                                                                             <  2:1511014/140‑1 (MQ=255)
  cGCTCAGGCCGATAACCACAAACAACACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGc                                                                                                                                           >  1:75728/1‑140 (MQ=255)
         gCCGATAACCACAAACAACACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGgca                                                                                                                                    >  2:916117/1‑140 (MQ=255)
         gCCGATAACCACAAACAACACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGgca                                                                                                                                    >  2:915967/1‑140 (MQ=255)
           cGATAACCACAAACACCACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGgcagc                                                                                                                                  <  2:1149117/140‑1 (MQ=255)
           cGATAACCACAAACAACACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGgcagc                                                                                                                                  <  2:742525/140‑1 (MQ=255)
            gATAACCACAAACAACACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGgcagca                                                                                                                                 >  2:273921/1‑140 (MQ=255)
              tAACCACAAACAACACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGt                                                                                                                               >  2:319777/1‑140 (MQ=255)
              tAACCACAAACAACACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGt                                                                                                                               >  1:400259/1‑140 (MQ=255)
                   acaatcaaCACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGttcttc                                                                                                                          >  1:609597/1‑140 (MQ=255)
                   acaaacaaCACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGttcttc                                                                                                                          >  1:1315753/1‑140 (MQ=255)
                   acaaacaaCACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGttcttc                                                                                                                          >  1:383532/1‑140 (MQ=255)
                   acaaacaaCACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGttcttc                                                                                                                          >  1:382756/1‑140 (MQ=255)
                       acaacaCCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTtcttcttct                                                                                                                      <  2:406179/140‑1 (MQ=255)
                        caacaCCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTg                                                                                                                     <  1:107955/140‑1 (MQ=255)
                            acCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGGTCTTCTTCTGGAat                                                                                                                 >  1:1336635/1‑140 (MQ=255)
                             ccAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAata                                                                                                                <  1:1143067/140‑1 (MQ=255)
                             ccAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAata                                                                                                                <  1:1142224/140‑1 (MQ=255)
                             ccAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAata                                                                                                                <  2:230654/140‑1 (MQ=255)
                                  cAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTAc                                                                                                           <  1:601291/140‑1 (MQ=255)
                                   aaGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTAcc                                                                                                          >  1:466589/1‑140 (MQ=255)
                                   aaGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTAcc                                                                                                          <  1:1020755/140‑1 (MQ=255)
                                   aaGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTAcc                                                                                                          >  1:267989/1‑140 (MQ=255)
                                     gTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAg                                                                                                        <  1:916117/140‑1 (MQ=255)
                                     gTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAg                                                                                                        <  1:915967/140‑1 (MQ=255)
                                       gAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGAc                                                                                                      <  1:835850/140‑1 (MQ=255)
                                          ggTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGAcgcg                                                                                                   >  1:864515/1‑140 (MQ=255)
                                           gTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCATCTTCTGGAATATGTACCAGACgcgc                                                                                                  >  2:559094/1‑140 (MQ=255)
                                           gTCATCAGCCAGCCGGCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACgcgc                                                                                                  >  2:559045/1‑140 (MQ=255)
                                             catcaGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCAt                                                                                                <  2:36972/140‑1 (MQ=255)
                                                 agccagccGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAaga                                                                                            <  2:869011/140‑1 (MQ=255)
                                                              gTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCcg                                                                               <  2:387874/140‑1 (MQ=255)
                                                                 ttttGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCgcgc                                                                            >  1:1588869/1‑140 (MQ=255)
                                                                 ttttGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACGTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCgcgc                                                                            >  2:721763/1‑140 (MQ=255)
                                                                     gcggcgTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGa                                                                        <  1:1512658/140‑1 (MQ=255)
                                                                     gcggcgTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGa                                                                        <  1:1572334/140‑1 (MQ=255)
                                                                            ttACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCGGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGa                                                                 <  2:889679/140‑1 (MQ=255)
                                                                            ttACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCGGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGa                                                                 <  2:889702/140‑1 (MQ=255)
                                                                                   gggAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGa                                                          >  2:864815/1‑140 (MQ=255)
                                                                                    ggAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGAt                                                         >  2:158713/1‑140 (MQ=255)
                                                                                       aaaaCTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATcgg                                                      <  2:351328/140‑1 (MQ=255)
                                                                                       aaaaCTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATcgg                                                      <  2:252057/140‑1 (MQ=255)
                                                                                       aaaaCTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATcgg                                                      <  2:464067/140‑1 (MQ=255)
                                                                                       aaaaCACTGCAACGCCGCAACTATTGAAGCCCGAAGTGGACCAAGCATTTCCAGCCCTCGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATcgg                                                      <  2:255069/140‑1 (MQ=255)
                                                                                          aCTCTGCATCGGCGCAACTTTTGCAGCCCGAGGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATcggcgg                                                   <  1:63113/140‑1 (MQ=255)
                                                                                          aCTCTGCATCGGCGCAACTATTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGACTATGTACCAGACGCGCATAAGATCGGATATATCCGCGCAGGCATGTTGAAGCAAGATcggcgg                                                   <  2:1262925/140‑1 (MQ=255)
                                                                                               gCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGcc                                              <  2:1011953/140‑1 (MQ=255)
                                                                                               gCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGcc                                              <  1:1075184/140‑1 (MQ=255)
                                                                                                    ggCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTAtc                                         >  2:1052737/1‑140 (MQ=255)
                                                                                                    ggCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTAtc                                         <  2:967945/140‑1 (MQ=255)
                                                                                                    ggCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTAtc                                         <  1:129743/140‑1 (MQ=255)
                                                                                                    ggCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTAtc                                         >  2:1051770/1‑140 (MQ=255)
                                                                                                        cAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCGGCCCTGGCAGCCGTTCTTCTTCTGGAATATGTACCAGACGCGCATACGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTAtctact                                     >  2:1691618/1‑140 (MQ=255)
                                                                                                        cAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCGGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATCTATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTGAtcttct                                     >  2:889517/1‑140 (MQ=255)
                                                                                                           cTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCa                                  <  1:1298673/140‑1 (MQ=255)
                                                                                                           cTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCa                                  <  1:1298587/140‑1 (MQ=255)
                                                                                                           cTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCa                                  <  1:1451599/140‑1 (MQ=255)
                                                                                                           cTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCa                                  <  1:1451660/140‑1 (MQ=255)
                                                                                                                   gCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTcgcg                          <  1:1229447/140‑1 (MQ=255)
                                                                                                                   gCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTcgcg                          <  1:1229258/140‑1 (MQ=255)
                                                                                                                   gCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAAGGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTcgcg                          <  1:1229237/140‑1 (MQ=255)
                                                                                                                     ccGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTa                        >  1:839196/1‑140 (MQ=255)
                                                                                                                                  tGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAATAGCGGATATATCCGCGCAGTAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCtt           >  2:1633309/1‑140 (MQ=255)
                                                                                                                                  tGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTCTCTTCTTCAATGTCGCGTACCGTGGAAACCtt           >  1:1001587/1‑140 (MQ=255)
                                                                                                                                  tGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATAGATCCGCGCAGGAATGTTGAAGCAAGATCGGAGGGTGCGTTGTCGTCTACAATGTCGCGTACCGGGGAAACCtt           >  2:1633025/1‑140 (MQ=255)
                                                                                                                                    cTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTCCCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCa         >  2:40617/1‑140 (MQ=255)
                                                                                                                                       ttCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGc      >  1:1006904/1‑140 (MQ=255)
                                                                                                                                          cAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAAGGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCaaa   >  2:202346/1‑140 (MQ=255)
                                                                                                                                           aGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGTGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAg  <  2:92972/140‑1 (MQ=255)
                                                                                                                                           aGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCTAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAg  <  2:92876/140‑1 (MQ=255)
                                                                                                                                           aGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAg  <  2:92858/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
ACCGCTCAGGCCGATAACCACAAACAACACCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAG  >  NC_000913/2642458‑2642736

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: