Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,020,838 T→C 13.5% V251A (GTG→GCG)  fliM → flagellar motor switching and energizing component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,020,8380TC13.5% 233.6 / 26.3 89V251A (GTG→GCG) fliMflagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base T (48/29);  new base C (5/7);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

ACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTC  >  NC_000913/2020690‑2020983
                                                                                                                                                    |                                                                                                                                                 
acGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGt                                                                                                                                                   >  4:393514/1‑149 (MQ=255)
acGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGt                                                                                                                                                   >  1:848859/1‑149 (MQ=255)
     gTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCAcc                                                                                                                                              <  3:932439/149‑1 (MQ=255)
       tCCATGTTGAAATTGGCAACCTGACCGGCCAAATTAATAACTTCCCGCCATTCAGCATTATCCGGCCGCCACGGGAATTTTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAg                                                                                                                                            <  3:489948/149‑1 (MQ=255)
         cATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGAGCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGt                                                                                                                                          >  4:436871/1‑149 (MQ=255)
          aggtgGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTg                                                                                                                                         <  4:896174/147‑1 (MQ=255)
           tgtgGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTgc                                                                                                                                        >  1:112109/1‑149 (MQ=255)
            gtgGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGt                                                                                                                                          <  1:228243/146‑1 (MQ=255)
            gtgGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGt                                                                                                                                          >  2:228243/1‑146 (MQ=255)
             tgGAGATTGGCAACCTGACCGGGGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttatc                                                                                                                                      <  4:617273/149‑5 (MQ=255)
             tgGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttatc                                                                                                                                      <  2:954903/149‑5 (MQ=255)
                 gaTTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCAtt                                                                                                                                  <  1:310486/149‑1 (MQ=255)
                   ttGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTca                                                                                                                                <  3:703497/149‑1 (MQ=255)
                     ggCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTcaca                                                                                                                              >  4:647068/1‑149 (MQ=255)
                     ggCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTcaca                                                                                                                              <  1:632381/149‑1 (MQ=255)
                      gCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACag                                                                                                                             >  1:434064/1‑149 (MQ=255)
                                cGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttatc                                                                                                                                      >  2:1007896/1‑126 (MQ=255)
                                cGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttatc                                                                                                                                      <  1:1007896/130‑5 (MQ=255)
                                  gCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACagctggagctg                                                                                                                    <  1:99560/146‑1 (MQ=255)
                                  gCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACagctggagctg                                                                                                                    >  2:99560/1‑146 (MQ=255)
                                     aaTTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGcc                                                                                                              >  4:189289/1‑149 (MQ=255)
                                       tttAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTgca                                                                                                                                       <  2:745613/122‑1 (MQ=255)
                                       tttAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTgca                                                                                                                                       >  1:745613/1‑122 (MQ=255)
                                       tttAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCaa                                                                                                            >  1:79851/1‑149 (MQ=255)
                                       tttAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttat                                                                                                                                       >  4:516198/1‑119 (MQ=255)
                                       tttAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttat                                                                                                                                       <  3:516198/122‑4 (MQ=255)
                                        ttAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAAc                                                                                                           >  4:906308/1‑149 (MQ=255)
                                        ttAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAAc                                                                                                           >  4:339116/1‑149 (MQ=255)
                                          aaTATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACtt                                                                                                         >  1:664383/1‑149 (MQ=255)
                                           atatCTGCCTGCCATTCAGCGTGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACttt                                                                                                        >  2:647139/1‑149 (MQ=255)
                                              tCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGcc                                                                                                     >  3:194106/1‑149 (MQ=255)
                                                 gcctgccATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGat                                                                                                  <  3:574823/149‑1 (MQ=255)
                                                  cctgccATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTTGAAAACTCGCGTAATGAAAATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGata                                                                                                 <  4:562476/149‑1 (MQ=255)
                                                        aTTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGc                                                                                           >  4:31935/1‑149 (MQ=255)
                                                           cAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATATCCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTAc                                                                                        <  4:622560/149‑1 (MQ=255)
                                                            aGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACg                                                                                       >  4:601912/1‑149 (MQ=255)
                                                            aGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACg                                                                                       >  4:484060/1‑149 (MQ=255)
                                                             gCATGATCGAGACGATACGGGAATTGTAGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGc                                                                                      <  4:7593/149‑1 (MQ=255)
                                                              cATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGcc                                                                                     >  3:714447/1‑149 (MQ=255)
                                                               aTGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCt                                                                                    >  2:211613/1‑149 (MQ=255)
                                                               aTGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttatc                                                                                                                                      <  3:78864/99‑5 (MQ=255)
                                                               aTGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttatc                                                                                                                                      >  4:78864/1‑95 (MQ=255)
                                                                  aTCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTc                                                                                 >  4:150163/1‑149 (MQ=255)
                                                                  aTCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTc                                                                                 <  4:257368/149‑1 (MQ=255)
                                                                     gAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTGGCTACGCCTGTCGCa                                                                              >  3:831227/1‑149 (MQ=255)
                                                                     gAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCa                                                                              >  2:1034369/1‑149 (MQ=255)
                                                                         cGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTg                                                                                                                                         <  4:333003/86‑1 (MQ=255)
                                                                         cGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTg                                                                                                                                         >  3:333003/1‑86 (MQ=255)
                                                                         cGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGAtt                                                                          <  1:303727/149‑1 (MQ=255)
                                                                         cGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGAtt                                                                          <  1:978347/149‑1 (MQ=255)
                                                                               gggAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTaaaa                                                                    <  3:321910/149‑1 (MQ=255)
                                                                               gggAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTaaaa                                                                    >  2:1094461/1‑149 (MQ=255)
                                                                                 gAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACt                                                                  >  3:55414/1‑149 (MQ=255)
                                                                                   aTTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGa                                                                >  3:874606/1‑149 (MQ=255)
                                                                                    ttgttgGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGaa                                                               >  2:610912/1‑149 (MQ=255)
                                                                                          gTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGTTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCgg                                                         >  3:25628/1‑149 (MQ=255)
                                                                                          gTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCgg                                                         >  3:700215/1‑149 (MQ=255)
                                                                                                 cgccgcTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTc                                                  <  3:429331/149‑1 (MQ=255)
                                                                                                 cgccgcTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTc                                                  >  3:378154/1‑149 (MQ=255)
                                                                                                       tGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCg                                            >  2:372681/1‑149 (MQ=255)
                                                                                                         gAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGAt                                          <  4:660479/149‑1 (MQ=255)
                                                                                                          aaaaCTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATa                                         <  4:874606/149‑1 (MQ=255)
                                                                                                          aaaaCTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATa                                         >  2:410187/1‑149 (MQ=255)
                                                                                                          aaaaCTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATa                                         <  2:900868/149‑1 (MQ=255)
                                                                                                              cTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGaaa                                     >  3:68694/1‑149 (MQ=255)
                                                                                                               tggcgTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGaaaa                                    >  4:741863/3‑149 (MQ=255)
                                                                                                                   gTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAAccc                                <  1:610912/149‑1 (MQ=255)
                                                                                                                   gTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAAccc                                <  2:434064/149‑1 (MQ=255)
                                                                                                                     aaTGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGa                              >  3:1046484/1‑149 (MQ=255)
                                                                                                                          aGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGca                         <  3:484060/149‑1 (MQ=255)
                                                                                                                           gATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGcat                        >  1:255659/1‑149 (MQ=255)
                                                                                                                            aTCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGCTAGAAAAACCCGATCGcatc                       >  1:225818/1‑149 (MQ=255)
                                                                                                                            aTCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGcatc                       >  1:773202/1‑149 (MQ=255)
                                                                                                                              cAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGcatcat                     <  2:457621/149‑1 (MQ=255)
                                                                                                                              cAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGcatcat                     >  3:4595/1‑149 (MQ=255)
                                                                                                                               aGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCatcatc                    <  4:831227/149‑1 (MQ=255)
                                                                                                                                      gcgcGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCAGGAGCTGGTCGCCAACTTTGCCGATATCGCGCTACTCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATg             >  4:485962/1‑149 (MQ=255)
                                                                                                                                          gATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGa         >  1:836915/1‑149 (MQ=255)
                                                                                                                                           tgcaCCTGGCGCACCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGAc        >  2:399284/4‑149 (MQ=255)
                                                                                                                                           tgcaCCTGGCGCACCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGAc        <  1:399284/146‑1 (MQ=255)
                                                                                                                                           aTAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGAc        <  2:395165/149‑1 (MQ=255)
                                                                                                                                              aCCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGc     >  3:552366/1‑149 (MQ=255)
                                                                                                                                                cTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGt   >  3:181145/1‑149 (MQ=255)
                                                                                                                                                cTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGt   >  3:259424/1‑149 (MQ=255)
                                                                                                                                                cTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGt   <  3:601912/149‑1 (MQ=255)
                                                                                                                                                cTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGt   >  4:213073/1‑149 (MQ=255)
                                                                                                                                                cTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGt   >  4:213074/1‑149 (MQ=255)
                                                                                                                                                cTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGt   >  4:975151/1‑149 (MQ=255)
                                                                                                                                                 tGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGc     <  1:857650/146‑1 (MQ=255)
                                                                                                                                                 tGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGc     >  2:857650/1‑146 (MQ=255)
                                                                                                                                                 tGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTc  >  1:886707/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTC  >  NC_000913/2020690‑2020983

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