Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,438,871 A→T 5.8% N308Y (AAT→TAT)  flk → putative flagella assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,438,8710AT 5.8% 415.6 / 8.5 142N308Y (AAT→TAT) flkputative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base A (63/70);  new base T (4/5);  total (67/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 144 total aligned reads displayed.

TATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAG  >  NC_000913/2438724‑2439018
                                                                                                                                                   |                                                                                                                                                   
tatGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCACCCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAat                                                                                                                                                   >  2:1854240/1‑150 (MQ=255)
 atcgTCTTACTGCGCAATCAACGCTTACTTCTGCCCAGGCCGAGGATATTCGTAATCAGCATTATCAAGGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTAGCataata                                                                                                                                                  >  2:560357/4‑150 (MQ=255)
  tGGTCTTTCTGCGCATTCATCGCTTTCTTCTGCCCAGGCCGAGGATTTTCTTAATCAGGTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAatat                                                                                                                                                 <  1:1271662/150‑1 (MQ=255)
  tGGTCTTTCTGCGCAATCATAGTATTCTTCATCCCAGGCCTAGGATATTGATAATCAGCTTTATCAACTGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTTCTTATTCCTTTTCCACCTATGATCGACACGTTGCAAAatat                                                                                                                                                 <  2:872779/150‑1 (MQ=255)
   gggcTTTCTGCGCAATCCTCGCGTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCTATGATGGACACGTTGCAAAATATg                                                                                                                                                <  2:235159/147‑1 (MQ=255)
        ttCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAAc                                                                                                                                           >  1:2120314/1‑150 (MQ=255)
          cTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCCCCGATGATGGACACGTTGCAAAATATGGCAAcgc                                                                                                                                         >  2:2186932/1‑150 (MQ=255)
           tGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCACCCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAAcgcg                                                                                                                                        >  2:2145935/1‑150 (MQ=255)
            gcgcAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGt                                                                                                                                       >  2:1376729/1‑150 (MQ=255)
             cgcAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCCATATTTTGCAACGTGTc                                                                                                                                      <  2:92010/150‑1 (MQ=255)
               cAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACGCGTTGCCATATTTTGCAACGTGTcca                                                                                                                                    <  2:818847/150‑2 (MQ=255)
               cAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACGCGTTGCCATATTTTGCAACGTGTcca                                                                                                                                    <  2:31796/150‑2 (MQ=255)
                aaTCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACGCGTTGCCATATTTTGCAACGTGTccat                                                                                                                                   <  1:1539913/150‑3 (MQ=255)
                aaTCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCcg                                                                                                                                   <  2:2175381/150‑1 (MQ=255)
                       gCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCtgtg                                                                                                                            >  1:1049524/1‑150 (MQ=255)
                        cTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTgg                                                                                                                           <  1:549059/150‑1 (MQ=255)
                        cTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGCGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTgg                                                                                                                           <  1:941303/150‑1 (MQ=255)
                          ttcttcTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGAt                                                                                                                         <  1:704945/150‑1 (MQ=255)
                            cttctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAATGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATAc                                                                                                                       <  2:1894304/150‑1 (MQ=255)
                              tctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCACCGGCGGGTTAAAGGGATGGAGCCGCGGGTTATGCAACCGCTGCTTACTCCTTTGCCACCGCGGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTg                                                                                                                     <  1:1404237/150‑1 (MQ=255)
                                ttcccAGGCAGATGATATTCTTAATCAGCTTTATCAACGGCGGGGTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGtt                                                                                                                   <  2:1559707/148‑1 (MQ=255)
                                tGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGtt                                                                                                                   <  2:2193291/150‑1 (MQ=255)
                                tGCCCAGGACGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGtt                                                                                                                   <  2:1756876/150‑1 (MQ=255)
                                 gCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTa                                                                                                                  >  2:204696/1‑150 (MQ=255)
                                 gCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGACCCGCGCTGTGGATACTGTTa                                                                                                                  >  2:1851209/1‑150 (MQ=255)
                                   ccAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGt                                                                                                                >  2:740474/1‑150 (MQ=255)
                                    cAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTc                                                                                                               >  1:183231/1‑150 (MQ=255)
                                      ggCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCTATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTcgc                                                                                                             <  2:293180/150‑1 (MQ=255)
                                      ggCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTcgc                                                                                                             <  1:1590443/150‑1 (MQ=255)
                                      ggCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTcgc                                                                                                             <  2:1919017/150‑1 (MQ=255)
                                         cGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGTTCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTGCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGTATACTGTTAGTCGCGAt                                                                                                          <  2:1115430/150‑1 (MQ=255)
                                         cGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGAt                                                                                                          <  1:119840/150‑1 (MQ=255)
                                             gATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATc                                                                                                      <  1:48468/150‑1 (MQ=255)
                                             gATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTAAGTCGCGATTATc                                                                                                      <  1:384317/150‑1 (MQ=255)
                                                 ttCTTAATCAGCTTTATCAACTGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATCGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGa                                                                                                  >  2:655824/1‑150 (MQ=255)
                                                 ttCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGa                                                                                                  >  1:486450/1‑150 (MQ=255)
                                                    ttATTCAGCTTTATCATCGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGc                                                                                               <  1:1835313/150‑1 (MQ=255)
                                                    ttAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGc                                                                                               <  1:998464/150‑1 (MQ=255)
                                                        tCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCTGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGt                                                                                           >  2:1886487/1‑150 (MQ=255)
                                                         cAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTc                                                                                          <  1:239723/150‑1 (MQ=255)
                                                           gCTTTATCAACGGCGGGGTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctg                                                                                        <  2:1233959/150‑1 (MQ=255)
                                                            cTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCTTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctcg                                                                                       >  1:407753/1‑148 (MQ=255)
                                                              ttATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctggct                                                                                     <  1:1002174/150‑1 (MQ=255)
                                                               tATCACCGGCGGGTTAAAGGGATTGATCCGCGTGGTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctggctg                                                                                    >  2:1478995/1‑150 (MQ=255)
                                                               tATCAACGGCTGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctggctg                                                                                    <  2:1375459/150‑1 (MQ=255)
                                                                 tCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGt                                                                                  <  2:12622/150‑1 (MQ=255)
                                                                 tCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGt                                                                                  >  2:69860/1‑150 (MQ=255)
                                                                  cAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCTCCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGtt                                                                                 >  1:342974/1‑150 (MQ=255)
                                                                    aCGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCg                                                                               <  1:1994107/150‑1 (MQ=255)
                                                                       gcggGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTa                                                                            <  2:820117/150‑1 (MQ=255)
                                                                         gggTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGATCTGGCTGGTTCGTTAAc                                                                          <  1:1071300/150‑1 (MQ=255)
                                                                          gtttAAAGGGAGTGATCCGCGTGTTATTCAACCGCTTCTTAATCCTTTTCCAGCGATGATGGACACGTGGCAAAATATGTCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGTTTCGTTAAcc                                                                         <  2:2039768/148‑1 (MQ=255)
                                                                          ggTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAAcc                                                                         >  2:53500/1‑150 (MQ=255)
                                                                          ggTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCATGATGCTGGTCTGGCTGGTTCGTTAAcc                                                                         >  1:1950542/1‑150 (MQ=255)
                                                                          ggTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTACCGCGCTGTGGAGACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAAcc                                                                         >  2:124971/1‑150 (MQ=255)
                                                                           gTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAAccc                                                                        >  2:145243/1‑150 (MQ=255)
                                                                             tAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCcg                                                                      <  2:223935/150‑1 (MQ=255)
                                                                              aaaGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTCTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCcga                                                                     >  1:1285048/1‑150 (MQ=255)
                                                                                    aTTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTATTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAAt                                                               <  1:654640/150‑1 (MQ=255)
                                                                                       gATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGCTGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGAc                                                            >  2:3545/1‑150 (MQ=255)
                                                                                       gATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGAc                                                            >  1:1778398/1‑150 (MQ=255)
                                                                                         gccGCGGGTTATGCAACCGCTGCTTAATCCTGTTCCACCGATGACTGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGACGCTGGTCTGGCTGGGTCGTTAACCCCGACGAAATGACAg                                                          <  2:417787/149‑1 (MQ=255)
                                                                                           cgcgTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCTACGCGTCCCGCGCTGTGGATACTGTTTGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTa                                                        <  1:1851209/150‑1 (MQ=255)
                                                                                            gcgTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCATGATGCTGGTCTGGCTGGATCGTTAACACCGACTAAATGACAGTaa                                                       >  2:787602/1‑149 (MQ=255)
                                                                                              gtgtTATGCAACCGCTGCTTAATCCTTTTACACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCg                                                     <  2:209765/150‑1 (MQ=255)
                                                                                                 ttATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAAACCCGACGAAATGACAGTATCGTGa                                                  >  1:34512/1‑150 (MQ=255)
                                                                                                  tATGCAACCGCTGCTTACTCCTTTTCCACCGATGATGGACTCGTTGCTATATATGGCAACGCGTCCCGCTCTTTGTATACTGTTAGTCGCGATTATCCTGATGCTGTTCTGGCGGGTTCGTTAACCCCGACGAAATGACAGTATCGTGAc                                                 <  2:221240/150‑1 (MQ=255)
                                                                                                  tATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGAc                                                 >  1:791562/1‑150 (MQ=255)
                                                                                                  tATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGAc                                                 <  2:2161853/150‑1 (MQ=255)
                                                                                                  tATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGAc                                                 <  2:1724509/150‑1 (MQ=255)
                                                                                                   aTGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCTTGACa                                                >  1:1201659/1‑150 (MQ=255)
                                                                                                     tcAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTTGATACTGTTAGTCGCTATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACaat                                              <  1:1883644/149‑1 (MQ=255)
                                                                                                       aaCCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACaataa                                            >  2:1432089/1‑150 (MQ=255)
                                                                                                          cGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGCCGAAATGACAGTATCGTGACAATAATAc                                         >  1:1125090/1‑150 (MQ=255)
                                                                                                          cGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATAc                                         <  2:1965653/150‑1 (MQ=255)
                                                                                                          cGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATAc                                         <  1:1860004/150‑1 (MQ=255)
                                                                                                            ctgctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCTGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATAccc                                       <  2:1919628/150‑1 (MQ=255)
                                                                                                            ctgctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATAccc                                       <  1:590684/150‑1 (MQ=255)
                                                                                                             tgctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCa                                      <  2:2120314/150‑1 (MQ=255)
                                                                                                              gctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAAAATGGCAACGCGTCCCGCGCTGTGGATACTGATAGTCGCGATTATCCTGATGCTGGTCACGCAGGATCGTAAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAg                                     >  1:127849/1‑150 (MQ=255)
                                                                                                                 tAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTATTCGCGATTATCCTTATGCTGCTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGAATCGTGACAATAATACCCAGCAc                                  >  1:668380/1‑150 (MQ=255)
                                                                                                                   aTCCTTTTCCACCGATGATGGCCACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCAccc                                >  1:1640625/1‑149 (MQ=255)
                                                                                                                     ccTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCGTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCa                              <  2:2154660/150‑1 (MQ=255)
                                                                                                                        tatCCACCGATCATGGAAACGATGCAAAATATGCCAACGCGTCCCGCGCTGTGGATACAGTTAGTCGCGATTATCCTGAAGCTGGTCTGGCTGGATCGTAAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCagag                           <  1:280806/148‑1 (MQ=255)
                                                                                                                           ccACCGATGCTGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATc                        >  1:939939/1‑150 (MQ=255)
                                                                                                                           ccACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATAGCCATCACCGCAGAGATc                        <  2:319683/150‑1 (MQ=255)
                                                                                                                            cacaGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGAGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCg                       <  1:1376729/150‑1 (MQ=255)
                                                                                                                               cGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGccc                    <  1:1129318/150‑1 (MQ=255)
                                                                                                                                gatgatGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGcccc                   >  1:1855530/1‑150 (MQ=255)
                                                                                                                                     gggACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCaa              >  1:796579/2‑150 (MQ=255)
                                                                                                                                      ggACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAg             >  1:1961340/1‑150 (MQ=255)
                                                                                                                                       gACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGa            >  1:2129948/1‑150 (MQ=255)
                                                                                                                                        acgcGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATACTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGaa           >  1:1934507/4‑150 (MQ=255)
                                                                                                                                        acacGTTGCAAAATATGGCAACGCGTCCCGCTCTGTGGCTACTGTGAGTCGCGATTATCCTGCTGCTGGTCGGGCTGGTTCGTTAGCCCCGACGAAATGCCAGAATCGTGACAATTATTACCAGCACCGAAGAGATCGCCCCGGCAAGaa           >  1:337773/1‑150 (MQ=255)
                                                                                                                                        acacGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGaa           <  2:746881/150‑1 (MQ=255)
                                                                                                                                         cacGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGaaa          >  2:1541096/1‑150 (MQ=255)
                                                                                                                                             ttGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATAcc      <  2:342974/150‑1 (MQ=255)
                                                                                                                                              tGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCg     >  2:479614/1‑150 (MQ=255)
                                                                                                                                                cAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGaa   >  1:1900909/1‑150 (MQ=255)
                                                                                                                                                 aaaaTATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAag  <  1:1740383/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAG  >  NC_000913/2438724‑2439018

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: