Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,965,089 A→G 5.4% V368A (GTG→GCG) ‡ flhB ← flagellin export apparatus, substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,965,0890AG 5.4% 368.3 / 6.6 110V368A (GTG→GCG) ‡flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (50/54);  new base G (3/3);  total (53/57)
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.99e-01

Only 100 of 112 total aligned reads displayed.

AGCACCATCATCGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACA  >  NC_000913/1964941‑1965230
                                                                                                                                                    |                                                                                                                                             
agcaCCATCATCGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCa                                                                                                                                               >  2:237205/1‑149 (MQ=255)
     catcatCGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGa                                                                                                                                          >  2:928575/1‑149 (MQ=255)
      atcatcGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGaa                                                                                                                                         <  3:60257/149‑1 (MQ=255)
      atcatcGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGaa                                                                                                                                         <  3:383834/149‑1 (MQ=255)
        catcGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGa                                                                                                                                       <  2:20756/149‑1 (MQ=255)
        catcGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGa                                                                                                                                       <  3:826017/149‑1 (MQ=255)
         atcGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGAt                                                                                                                                      <  2:37306/149‑1 (MQ=255)
         atcGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGAt                                                                                                                                      >  1:614221/1‑149 (MQ=255)
         atcGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGAt                                                                                                                                      <  2:247249/149‑1 (MQ=255)
          tcGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATg                                                                                                                                     <  1:507015/149‑1 (MQ=255)
          tcGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATg                                                                                                                                     >  3:919937/1‑149 (MQ=255)
           cGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGTGCCGGAAgccct                                                                                                                                    <  2:372211/149‑5 (MQ=255)
           cGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGa                                                                                                                                    <  1:480798/149‑1 (MQ=255)
            gACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGTGCCGGAAgccctg                                                                                                                                   <  2:846685/149‑6 (MQ=255)
             aCAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGt                                                                                                                                  >  1:82823/1‑149 (MQ=255)
              cAAGATCAACATGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGTCACCGGAAGATGAGTa                                                                                                                                 <  2:831623/149‑1 (MQ=255)
              cAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTa                                                                                                                                 <  3:565379/149‑1 (MQ=255)
                  aTCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGtt                                                                                                                             <  3:298642/149‑1 (MQ=255)
                  aTCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGtt                                                                                                                             >  4:569968/1‑149 (MQ=255)
                  aTCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGtt                                                                                                                             <  1:639650/149‑1 (MQ=255)
                  aTCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGtt                                                                                                                             >  3:502240/1‑149 (MQ=255)
                     aaCAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTa                                                                                                                          <  1:525486/149‑1 (MQ=255)
                       cAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACa                                                                                                                        >  1:250351/1‑149 (MQ=255)
                        aGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAg                                                                                                                       >  1:82407/1‑149 (MQ=255)
                          gATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAggg                                                                                                                     <  3:611222/149‑1 (MQ=255)
                           aTCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGgc                                                                                                                    >  2:498935/1‑149 (MQ=255)
                           aTCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGgc                                                                                                                    >  4:556766/1‑149 (MQ=255)
                            tCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTCGGTTGTACAGGgcg                                                                                                                   <  3:69744/149‑1 (MQ=255)
                            tCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGgcg                                                                                                                   >  4:865976/1‑149 (MQ=255)
                             cAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGgcgc                                                                                                                  >  1:779164/1‑149 (MQ=255)
                              aaaaTCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCt                                                                                                                 >  4:641524/1‑149 (MQ=255)
                                    ggTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTccac                                                                                                           <  3:319155/149‑1 (MQ=255)
                                         ggCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCATCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCAccagcc                                                                                                      <  4:536517/149‑1 (MQ=255)
                                            aaGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGa                                                                                                   >  4:309795/1‑149 (MQ=255)
                                             aGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGAc                                                                                                  <  3:262264/149‑1 (MQ=255)
                                                aTCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGcc                                                                                               <  3:92882/149‑1 (MQ=255)
                                                 tCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCa                                                                                              <  1:768255/149‑1 (MQ=255)
                                                    gCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGAGTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCAGAAGATGAGTAGATTGTACAGGGCGCTGTCCACTAGCCAGACTCCAGCg                                                                                           >  1:725497/1‑149 (MQ=255)
                                                     ccATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGt                                                                                          >  3:832631/1‑149 (MQ=255)
                                                      cATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGtt                                                                                         >  4:16179/1‑149 (MQ=255)
                                                       aTTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGttt                                                                                        <  1:237205/149‑1 (MQ=255)
                                                        ttGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTc                                                                                       <  3:556766/149‑1 (MQ=255)
                                                          gtgtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAg                                                                                     >  2:629860/1‑149 (MQ=255)
                                                           tgtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAgg                                                                                                                               <  4:834702/106‑1 (MQ=255)
                                                           tgtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAgg                                                                                                                               >  3:834702/1‑106 (MQ=255)
                                                            gtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGtt                                                                                   <  2:779164/149‑1 (MQ=255)
                                                            gtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGtt                                                                                   <  4:129511/149‑1 (MQ=255)
                                                            gtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGtt                                                                                   <  3:51173/149‑1 (MQ=255)
                                                            gtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGtt                                                                                   <  2:614221/149‑1 (MQ=255)
                                                              cGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGc                                                                                 <  1:653622/149‑1 (MQ=255)
                                                               gATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGcc                                                                                >  3:93542/1‑149 (MQ=255)
                                                               gATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGcc                                                                                <  3:743962/149‑1 (MQ=255)
                                                                aTTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCa                                                                               >  4:670783/1‑149 (MQ=255)
                                                                       gTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAg                                                                        >  4:499377/1‑149 (MQ=255)
                                                                       gTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAg                                                                        <  1:891703/149‑1 (MQ=255)
                                                                        tttGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAgg                                                                       <  3:809339/149‑1 (MQ=255)
                                                                             gggCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGc                                                                  >  1:265097/1‑149 (MQ=255)
                                                                                   gcgcAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGAc                                                                                                  <  4:459760/111‑1 (MQ=255)
                                                                                   gcgcAGCATCGCGGCCAGATTACTCATGGGTCGGCTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTcc                                                            >  2:609145/1‑149 (MQ=255)
                                                                                     gcagcaTCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGc                                                          >  3:207221/1‑149 (MQ=255)
                                                                                               ggCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCgcggcg                                                <  1:629860/149‑1 (MQ=255)
                                                                                                gCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGt                                               >  1:371970/1‑149 (MQ=255)
                                                                                                 ccAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGAc                                                                                                  >  4:711431/1‑97 (MQ=255)
                                                                                                 ccAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGAc                                                                                                  <  3:711431/97‑1 (MQ=255)
                                                                                                 ccAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGcc                                                                                <  4:674966/115‑1 (MQ=255)
                                                                                                  cAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTAc                                             >  4:299665/1‑149 (MQ=255)
                                                                                                   aGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACa                                            <  2:82823/149‑1 (MQ=255)
                                                                                                    gATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTccac                                                                                                           >  1:49347/1‑85 (MQ=255)
                                                                                                    gATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTccac                                                                                                           <  2:49347/85‑1 (MQ=255)
                                                                                                          tCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGAcc                                     >  4:709491/1‑149 (MQ=255)
                                                                                                          tCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGAcc                                     >  1:175636/1‑149 (MQ=255)
                                                                                                          tCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGAcc                                     >  2:185642/1‑149 (MQ=255)
                                                                                                                 tCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGAtt                              >  3:411151/1‑149 (MQ=255)
                                                                                                                 tCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGAtt                              >  4:767094/1‑149 (MQ=255)
                                                                                                                     tttCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttgtt                          <  3:605140/149‑1 (MQ=255)
                                                                                                                     tttCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttgtt                          <  2:52192/149‑1 (MQ=255)
                                                                                                                     tttCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttgtt                          <  4:832631/149‑1 (MQ=255)
                                                                                                                     tttCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATtcgtt                          <  2:82407/149‑1 (MQ=255)
                                                                                                                       tctcGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGa                        <  4:207221/149‑1 (MQ=255)
                                                                                                                       tctcGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGa                        <  3:865976/149‑1 (MQ=255)
                                                                                                                        ctcGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGAc                       <  1:170705/149‑1 (MQ=255)
                                                                                                                         tcGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGAcc                      <  4:90630/149‑1 (MQ=255)
                                                                                                                           gTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCaa                    >  3:257456/1‑149 (MQ=255)
                                                                                                                           gTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCaa                    >  2:104151/1‑149 (MQ=255)
                                                                                                                           gTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCAAGACCCAGGCCAGCACTTCCAACACCGCGGCGTACAGTTGACCCGAGATTTCTTGACCaa                    >  2:126597/1‑149 (MQ=255)
                                                                                                                             taataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAAtc                  <  2:250351/149‑1 (MQ=255)
                                                                                                                               ataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAAtctc                <  1:899503/149‑1 (MQ=255)
                                                                                                                               ataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAAtctc                <  2:265097/149‑1 (MQ=255)
                                                                                                                                taaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGTGATTTGTTGACCAATCTcc               >  4:913002/1‑149 (MQ=255)
                                                                                                                                taaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTcc               >  3:427346/1‑149 (MQ=255)
                                                                                                                                taaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTCTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTcc               >  4:668524/1‑149 (MQ=255)
                                                                                                                                 aaaaTCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCg              >  2:7525/1‑149 (MQ=255)
                                                                                                                                       cAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTc        <  2:554593/149‑1 (MQ=255)
                                                                                                                                        tcatcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttg                            >  2:580166/5‑128 (MQ=255)
                                                                                                                                        tcatcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttg                            <  1:580166/124‑1 (MQ=255)
                                                                                                                                        aGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCg       >  2:824167/1‑149 (MQ=255)
                                                                                                                                        aGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCg       >  3:47485/1‑149 (MQ=255)
                                                                                                                                           tcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATa    >  1:541883/2‑149 (MQ=255)
                                                                                                                                           tcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATa    >  3:379021/2‑149 (MQ=255)
                                                                                                                                             ttCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACa  >  1:271280/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
AGCACCATCATCGACAAGATCAACAGGATCAAAATCGGTCCGGCAAGGATCTGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACA  >  NC_000913/1964941‑1965230

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