Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,020,841 G→A 14.6% R252H (CGC→CAC)  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,8410GA14.6% 227.1 / 28.3 89R252H (CGC→CAC) fliMflagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base G (48/28);  new base A (5/8);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

GTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTG  >  NC_000913/2020695‑2020989
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                                    <  1:632381/149‑1 (MQ=255)
                ggCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTcaca                                                                                                                                    >  4:647068/1‑149 (MQ=255)
                 gCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACag                                                                                                                                   >  1:434064/1‑149 (MQ=255)
                           cGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttatc                                                                                                                                            <  1:1007896/130‑5 (MQ=255)
                           cGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttatc                                                                                                                                            >  2:1007896/1‑126 (MQ=255)
                             gCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACagctggagctg                                                                                                                          >  2:99560/1‑146 (MQ=255)
                             gCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACagctggagctg                                                                                                                          <  1:99560/146‑1 (MQ=255)
                                aaTTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGcc                                                                                                                    >  4:189289/1‑149 (MQ=255)
                                  tttAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTgca                                                                                                                                             >  1:745613/1‑122 (MQ=255)
                                  tttAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTgca                                                                                                                                             <  2:745613/122‑1 (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:339116/1‑149 (MQ=255)
                                   ttAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAAc                                                                                                                 >  4:906308/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:484060/1‑149 (MQ=255)
                                                       aGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACg                                                                                             >  4:601912/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                                                                                                                                            >  4:78864/1‑95 (MQ=255)
                                                          aTGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGCGCACCAGGttatc                                                                                                                                            <  3:78864/99‑5 (MQ=255)
                                                             aTCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTc                                                                                       <  4:257368/149‑1 (MQ=255)
                                                             aTCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTc                                                                                       >  4:150163/1‑149 (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:378154/1‑149 (MQ=255)
                                                                                            cgccgcTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTc                                                        <  3:429331/149‑1 (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:900868/149‑1 (MQ=255)
                                                                                                     aaaaCTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATa                                               >  2:410187/1‑149 (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              <  1:399284/146‑1 (MQ=255)
                                                                                                                                      tgcaCCTGGCGCACCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGAc              >  2:399284/4‑149 (MQ=255)
                                                                                                                                      aTAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGAc              <  2:395165/149‑1 (MQ=255)
                                                                                                                                         aCCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGc           >  3:552366/1‑149 (MQ=255)
                                                                                                                                           cTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGt         >  4:975151/1‑149 (MQ=255)
                                                                                                                                           cTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGt         >  4:213074/1‑149 (MQ=255)
                                                                                                                                           cTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGt         >  4:213073/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)
                                                                                                                                            tGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGc           <  1:857650/146‑1 (MQ=255)
                                                                                                                                            tGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGc           >  2:857650/1‑146 (MQ=255)
                                                                                                                                            tGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTc        >  1:886707/1‑149 (MQ=255)
                                                                                                                                                  gcCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTg  <  4:4595/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTG  >  NC_000913/2020695‑2020989

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: