Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,003,932 A→G S21G (AGC→GGC)  fliD → flagellar filament capping protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,003,9320AG99.1% 394.9 / NA 105S21G (AGC→GGC) fliDflagellar filament capping protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (47/57);  total (47/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

AAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTT  >  NC_000913/2003784‑2004078
                                                                                                                                                    |                                                                                                                                                  
aaaaaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGc                                                                                                                                                  <  4:132400/151‑1 (MQ=255)
aaaaaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGc                                                                                                                                                  >  2:925383/1‑151 (MQ=255)
 aaaaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGcc                                                                                                                                                 >  2:1492369/1‑151 (MQ=255)
  aaaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCt                                                                                                                                                <  4:164231/151‑1 (MQ=255)
  aaaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCt                                                                                                                                                >  4:1462259/1‑151 (MQ=255)
   aaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTc                                                                                                                                               >  2:505165/1‑151 (MQ=255)
     aaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCAc                                                                                                                                             >  4:645273/1‑151 (MQ=255)
       aCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCAccg                                                                                                                                           >  4:910507/1‑151 (MQ=255)
       aCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCAccg                                                                                                                                           >  4:1232450/1‑151 (MQ=255)
        cAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCAccgc                                                                                                                                          >  2:841127/1‑151 (MQ=255)
         aaTTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCAccgcc                                                                                                                                         <  1:1439/151‑1 (MQ=255)
          aTTAAACGTAAACTTTGCGCAATTCAGTCCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCAccgccg                                                                                                                                        <  4:1198757/151‑1 (MQ=255)
            tAAACGTAAACTTTGCGAAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACcgccgc                                                                                                                                       <  2:1389121/150‑1 (MQ=255)
               aCGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCaa                                                                                                                                   <  2:28210/151‑1 (MQ=255)
               aCGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCa                                                                                                                                    >  3:21085/1‑150 (MQ=255)
               aCGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCa                                                                                                                                    >  4:280153/1‑150 (MQ=255)
                 gTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCaaa                                                                                                                                  >  3:343963/1‑150 (MQ=255)
                 gTAAACTTCGCGCAATTCAGAACGATATCTGCGGTATTCGTTTTACGTGGCGAAAGATAAAAGGAAATCGTATGGCAAGTATTTCTTCGCTGGGCTTCGGGTCAGGTCTGGATTTAAGCGCCATCCTTCATGGCCCCACCGCCGCGCaaaa                                                                                                                                 <  1:41855/151‑1 (MQ=255)
                   aaaCTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCaaaaa                                                                                                                                <  4:551413/150‑1 (MQ=255)
                    aaCTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGc                                                                                                                              <  2:694941/151‑1 (MQ=255)
                     aCTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCg                                                                                                                             <  2:1600705/151‑1 (MQ=255)
                       tttGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGAc                                                                                                                           >  2:438585/1‑151 (MQ=255)
                        ttGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACg                                                                                                                          >  4:196562/1‑151 (MQ=255)
                         tGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACg                                                                                                                          <  2:576261/150‑1 (MQ=255)
                          gcgcAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGc                                                                                                                         <  4:282339/150‑1 (MQ=255)
                            gcAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTaa                                                                                                                      <  4:1285820/151‑1 (MQ=255)
                                  cAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAAccccc                                                                                                                 >  4:114410/1‑150 (MQ=255)
                                  cAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAAcccca                                                                                                                 <  2:626273/150‑2 (MQ=255)
                                        gATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCa                                                                                                           >  4:863325/1‑150 (MQ=255)
                                        gATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCa                                                                                                           >  4:303378/1‑150 (MQ=255)
                                         aTAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCaaa                                                                                                         <  2:864228/151‑1 (MQ=255)
                                          tAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAAt                                                                                                        >  2:74312/1‑151 (MQ=255)
                                           aaCCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAAt                                                                                                        <  4:205954/150‑1 (MQ=255)
                                           aaCCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAAt                                                                                                        <  2:1355596/150‑1 (MQ=255)
                                               ccGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCa                                                                                                   <  1:255567/151‑1 (MQ=255)
                                                 gttATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTTTTCTTTGGGAGAAGGGTCCGGCCCCCACTAAATTTCCATCCTTGCCAGCCTCTCTGCCCCGCAAAAAGCGTCGCTAAccc                                                                                                                   <  2:696888/131‑1 (MQ=255)
                                                    aTTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAAtcgt                                                                                              <  2:138485/151‑1 (MQ=255)
                                                    aTTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAAtcgt                                                                                              <  4:136933/151‑1 (MQ=255)
                                                       cGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATcgtcgt                                                                                           <  2:1396428/151‑1 (MQ=255)
                                                        gTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGtt                                                                                          >  2:330160/1‑151 (MQ=255)
                                                          tttACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGttt                                                                                         >  4:223461/1‑150 (MQ=255)
                                                             aCGTGTCGAAAGATAAAAGGAAAGCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTAcc                                                                                      <  2:794299/150‑1 (MQ=255)
                                                                  tCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTaa                                                                                 >  3:413019/1‑150 (MQ=255)
                                                                   tgAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAAc                                                                               >  1:45242/2‑151 (MQ=255)
                                                                   cGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTaaa                                                                                <  4:192746/150‑1 (MQ=255)
                                                                   cGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAAc                                                                               <  3:113355/151‑1 (MQ=255)
                                                                      aaGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACtt                                                                             >  2:1170304/1‑150 (MQ=255)
                                                                       aGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAg                                                                           <  2:1551568/151‑1 (MQ=255)
                                                                        gATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAgc                                                                          >  4:1096200/1‑151 (MQ=255)
                                                                         aTAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAgc                                                                          <  2:1462700/150‑1 (MQ=255)
                                                                           aaaaGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAgcgc                                                                        <  2:712291/150‑1 (MQ=255)
                                                                                gAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACgg                                                                  <  3:312041/151‑1 (MQ=255)
                                                                                 aaaTCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGt                                                                 <  4:1524596/151‑1 (MQ=255)
                                                                                   aTCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTAc                                                               >  4:92190/1‑151 (MQ=255)
                                                                                     cGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACg                                                              <  2:597538/150‑1 (MQ=255)
                                                                                        aTGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTg                                                           >  2:666605/1‑150 (MQ=255)
                                                                                        aTGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGa                                                          <  4:1054232/151‑1 (MQ=255)
                                                                                          ggCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGaa                                                         <  2:883755/150‑1 (MQ=255)
                                                                                          ggCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGaa                                                         <  4:1318079/150‑1 (MQ=255)
                                                                                           gCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGaaaa                                                       <  2:19845/151‑1 (MQ=255)
                                                                                             aaGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAgc                                                     >  4:1401794/1‑151 (MQ=255)
                                                                                             aaGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAgc                                                     >  3:307328/1‑151 (MQ=255)
                                                                                             aaGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAgc                                                     >  2:943705/1‑151 (MQ=255)
                                                                                              aGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAgcg                                                    <  2:398048/151‑1 (MQ=255)
                                                                                                  tttCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCt                                                >  1:56913/1‑151 (MQ=255)
                                                                                                  tttCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCt                                                >  4:124178/1‑151 (MQ=255)
                                                                                                    tCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCt                                                <  2:245761/149‑1 (MQ=255)
                                                                                                    tCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTg                                               >  2:1298254/1‑150 (MQ=255)
                                                                                                      aTCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTgacg                                            >  4:1064030/1‑151 (MQ=255)
                                                                                                       tCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTgacg                                            >  1:361843/1‑150 (MQ=255)
                                                                                                        cGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTgacga                                           <  2:685543/150‑1 (MQ=255)
                                                                                                        cGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTgacga                                           <  1:413002/150‑1 (MQ=255)
                                                                                                         gCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACt                                         >  2:1504561/1‑151 (MQ=255)
                                                                                                            tggAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTc                                      >  2:1545195/2‑151 (MQ=255)
                                                                                                            gggAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACttt                                       <  4:636460/150‑1 (MQ=255)
                                                                                                               aGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTcc                                     >  2:988800/1‑149 (MQ=255)
                                                                                                               aGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAg                                   <  4:39009/151‑1 (MQ=255)
                                                                                                                gTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTcc                                     >  2:1407566/1‑148 (MQ=255)
                                                                                                                  cGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGAcc                                >  1:426699/1‑151 (MQ=255)
                                                                                                                  cGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGAcc                                <  2:1460432/151‑1 (MQ=255)
                                                                                                                   gggTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCCCCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGAcc                                <  4:1506674/150‑1 (MQ=255)
                                                                                                                   gggTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCCCCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGAcc                                <  2:1534803/150‑1 (MQ=255)
                                                                                                                    ggTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGa                                                          >  1:199979/1‑123 (MQ=255)
                                                                                                                    ggTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGa                                                          <  3:199993/123‑1 (MQ=255)
                                                                                                                     gTCAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGcc                             >  4:1470739/1‑151 (MQ=255)
                                                                                                                       cAGGTCTGGATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCa                            <  4:1409326/150‑1 (MQ=255)
                                                                                                                              ggATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACt                       <  1:206872/148‑1 (MQ=255)
                                                                                                                              ggATTTAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGc                     <  2:644402/150‑1 (MQ=255)
                                                                                                                                aTTTAAGTCCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAACCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCAt                   <  2:581691/150‑1 (MQ=255)
                                                                                                                                 tttAAGTTCCATCCTTGATGGCCTCCCCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCAt                   <  2:625083/149‑1 (MQ=255)
                                                                                                                                 tttAAGTTCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCAtt                  <  4:818136/150‑1 (MQ=255)
                                                                                                                                       ttCCATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTaa            <  2:995500/150‑1 (MQ=255)
                                                                                                                                         ccATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAACCTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAg          <  2:400090/150‑1 (MQ=255)
                                                                                                                                          cATCCTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGcc        <  4:201402/151‑1 (MQ=255)
                                                                                                                                             ccTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGa      <  4:926574/150‑1 (MQ=255)
                                                                                                                                             ccTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGAt     >  4:459717/1‑151 (MQ=255)
                                                                                                                                              cTTGATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATc    >  4:1430183/1‑151 (MQ=255)
                                                                                                                                                 gATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAACCTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCtt  <  1:209859/150‑1 (MQ=255)
                                                                                                                                                 gATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATc    <  2:1640608/148‑1 (MQ=255)
                                                                                                                                                 gATGGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCtt  <  4:849533/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAGCAATCGTCGTTTACCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTT  >  NC_000913/2003784‑2004078

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑