Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,966,023 A→T 6.2% S57T (TCG→ACG)  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,966,0230AT 6.2% 232.7 / 4.3 91S57T (TCG→ACG) flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (42/42);  new base T (2/4);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

GAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTT  >  NC_000913/1965875‑1966170
                                                                                                                                                    |                                                                                                                                                   
gaAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGa                                                                                                                                                     >  2:2094/1‑149 (MQ=255)
    tGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGCTAATACTGTGATCAAAATGCAGCCCAGCGGAGCGCATGCCCCCCAAt                                                                                                                                                 >  2:1370519/1‑149 (MQ=255)
     gAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCACCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATc                                                                                                                                                <  1:324620/149‑1 (MQ=255)
       gcgcCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGa                                                                                                                                              <  1:1272788/149‑1 (MQ=255)
        cgcCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGAc                                                                                                                                             <  1:101032/149‑1 (MQ=255)
         gcCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACg                                                                                                                                            <  1:451236/149‑1 (MQ=255)
           caccaGCACCACGCCGCTAATCAGCGGAAGCAGCGCCAGCAAGGCTTCTCTGATAAGCAGAATAATCTGCCCGAGGATCAGATTAGGGTCATTGATAATACTGTGATCAAAATGCCGCCCAGCGGAGAGCATGCCCTACAATCGACggg                                                                                                                                          >  2:483575/1‑149 (MQ=255)
            accaccACTACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTTCCCGAGGATCAGATTCGGGTCATTGATAATACTCTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGc                                                                                                                                         <  2:702265/149‑1 (MQ=255)
             ccaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGcc                                                                                                                                        >  2:664317/1‑149 (MQ=255)
              caGCACCCCGCCGCTAATCAGCGGCAGCAGCGCCTGCATGCCTGCTCTGATCACCAGAATAATCTGCCCGAGGATCAGATTCGGGGCATTGATAATACTGTGATCAAACTCCCGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGcat                                                                                                                                       <  1:417543/149‑3 (MQ=255)
               aGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGcatg                                                                                                                                      <  1:758313/149‑4 (MQ=255)
                 gctgctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGCTCAAAATGAAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCg                                                                                                                                    >  1:727911/7‑149 (MQ=255)
                  accacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGa                                                                                                                                   >  2:45482/1‑149 (MQ=255)
                   ccacGCCGCCAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAACGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGac                                                                                                                                  <  1:1559290/149‑1 (MQ=255)
                     acGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGacac                                                                                                                                <  2:964563/149‑1 (MQ=255)
                     acGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGacac                                                                                                                                >  1:1223430/1‑149 (MQ=255)
                      cgccgcTAATCATCGGCAGCTGCGCCATCATTGCTTCTCTGATCTGCAGAATATTCTGCCCGCGGATCAGATTCGTGTCATTGA‑AATACTGTGATCAAAATGCTGCGCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGacacac                                                                                                                              <  2:1454592/149‑1 (MQ=255)
                       gccgcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGacacac                                                                                                                              >  2:224150/1‑149 (MQ=255)
                        ccgcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATTATAATGTGATCAAAATGCAGCCCAGCGGATAGCATGCCCGACAATCGACGGGCCAGCGACACAcc                                                                                                                             >  2:1312684/1‑149 (MQ=255)
                                 aGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCa                                                                                                                    <  2:740347/149‑1 (MQ=255)
                                  gCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAACCCAg                                                                                                                   >  1:1040891/1‑149 (MQ=255)
                                   cGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGa                                                                                                                  <  1:254587/149‑1 (MQ=255)
                                      cagcagCGCCAGCATGGCTTCTATGATCAGCAGAATAATCTGCCCGAGGATCAGATTCTGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCAGACAATCGACGGGCCAGCGACACACCGCCAAACCATATaa                                                                                                               <  1:665000/149‑1 (MQ=255)
                                       agcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAac                                                                                                              >  2:16600/1‑149 (MQ=255)
                                       agcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGAGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAac                                                                                                              >  2:1141236/1‑149 (MQ=255)
                                          agcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGAGGGGCCAGCGACACACCGCCAAACCAGATAACACt                                                                                                           >  2:541498/1‑149 (MQ=255)
                                          agcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCAGGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACt                                                                                                           <  2:1337953/149‑1 (MQ=255)
                                          agcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGCTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCACGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAaccaa                                                                                                           >  2:1464351/1‑146 (MQ=255)
                                          agcGCCAGCAGGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACt                                                                                                           <  2:477370/149‑1 (MQ=255)
                                           gcgcCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTa                                                                                                          <  1:319554/149‑1 (MQ=255)
                                             gcCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTacc                                                                                                        >  1:906902/1‑147 (MQ=255)
                                                aGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGTAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGcc                                                                                                     >  2:193428/1‑149 (MQ=255)
                                                    tGGCTTCCCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCccacc                                                                                                 >  2:995318/1‑149 (MQ=255)
                                                     ggCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACCCCGCCACACCAGATAACACTAACGCccacca                                                                                                >  2:605847/1‑149 (MQ=255)
                                                     ggCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCAGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCccccca                                                                                                >  2:1460058/1‑149 (MQ=255)
                                                     ggCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGACCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCccacca                                                                                                >  2:599040/1‑149 (MQ=255)
                                                          ctctGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCaaa                                                                                           >  1:1568968/1‑149 (MQ=255)
                                                              gATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATca                                                                                       <  2:1568968/149‑1 (MQ=255)
                                                              gATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGACCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATca                                                                                       <  2:689522/149‑1 (MQ=255)
                                                               atcatcAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATcag                                                                                      <  1:366574/149‑1 (MQ=255)
                                                               aTCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCCGCGACACACCGCCAACCCAGATAACACTACCGCCCACCAGCAAAATcag                                                                                      >  1:420338/1‑149 (MQ=255)
                                                                 aagcagAATAATCTGCCCGAGCAGCAGATTCGTGTCATTGATAAAACTGTGATAAAAATGCAGCCCAGCGGAGCGCATGCCCGTCAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTATCGCCCTCCAGCAAAATcagca                                                                                    <  2:1019807/148‑1 (MQ=255)
                                                                    cagAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACCCCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTg                                                                                 >  2:43932/1‑149 (MQ=255)
                                                                    cagAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTg                                                                                 >  2:1240893/1‑149 (MQ=255)
                                                                    cagAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTg                                                                                 >  2:1144369/1‑149 (MQ=255)
                                                                      gAATAATCTGCCCGAGGATAAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCACATAACACTAACGCCCACCAGCAAATACAGCAGTGAg                                                                               >  2:105409/1‑149 (MQ=255)
                                                                        ataatCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGt                                                                             <  2:746769/149‑1 (MQ=255)
                                                                         taatCTGACCGAGGATCAGAATCGGGTCATAGATAATACAGTGATCAAAAAGCAGCCCAGCGGAGAGCATGCCCTACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTCACGCCCACCAGCATAATCACCAGTGAGCTc                                                                            >  2:243063/1‑149 (MQ=255)
                                                                           atCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAg                                                                          >  1:445726/1‑149 (MQ=255)
                                                                             cTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGtt                                                                        <  1:1427115/149‑1 (MQ=255)
                                                                                cccGATTATCAGATTCTGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCAATCGACTGGCCAGCGACACACCTCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTTAGGTCAGTTCAc                                                                     <  1:483575/149‑1 (MQ=255)
                                                                                cccGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCAc                                                                     <  1:950652/149‑1 (MQ=255)
                                                                                 ccGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACATCCAAACAAGATAACAATAACTCCAACCAGAAAAATCAGCACTGAGGTCAGTTAACg                                                                    >  2:311976/1‑149 (MQ=255)
                                                                                  cGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACgg                                                                   <  1:1187969/149‑1 (MQ=255)
                                                                                  cGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACgg                                                                   >  2:142454/1‑149 (MQ=255)
                                                                                          cgATTCGGGCCATTGATCATACTGGGATCAAAATGCCGCCCACCGGAGAGCATGCCCGACACTCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGGCAGTTCACGGGAACGCgg                                                           <  1:122933/148‑1 (MQ=255)
                                                                                           gATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCggg                                                          >  1:640461/1‑149 (MQ=255)
                                                                                           gATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCggg                                                          >  1:1317515/1‑149 (MQ=255)
                                                                                            aTTCGGGGCATTGATAATACTGTGATCAAACTGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCACTTCACGGGAACGCGGGa                                                         <  2:746044/149‑1 (MQ=255)
                                                                                               cGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGAttt                                                      <  1:1434397/149‑1 (MQ=255)
                                                                                                 gggcATTGATAATACTGTTATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGGCAGTTCACGGGAACGCGGGATTTGc                                                    <  2:160561/146‑1 (MQ=255)
                                                                                                 gggcATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTACCGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGc                                                    <  2:1040891/146‑1 (MQ=255)
                                                                                                   tCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGccc                                                  <  1:595625/149‑1 (MQ=255)
                                                                                                     aTTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATACCACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCct                                                >  2:255624/1‑149 (MQ=255)
                                                                                                      ttGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACCCCGCCAAACCAGATAACACTAACGCCCCCCAGCAAGATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCctc                                               >  1:1530466/1‑149 (MQ=255)
                                                                                                      ttGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTCACGCCCCCCAGCAAACGCAGCAGTGAGGGCAGTTCACGGGAACGCGGGCTTTGCCCctc                                               >  2:1524164/1‑149 (MQ=255)
                                                                                                      ttGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCctc                                               <  2:809969/149‑1 (MQ=255)
                                                                                                         ataataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACCCCGCCAAACCTGATAACACTAACGCCCACCCGCAAAATCAGCAGTGAGGTCAGTTCCCGGGAACGCGGGATTTGCGCCTCTTc                                            >  1:298105/1‑149 (MQ=255)
                                                                                                          taataCTGTTATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTcc                                           <  2:1105890/149‑1 (MQ=255)
                                                                                                          taataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACTCTACCGCCCACCAGCAAAATCAGCAGTGAGGGCAGTTCACGGGAACGCGGGATTTGCCCCTCTTcc                                           >  2:540301/1‑149 (MQ=255)
                                                                                                               cTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCCAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGCTTTGCCCCTCTGCCcgcgg                                      >  2:1324315/1‑148 (MQ=255)
                                                                                                                  tgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCttt                                   <  2:34730/149‑1 (MQ=255)
                                                                                                                    aTCAAAATGCAGCCCCCCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCttttt                                 <  1:135312/149‑1 (MQ=255)
                                                                                                                    aTCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGTCCGCCGAATTTGCCCCTCTTCCCGCGCgtttt                                 >  2:857139/1‑149 (MQ=255)
                                                                                                                     tCAAATTGCAGCCCCCCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCATAATCAGCTGTGAGGTCAGTTCACGGGTACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTc                                <  2:1066219/149‑1 (MQ=255)
                                                                                                                     tCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTc                                <  1:142454/149‑1 (MQ=255)
                                                                                                                      cAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCt                               >  1:1346913/1‑149 (MQ=255)
                                                                                                                       aaaaTGCAGCCCAGCGGAGAGCATGCCCGACAATCGACTGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCGACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTa                              <  1:300406/149‑1 (MQ=255)
                                                                                                                               gCCCAGCGGAGAGCATGTCCGACAATCGTCGGGCCAGCGACACTCCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGa                      <  1:862268/149‑1 (MQ=255)
                                                                                                                                  cAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGgt                   <  1:664317/149‑1 (MQ=255)
                                                                                                                                   aGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTACCGCGCTGTTTCTAGTCGGTGAGgtg                  >  2:1529466/1‑149 (MQ=255)
                                                                                                                                      ggAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTgg               <  1:2094/149‑1 (MQ=255)
                                                                                                                                       gagaGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTggg              >  1:402326/1‑149 (MQ=255)
                                                                                                                                         tcgctggCCCGTCGATTGTCGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTggggg            >  2:1258958/7‑149 (MQ=255)
                                                                                                                                          agCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCAGCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTATTCCCGCGCTTTTTCTAGTCGGTGAGGGGTGGGGGc           <  2:1121368/149‑1 (MQ=255)
                                                                                                                                           gctggCCCGTCGATTGTCGGGCCAGCGACACACCGCCAATCCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCGTCCCGCGCGTTTTCTAGTCGGGGAGGTGTGGGGGcg          >  1:1349862/5‑148 (MQ=255)
                                                                                                                                               gCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTg      <  2:1533922/149‑1 (MQ=255)
                                                                                                                                                cccGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGt     >  2:486793/1‑149 (MQ=255)
                                                                                                                                                 ccGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGtt    <  1:45482/149‑1 (MQ=255)
                                                                                                                                                  cGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCACTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGttt   <  1:16600/149‑1 (MQ=255)
                                                                                                                                                   gACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGtttt  >  2:90237/1‑149 (MQ=255)
                                                                                                                                                   gACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGtttt  <  1:193428/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTT  >  NC_000913/1965875‑1966170

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑