Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 393,298 T→C pseudogene (329/1449 nt) yaiT → pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913393,2980TC100.0% 299.2 / NA 76pseudogene (329/1449 nt)yaiTpseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/41);  total (35/41)

CGCTGATCTGGTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGTAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTCTGCTGCTAACAA  >  NC_000913/393155‑393442
                                                                                                                                               |                                                                                                                                                
cGCTGATCTGGTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAAt                                                                                                                                           <  2:765862/151‑1 (MQ=255)
cGCTGATCTGGTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAAt                                                                                                                                           <  1:138193/151‑1 (MQ=255)
cGCTGATCTGGTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAAt                                                                                                                                           >  2:261172/1‑151 (MQ=255)
   tGATCTGGTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAAtgtt                                                                                                                                        >  2:723617/1‑151 (MQ=255)
   tGATCTGGTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAAtgtt                                                                                                                                        >  2:698123/1‑151 (MQ=255)
     aTCTGGTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAAtgttgt                                                                                                                                      >  2:33536/1‑151 (MQ=255)
      tCTGGTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATgttgtt                                                                                                                                     >  2:894018/1‑151 (MQ=255)
       cTGGTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGttgttg                                                                                                                                    >  2:365466/1‑151 (MQ=255)
         ggTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAAtgttgt                                                                                                                                      <  1:358911/147‑1 (MQ=255)
         ggTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATgttgtt                                                                                                                                     <  1:984372/148‑1 (MQ=255)
         ggTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTtgc                                                                                                                                   <  2:682041/150‑1 (MQ=255)
             aaCGACCGTACCTGGGATCCTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTGCTGAACTGGATAACAGCTCTGCAGACAATGTTGTtgctg                                                                                                                                 >  2:830223/1‑148 (MQ=255)
             aaCGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCaa                                                                                                                              <  2:526753/151‑1 (MQ=255)
              aCGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAAc                                                                                                                             >  2:979870/1‑151 (MQ=255)
              aCGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAAc                                                                                                                             >  2:309248/1‑151 (MQ=255)
                  ccGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGt                                                                                                                         <  1:736349/151‑1 (MQ=255)
                   cGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTa                                                                                                                        <  1:1050322/151‑1 (MQ=255)
                    gTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTaa                                                                                                                       <  2:912063/151‑1 (MQ=255)
                          gggATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACa                                                                                                                  >  1:133128/1‑150 (MQ=255)
                            gATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACaa                                                                                                                 >  2:558163/1‑149 (MQ=255)
                             aTACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTTCaaa                                                                                                                <  2:414156/149‑1 (MQ=255)
                               aCTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTc                                                                                                            <  1:75638/151‑1 (MQ=255)
                                  tCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGt                                                                                                          <  1:482087/150‑1 (MQ=255)
                                  tCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGt                                                                                                          >  1:209569/1‑150 (MQ=255)
                                   ccAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTAt                                                                                                        >  2:775465/1‑151 (MQ=255)
                                     aaGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAGCCGGTAACTACAAAGTTCGTAt                                                                                                        <  2:496004/149‑1 (MQ=255)
                                     aaGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCg                                                                                                      <  1:779084/151‑1 (MQ=255)
                                          tAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAAc                                                                                                 >  1:908101/1‑151 (MQ=255)
                                             cTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGc                                                                                               >  1:837861/1‑150 (MQ=255)
                                               acggttacggtATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCa                                                                                              <  2:675754/149‑1 (MQ=255)
                                                  gttacggtATTGTTGCGATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAActgg                                                                                         <  2:1020812/151‑1 (MQ=255)
                                                     acggtATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAActggcgc                                                                                      >  2:1020019/1‑151 (MQ=255)
                                                       ggtATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAActggcgctg                                                                                    >  1:1060099/1‑151 (MQ=255)
                                                       ggtATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAActggcgctg                                                                                    <  2:800356/151‑1 (MQ=255)
                                                        gtATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAActggcgctgg                                                                                   <  2:468071/151‑1 (MQ=255)
                                                                gCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATc                                                                            <  1:655010/150‑1 (MQ=255)
                                                                 cTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGc                                                                          >  1:362054/1‑151 (MQ=255)
                                                                   aTGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTg                                                                        <  2:2037/151‑1 (MQ=255)
                                                                   aTGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTg                                                                        >  2:206876/1‑151 (MQ=255)
                                                                    tGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTg                                                                        >  1:556453/1‑150 (MQ=255)
                                                                       aCTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGAtt                                                                     >  2:630618/1‑150 (MQ=255)
                                                                       aCTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGAtt                                                                     >  2:551085/1‑150 (MQ=255)
                                                                       aCTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGAtt                                                                     >  1:506133/1‑150 (MQ=255)
                                                                          ctGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACa                                                                  >  2:568567/1‑150 (MQ=255)
                                                                            gATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACa                                                                  <  1:752180/148‑1 (MQ=255)
                                                                             aTGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGa                                                              <  1:554576/151‑1 (MQ=255)
                                                                               ggTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATa                                                            <  2:79857/151‑1 (MQ=255)
                                                                                  cACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATaaa                                                          <  2:133123/150‑1 (MQ=255)
                                                                                  cACCTGACTATCAACGGTAACGGCAACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATaa                                                           >  2:69815/1‑149 (MQ=255)
                                                                                    ccTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGa                                                        >  1:68861/1‑150 (MQ=255)
                                                                                       gACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAAt                                                     <  2:87466/150‑1 (MQ=255)
                                                                                       gACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAAt                                                     <  1:665588/150‑1 (MQ=255)
                                                                                            tCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCt                                                >  2:840184/1‑150 (MQ=255)
                                                                                               aCGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGt                                            <  1:739262/151‑1 (MQ=255)
                                                                                               aCGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGt                                            <  2:452379/151‑1 (MQ=255)
                                                                                                  gTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTaaa                                         <  2:26252/151‑1 (MQ=255)
                                                                                                        gCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGt                                   <  2:281124/151‑1 (MQ=255)
                                                                                                         cGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTc                                  <  1:730736/151‑1 (MQ=255)
                                                                                                                gACAACGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCa                           <  2:233289/151‑1 (MQ=255)
                                                                                                                   aaCGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAAcg                        >  1:97359/1‑151 (MQ=255)
                                                                                                                    aCGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAAcgc                       <  1:880858/151‑1 (MQ=255)
                                                                                                                    aCGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAAcgc                       <  1:125391/151‑1 (MQ=255)
                                                                                                                     cGGTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAAcgcg                      >  1:460817/1‑151 (MQ=255)
                                                                                                                      ggTACTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGa                     <  2:155410/151‑1 (MQ=255)
                                                                                                                         aCTGAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCt                  >  2:98093/1‑151 (MQ=255)
                                                                                                                            gAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCtt                 <  2:330078/149‑1 (MQ=255)
                                                                                                                            gAACTGGATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTtc                <  2:1030527/150‑1 (MQ=255)
                                                                                                                                 ggATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTctgct          <  2:698714/151‑1 (MQ=255)
                                                                                                                                 ggATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTctgct          >  1:910546/1‑151 (MQ=255)
                                                                                                                                 ggATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTctg            <  1:943043/149‑1 (MQ=255)
                                                                                                                                  gATAACAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTctgctg         >  1:163846/1‑151 (MQ=255)
                                                                                                                                      aCAGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTCTGCTGCTa      >  2:115374/1‑150 (MQ=255)
                                                                                                                                        aGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTCTGCTGCTAACa   <  1:915073/151‑1 (MQ=255)
                                                                                                                                        aGCTCTGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTCTGCTGCTAACa   <  2:846367/151‑1 (MQ=255)
                                                                                                                                           tctGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTCTGCTGCTAACaa  >  1:91428/1‑149 (MQ=255)
                                                                                                                                           tctGCAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTCTGCTGCTAACaa  >  1:134663/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CGCTGATCTGGTGAACGACCGTACCTGGGATACTTCCAAGTCTAACTACGGTTACGGTATTGTTGCTATGAACTCTGATGGTCACCTGACTATCAACGGTAACGGCGACGTAGACAACGGTACTGAACTGGATAACAGCTCTGTAGACAATGTTGTTGCTGCAACCGGTAACTACAAAGTTCGTATCGACAACGCAACTGGCGCTGGCGCTATCGCTGATTACAAAGATAAAGAAATTATCTACGTAAACGACGTCAACAGCAACGCGACCTTCTCTGCTGCTAACAA  >  NC_000913/393155‑393442

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