Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,634,642 C→T G82S (GGC→AGC)  tfaQ ← Qin prophage; putative tail fiber assembly protein TfaQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,634,6420CT100.0% 195.5 / NA 86G82S (GGC→AGC) tfaQQin prophage; putative tail fiber assembly protein TfaQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/42);  total (44/42)

CGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACC  >  NC_000913/1634506‑1634777
                                                                                                                                        |                                                                                                                                       
cGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCATTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGa                                                                                                                                      >  2:596416/1‑140 (MQ=17)
 ggATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGaa                                                                                                                                     >  2:679344/1‑140 (MQ=17)
        gAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGt                                                                                                                              >  2:309100/1‑140 (MQ=17)
              tttttCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATaaa                                                                                                                        <  1:309100/140‑1 (MQ=18)
                 ttCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCg                                                                                                                     <  2:387213/140‑1 (MQ=18)
                         ccGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCATTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTAccc                                                                                                             <  1:596416/140‑1 (MQ=16)
                                  tCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCt                                                                                                    <  2:52656/140‑1 (MQ=18)
                                   cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAAGGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCtt                                                                                                   >  2:179967/1‑140 (MQ=17)
                                   cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTAccc                                                                                                             <  2:385808/130‑1 (MQ=18)
                                   cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTAccc                                                                                                             >  1:385808/1‑130 (MQ=17)
                                   cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCtt                                                                                                   >  1:893221/1‑140 (MQ=17)
                                         ggCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgat                                                                                                         >  2:497217/1‑128 (MQ=17)
                                         ggCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgat                                                                                                         <  1:497217/128‑1 (MQ=18)
                                         ggCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACga                                                                                             <  2:274634/140‑1 (MQ=17)
                                          gCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACga                                                                                             <  2:832928/139‑1 (MQ=17)
                                          gCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACga                                                                                             >  1:832928/1‑139 (MQ=17)
                                            tgtgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGAt                                                                                          <  1:344682/140‑1 (MQ=17)
                                              tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGc                                                                                        >  2:232881/1‑140 (MQ=17)
                                              tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGc                                                                                        >  2:320684/1‑140 (MQ=17)
                                              tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGc                                                                                        >  1:518549/1‑140 (MQ=17)
                                               gCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAAGGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGcc                                                                                       <  1:179967/140‑1 (MQ=17)
                                                ccGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCa                                                                                      >  1:192863/1‑140 (MQ=17)
                                                ccGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCa                                                                                      >  1:351689/1‑140 (MQ=17)
                                                ccGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCa                                                                                      >  1:119986/1‑140 (MQ=17)
                                                   ttCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGa                                                                                   <  1:320684/140‑1 (MQ=17)
                                                      cACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGa                                                                                <  2:518549/140‑1 (MQ=17)
                                                        cTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACt                                                                              >  1:689974/1‑140 (MQ=17)
                                                           cTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat                                                                           >  1:475613/1‑140 (MQ=17)
                                                             gAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcatca                                                                         <  1:232881/140‑1 (MQ=17)
                                                              aaaCTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCAc                                                                        >  2:307032/1‑140 (MQ=17)
                                                              aaaCTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCAc                                                                        >  1:109130/1‑140 (MQ=17)
                                                                   ccccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTg                                                                   >  1:613987/1‑140 (MQ=17)
                                                                    cccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGa                                                                                   >  1:7560/1‑123 (MQ=17)
                                                                    cccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGa                                                                                   <  2:7560/123‑1 (MQ=17)
                                                                       ttCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGaaaa                                                               <  2:689974/140‑1 (MQ=17)
                                                                       ttCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGaaaa                                                               <  1:821392/140‑1 (MQ=17)
                                                                            gCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGcc                                                          >  2:661597/1‑140 (MQ=17)
                                                                            gCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGcc                                                          >  1:129063/1‑140 (MQ=17)
                                                                             cGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCa                                                         >  2:110987/1‑140 (MQ=17)
                                                                              gATAAGCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCAc                                                        >  2:419584/1‑140 (MQ=17)
                                                                              gATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCAc                                                        >  2:333673/1‑140 (MQ=17)
                                                                                     aGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagcca                                                 <  1:333673/140‑1 (MQ=17)
                                                                                      ggTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagccag                                                >  2:58803/1‑140 (MQ=17)
                                                                                      ggTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagccag                                                >  2:268500/1‑140 (MQ=17)
                                                                                       gTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagccagc                                               <  1:479754/140‑1 (MQ=17)
                                                                                           cATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagcca                                                 >  1:867405/1‑134 (MQ=17)
                                                                                           cATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagcca                                                 <  2:867405/134‑1 (MQ=17)
                                                                                           cATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGa                                           <  1:661597/140‑1 (MQ=17)
                                                                                           cATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGa                                           >  1:446336/1‑140 (MQ=17)
                                                                                            aTTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGaa                                          >  2:430637/1‑140 (MQ=17)
                                                                                                tCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCTGGAataag                                      >  1:33060/1‑138 (MQ=12)
                                                                                                tCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAatat                                       >  2:629161/1‑139 (MQ=17)
                                                                                                tCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAatat                                       <  1:629161/139‑1 (MQ=17)
                                                                                                 ccGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCt                                     >  1:863411/1‑140 (MQ=17)
                                                                                                 ccGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCt                                     >  2:736076/1‑140 (MQ=17)
                                                                                                 ccGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCt                                     >  2:829170/1‑140 (MQ=17)
                                                                                                  cGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCTGGAataaggg                                    <  2:33060/140‑5 (MQ=255)
                                                                                                   ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                                   <  2:119986/140‑1 (MQ=17)
                                                                                                   ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                                   <  2:192863/140‑1 (MQ=17)
                                                                                                   ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                                   <  2:637463/140‑1 (MQ=17)
                                                                                                   ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                                   <  1:268500/140‑1 (MQ=17)
                                                                                                     tAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                                   >  1:98937/1‑138 (MQ=18)
                                                                                                     tAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                                   <  2:98937/138‑1 (MQ=17)
                                                                                                      aaCGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCgg                                >  1:857370/1‑140 (MQ=18)
                                                                                                         ggACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGc                             <  2:446336/140‑1 (MQ=17)
                                                                                                         ggACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGc                             <  2:109130/140‑1 (MQ=17)
                                                                                                         gagacaggTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGc                             <  2:374025/133‑1 (MQ=12)
                                                                                                           aCCGAGTTCAGAAATAAATAACGCGTCGCTGGGAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCaa                           <  2:771662/140‑1 (MQ=17)
                                                                                                             cGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAata                         <  1:307032/140‑1 (MQ=17)
                                                                                                               aGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATc                       >  2:917721/1‑140 (MQ=17)
                                                                                                                gTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCg                      <  1:530174/140‑1 (MQ=17)
                                                                                                                 ttCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCgg                     <  2:129063/140‑1 (MQ=17)
                                                                                                                  tCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTa                   <  1:679344/140‑1 (MQ=12)
                                                                                                                         taaataaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGttt             <  1:874869/140‑1 (MQ=17)
                                                                                                                         taaataaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGttt             <  1:430637/140‑1 (MQ=17)
                                                                                                                           aataaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATACCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGc           >  1:186106/1‑140 (MQ=16)
                                                                                                                                cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTAc                  <  2:325836/128‑1 (MQ=17)
                                                                                                                                cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTAc                  >  1:325836/1‑128 (MQ=18)
                                                                                                                                cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTgg        <  2:124884/138‑1 (MQ=17)
                                                                                                                                cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTgg        >  1:124884/1‑138 (MQ=17)
                                                                                                                                cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa      <  1:395955/140‑1 (MQ=17)
                                                                                                                                cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa      <  1:419584/140‑1 (MQ=17)
                                                                                                                                cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa      <  1:829170/140‑1 (MQ=17)
                                                                                                                                cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa      >  2:127923/1‑140 (MQ=18)
                                                                                                                                  cgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGa    >  1:119392/1‑140 (MQ=17)
                                                                                                                                    tCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAGGCCAGCCGGAATATCTGGCGGTGCAATATCGGCACTGTTTGCTGGCAGAcc  >  2:791998/1‑140 (MQ=9)
                                                                                                                                        |                                                                                                                                       
CGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACC  >  NC_000913/1634506‑1634777

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