Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,977,227 T→A 7.1% intergenic (‑88/+39) motA ← / ← flhC proton conductor component of flagella motor/flagellar class II regulon transcriptional activator, with FlhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,977,2270TA 7.1% 419.9 / 17.4 126intergenic (‑88/+39)motA/flhCproton conductor component of flagella motor/flagellar class II regulon transcriptional activator, with FlhD
Reads supporting (aligned to +/- strand):  ref base T (58/59);  new base A (5/4);  total (63/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

Only 100 of 137 total aligned reads displayed.

TCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGC  >  NC_000913/1977079‑1977375
                                                                                                                                                     |                                                                                                                                                    
tcatcaAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgt                                                                                                                                                      >  4:273671/1‑149 (MQ=255)
 catcaAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgtc                                                                                                                                                     <  4:699289/149‑1 (MQ=255)
 catcaAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgtc                                                                                                                                                     <  4:642569/149‑1 (MQ=255)
 catcaAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgtc                                                                                                                                                     <  1:378419/149‑1 (MQ=255)
  atcaAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgtcg                                                                                                                                                    >  3:382347/1‑149 (MQ=255)
  atcaAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgtcg                                                                                                                                                    >  4:1042638/1‑149 (MQ=255)
   tcaAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgtcgt                                                                                                                                                   >  1:425876/1‑149 (MQ=255)
   tcaAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgtcgt                                                                                                                                                   >  2:215799/1‑149 (MQ=255)
     aaaTAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTa                                                                                                                                                 <  4:382347/149‑1 (MQ=255)
     aaaTAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTa                                                                                                                                                 >  2:425175/1‑149 (MQ=255)
           ccgccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgctg                                                                                                                                           >  2:347040/1‑149 (MQ=255)
            cgccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgctgc                                                                                                                                          <  3:382077/149‑1 (MQ=255)
             gccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgctgct                                                                                                                                         >  3:59131/1‑149 (MQ=255)
             gccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgctgct                                                                                                                                         <  3:273671/149‑1 (MQ=255)
                 aaaaCTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgcttttt                                                                                                                                     <  1:280403/149‑5 (MQ=255)
                       gTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGATGCTGGAATGTTgc                                                                                                                               <  4:572645/149‑1 (MQ=255)
                       gTACCGAGAACAACCAGGTAACCTAAGAAGATAAGCACGACATCATCCTTCCACGGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTgc                                                                                                                               <  2:374262/149‑1 (MQ=255)
                        tACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGtt                                                                                                                                                  >  3:338193/1‑129 (MQ=255)
                        tACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTgcg                                                                                                                              >  1:32923/1‑149 (MQ=255)
                        tACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTgcg                                                                                                                              >  3:783575/1‑149 (MQ=255)
                          ccGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgc                                                                                                                                          <  2:614425/135‑1 (MQ=255)
                          ccGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgc                                                                                                                                          >  1:614425/1‑135 (MQ=255)
                               aacaacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAcc                                                                                                                       >  1:136050/1‑149 (MQ=255)
                                 caacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                                                                                                     >  3:773514/1‑149 (MQ=255)
                                  aacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgc                                                                                                                                          <  4:256511/127‑1 (MQ=255)
                                  aacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgc                                                                                                                                          >  3:256511/1‑127 (MQ=255)
                                  aacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGaa                                                                                                                                     >  1:330524/1‑132 (MQ=255)
                                  aacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGaa                                                                                                                                     <  2:330524/132‑1 (MQ=255)
                                   acCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTAt                                                                                                                   <  1:811603/149‑1 (MQ=255)
                                     cAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa                                                                                                                 <  1:425175/149‑1 (MQ=255)
                                       ggTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGt                                                                                                               <  2:764897/149‑1 (MQ=255)
                                       ggTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGt                                                                                                               >  2:817062/1‑149 (MQ=255)
                                           aCCTAATAAGATAAGAACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaaa                                                                                                           <  1:448361/149‑1 (MQ=255)
                                             cTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCTTCACCGTATCAGTTAAACa                                                                                                         >  2:193121/1‑149 (MQ=255)
                                              taataaGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAg                                                                                                        >  3:104105/1‑149 (MQ=255)
                                               aacaagataagCACGACATCATCCTTCAACTATTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCAATTGGGGCAAAAAAAAGCAGAGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGc                                                                                                       <  4:709096/146‑1 (MQ=255)
                                                ataagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                                                                      <  4:378267/149‑1 (MQ=255)
                                                 taagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                                                                                     <  3:839424/149‑1 (MQ=255)
                                                     ataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAc                                                                                                 <  2:425876/149‑1 (MQ=255)
                                                     ataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAc                                                                                                 >  3:158129/1‑149 (MQ=255)
                                                       aagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActc                                                                                               <  1:616286/149‑1 (MQ=255)
                                                         gCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctc                                                                                             >  2:401917/1‑149 (MQ=255)
                                                           aCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTg                                                                                           <  2:855894/149‑1 (MQ=255)
                                                                atcatcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCAt                                                                                      >  2:213229/1‑149 (MQ=255)
                                                                 tcatcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATc                                                                                     >  1:703911/1‑149 (MQ=255)
                                                                  catcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATcc                                                                                    <  2:586671/149‑1 (MQ=255)
                                                                   atcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCca                                                                                   >  1:984350/1‑149 (MQ=255)
                                                                   atcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCca                                                                                   >  1:624281/1‑149 (MQ=255)
                                                                     ccTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagc                                                                                 >  3:803373/1‑149 (MQ=255)
                                                                       ttCCACTGTTGGCCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagcag                                                                               >  1:627333/1‑149 (MQ=255)
                                                                        tCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGt                                                                              >  4:822949/1‑149 (MQ=255)
                                                                         ccACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGtt                                                                             <  1:48767/149‑1 (MQ=255)
                                                                          cACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtg                                                                            >  1:829771/1‑149 (MQ=255)
                                                                           aCTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtgt                                                                           <  4:15026/149‑1 (MQ=255)
                                                                            cTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtgtg                                                                          <  3:1024506/149‑1 (MQ=255)
                                                                             tGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTgg                                                                         <  2:260298/149‑1 (MQ=255)
                                                                             tGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTgg                                                                         >  3:464312/1‑149 (MQ=255)
                                                                             tGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTgg                                                                         >  4:851647/1‑149 (MQ=255)
                                                                                     aTGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGC‑AAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc                                                                >  4:641979/1‑149 (MQ=255)
                                                                                      tGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc                                                                <  1:364909/149‑1 (MQ=255)
                                                                                           ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAg                                                           >  1:295121/1‑149 (MQ=255)
                                                                                           ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAg                                                           >  4:563459/1‑149 (MQ=255)
                                                                                           ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAg                                                           >  4:833357/1‑149 (MQ=255)
                                                                                                ttCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAACAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGGCTGTACTCTCTGTTCATACAGACGTTGTGGGATAATATCGGCAGGATTc                                                      >  1:378903/1‑149 (MQ=255)
                                                                                                 tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCt                                                     <  1:101446/149‑1 (MQ=255)
                                                                                                     tcgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGa                                                 <  1:130873/149‑1 (MQ=255)
                                                                                                     tcgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGa                                                 >  1:669399/1‑149 (MQ=255)
                                                                                                      cgtcAGGCGTTAACGCGCGATTGGGGC‑AAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGaaa                                               <  3:641979/149‑1 (MQ=255)
                                                                                                        tcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAg                                              >  2:782699/1‑149 (MQ=255)
                                                                                                         cAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGt                                             <  3:165061/149‑1 (MQ=255)
                                                                                                          aGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGtt                                            >  2:843115/1‑149 (MQ=255)
                                                                                                           ggCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGttt                                           >  1:901297/1‑149 (MQ=255)
                                                                                                                   cgcgcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCttt                                   <  4:146211/149‑1 (MQ=255)
                                                                                                                    gcgcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                                                                      >  1:573669/1‑81 (MQ=255)
                                                                                                                     cgcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCttttt                                 <  4:803373/149‑1 (MQ=255)
                                                                                                                        gATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACt                              >  1:817028/1‑149 (MQ=255)
                                                                                                                               gCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCggg                       >  4:903062/1‑149 (MQ=255)
                                                                                                                                cAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctct                                                                                            <  1:970908/79‑1 (MQ=255)
                                                                                                                                cAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctct                                                                                            >  2:970908/1‑79 (MQ=255)
                                                                                                                                cAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGa                      >  1:969334/1‑149 (MQ=255)
                                                                                                                                cAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGa                      <  2:945592/149‑1 (MQ=255)
                                                                                                                                cAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGa                      >  2:384702/1‑149 (MQ=255)
                                                                                                                                 aaaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGa                      <  3:583772/148‑1 (MQ=255)
                                                                                                                                 aaaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGa                      >  4:583772/1‑148 (MQ=255)
                                                                                                                                 aaaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGAt                     >  1:424458/1‑149 (MQ=255)
                                                                                                                                 aaaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAAAAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGAt                     >  2:248319/1‑149 (MQ=255)
                                                                                                                                  aaaaaaaGCAGTGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATg                    >  2:726979/1‑149 (MQ=255)
                                                                                                                                  aaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATg                    <  4:732432/149‑1 (MQ=255)
                                                                                                                                     aaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCg                 >  1:898394/1‑149 (MQ=255)
                                                                                                                                       cagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGttt                                           >  2:719062/2‑121 (MQ=255)
                                                                                                                                       cagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGttt                                           <  1:719062/120‑1 (MQ=255)
                                                                                                                                       aaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGt               <  2:32923/149‑1 (MQ=255)
                                                                                                                                        agcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGtt              >  4:475863/1‑149 (MQ=255)
                                                                                                                                        agcagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGtt              <  1:705047/149‑1 (MQ=255)
                                                                                                                                        agcagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGtt              <  2:969334/149‑1 (MQ=255)
                                                                                                                                         gcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTg             >  2:992183/1‑149 (MQ=255)
                                                                                                                                                t‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAg      <  2:607844/149‑1 (MQ=255)
                                                                                                                                                  aCGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGc     <  4:102053/149‑1 (MQ=255)
                                                                                                                                                    gtcgtTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTg   >  2:616621/1‑149 (MQ=255)
                                                                                                                                                     tcgtTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGc  >  1:226510/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGC  >  NC_000913/1977079‑1977375

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