Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,964,425 G→C 5.3% A209G (GCC→GGC)  flhA ← putative flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,964,4250GC 5.3% 386.3 / 11.3 113A209G (GCC→GGC) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base G (58/48);  new base C (2/4);  total (61/52)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01

Only 100 of 115 total aligned reads displayed.

GCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCA  >  NC_000913/1964281‑1964572
                                                                                                                                                |                                                                                                                                                   
gcCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATc                                                                                                                                                 >  3:162044/1‑149 (MQ=255)
 ccGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATcg                                                                                                                                                >  3:830807/1‑149 (MQ=255)
     aaTTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATcgccg                                                                                                                                             <  2:298259/148‑1 (MQ=255)
     aaTTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCACATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCACGGCGATGCCATcgccg                                                                                                                                             <  1:282384/148‑1 (MQ=255)
       tttGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCcgcg                                                                                                                                           <  2:97551/148‑1 (MQ=255)
            gCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACaa                                                                                                                                      <  2:71727/148‑1 (MQ=255)
             ccaccaGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACaaa                                                                                                                                     <  2:252867/148‑1 (MQ=255)
             ccaccaGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACaaa                                                                                                                                     <  1:391986/148‑1 (MQ=255)
              caccaGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAact                                                                                                                                   >  3:343743/1‑149 (MQ=255)
              caccaGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAact                                                                                                                                   >  4:597623/1‑149 (MQ=255)
               accaGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAactt                                                                                                                                  >  3:286609/1‑149 (MQ=255)
                   gACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAacttact                                                                                                                               >  2:286026/1‑148 (MQ=255)
                      cGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGcc                                                                                                                           >  4:363115/1‑149 (MQ=255)
                       gTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGcc                                                                                                                           <  1:824801/148‑1 (MQ=255)
                          gCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGt                                                                                                                       >  3:436669/1‑149 (MQ=255)
                           ccAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGcccc                                                                                                                         <  3:293543/146‑1 (MQ=255)
                           ccAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGcccc                                                                                                                         >  4:293543/1‑146 (MQ=255)
                           ccAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGt                                                                                                                       >  1:274874/1‑148 (MQ=255)
                           ccAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTc                                                                                                                      <  4:830807/149‑1 (MQ=255)
                            cAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTc                                                                                                                      <  2:654326/148‑1 (MQ=255)
                              aTGGTTAATAGAGTATAACTTTCCGCCGCGTGCCCCTTGCTAATACCATGTTGAAGCTCCCCGACCAGCAACCCGCCGACAATGGTGATTACAATGATACGTGTTCCGGCGATGTCTTCTCACCGGACACACGCACTTGCCCCATcaac                                                                                                                   >  4:638691/1‑146 (MQ=255)
                              aTGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGt                                                                                                                       >  1:494628/1‑145 (MQ=255)
                              aTGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGt                                                                                                                       <  2:494628/145‑1 (MQ=255)
                                 gtaaaTAGCGTATAACTTTCCGCCGCATGTCCCATGCTCATGCCATGTTGCAGCACGCCGAACAACAACACGCCGACAATGTTACTGCCCATGATGAGGATCCCGGCCATGGCATCGCCGCGAACAAACTTCCTTGCCCCGTCCATTGa                                                                                                                <  4:48617/146‑1 (MQ=255)
                                 gTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATTTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTg                                                                                                                 <  1:286026/148‑1 (MQ=255)
                                 gTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGa                                                                                                                <  4:752464/149‑1 (MQ=255)
                                  tCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGa                                                                                                                >  1:217879/1‑148 (MQ=255)
                                    aaTAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGcc                                                                                                             >  3:613227/1‑149 (MQ=255)
                                       aGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGc                                                                                                                                                    <  2:164671/107‑1 (MQ=255)
                                       aGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGc                                                                                                                                                    >  1:164671/1‑107 (MQ=255)
                                             taACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAAt                                                                                                     <  1:428244/148‑1 (MQ=255)
                                             taACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATc                                                                                                    >  4:874518/1‑149 (MQ=255)
                                              aaCTTTCCGCCGCGTGTCCCATGCTCTTGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACACTGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAAAAAACTTACTTGCCCCGTCCAATGAGCCGTAAAAATc                                                                                                    <  2:456815/148‑1 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGc                                                                                                 >  1:758969/1‑148 (MQ=255)
                                                  ttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCt                                                                                                >  1:945453/1‑148 (MQ=255)
                                                     cgccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCt                                                                                            <  4:447768/149‑1 (MQ=255)
                                                     cgccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCt                                                                                            <  4:436669/149‑1 (MQ=255)
                                                     cgccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCt                                                                                            <  3:597623/149‑1 (MQ=255)
                                                       ccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCGGTAAAAATCGGCTTCCTGa                                                                                          >  4:651934/1‑149 (MQ=255)
                                                       ccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGa                                                                                          <  3:382017/149‑1 (MQ=255)
                                                       ccgcGTGTCCCATGCTCATGCCATCTTGCAGCACGCCGACCAGCAACCCGCCGATAACGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCCCGAACAAACTTACTTTCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTg                                                                                           <  2:513994/148‑1 (MQ=255)
                                                           gtgtCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCa                                                                                      >  4:666435/1‑149 (MQ=255)
                                                            tgtCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCAc                                                                                     <  4:343743/149‑1 (MQ=255)
                                                               cccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACtt                                                                                   >  1:718702/1‑148 (MQ=255)
                                                                  aTGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCgg                                                                                >  1:759314/1‑148 (MQ=255)
                                                                  aTGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGa                                                                               >  3:551099/1‑149 (MQ=255)
                                                                    gCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAgc                                                                             <  3:910905/149‑1 (MQ=255)
                                                                    gCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAgc                                                                             <  3:488356/149‑1 (MQ=255)
                                                                      tCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAgcgg                                                                           >  3:608798/1‑149 (MQ=255)
                                                                       cATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATcg                                                                                                                                                >  1:52881/1‑79 (MQ=255)
                                                                       cATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAgcgg                                                                           <  2:244502/148‑1 (MQ=255)
                                                                         tGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGt                                                                        <  3:14347/149‑1 (MQ=255)
                                                                          gCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGtt                                                                       >  3:208580/1‑149 (MQ=255)
                                                                          gCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGt                                                                        >  2:780847/1‑148 (MQ=255)
                                                                            cATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGtttt                                                                     >  4:879180/1‑149 (MQ=255)
                                                                              tGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGttttt                                                                    <  1:780847/148‑1 (MQ=255)
                                                                               gTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGTCGCGACCAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGtttttt                                                                   <  2:803467/148‑1 (MQ=255)
                                                                               gTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGtttttt                                                                   >  1:131236/1‑148 (MQ=255)
                                                                                   cagcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGcc                                                               <  1:365667/148‑1 (MQ=255)
                                                                                        cGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATc                                                          <  3:141537/148‑1 (MQ=255)
                                                                                        cGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATc                                                          >  4:141537/1‑148 (MQ=255)
                                                                                           cGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTc                                                       <  1:136009/148‑1 (MQ=255)
                                                                                               cAGCAACCCGCCGACAATGTTAATGACCATAATGAGGAACCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCTCCTCATCTTTACca                                                   <  2:142833/148‑1 (MQ=255)
                                                                                                 gCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATc                                                                                                    >  4:331959/1‑97 (MQ=255)
                                                                                                 gCAACCCGCCGACAATGTTAATGACCAGGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATc                                                                                                    <  3:331959/97‑1 (MQ=255)
                                                                                                  cAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACcaatc                                                >  1:914698/1‑148 (MQ=255)
                                                                                                    aCCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACcaatcaat                                             <  3:797881/149‑1 (MQ=255)
                                                                                                     cccGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACcaatcaat                                             >  1:169050/1‑148 (MQ=255)
                                                                                                      ccgccgACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATcc                                           <  4:823455/149‑1 (MQ=255)
                                                                                                      ccgccgACAATGTTAATGACCATGATGAGGATCCCGGCGATGGAATCGCCGCGACCAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATcc                                           <  4:665751/149‑1 (MQ=255)
                                                                                                        gccgACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCtt                                                                                               >  2:701600/1‑95 (MQ=255)
                                                                                                        gccgACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCtt                                                                                               <  1:701600/95‑1 (MQ=255)
                                                                                                              aaTGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTaa                                   >  3:844392/1‑149 (MQ=255)
                                                                                                              aaTGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTaa                                   <  4:162044/149‑1 (MQ=255)
                                                                                                               aTGTTTATTACCATGAGCAGTCTGCCGGAGATGGCATGGGCGCGAACAGACTTACTTGCCCAGTCCATTGAGCCGTAAAAATAGGCTTCCTGCGTCATTTCGGAGCTGCGTTTTTCCGTCTCATCTTCACCAATCAATCCTGCGttga                                   >  2:607447/1‑146 (MQ=255)
                                                                                                               aTGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTaa                                   <  2:556657/148‑1 (MQ=255)
                                                                                                               aTGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAg                                  <  4:890156/149‑1 (MQ=255)
                                                                                                                 gTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTAACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAgg                                 >  1:655648/1‑148 (MQ=255)
                                                                                                                        aCCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTCTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTc                         >  4:243202/1‑149 (MQ=255)
                                                                                                                         ccATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTc                         >  1:258681/1‑148 (MQ=255)
                                                                                                                         ccATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCa                        >  3:799592/1‑149 (MQ=255)
                                                                                                                          cATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGGTAAGGTCGGCGTCaa                       >  4:905135/1‑149 (MQ=255)
                                                                                                                           atgatgAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCaa                       >  2:603482/1‑148 (MQ=255)
                                                                                                                            tgatgaGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAAt                      >  1:84839/1‑148 (MQ=255)
                                                                                                                            tgatgaGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAAt                      <  1:182542/148‑1 (MQ=255)
                                                                                                                                  ggATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTCTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATc               >  4:213575/1‑149 (MQ=255)
                                                                                                                                   gATCCTGGCGATGGCATCGCGGCGAACAAACTAACTTGCCCCGTCCATGGAGCCGTAAAAAACGGCTTCCTGAGTCACTTCGGAGCGGCGATGTATGGTATCATCTTCATCAATCAAGTCGGCGATAAGCGCGGTGTCAATCGACATca              >  4:700224/1‑148 (MQ=255)
                                                                                                                                    aTCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTg             >  3:67932/1‑149 (MQ=255)
                                                                                                                                    aTCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGCCACTTTGGAGCGGCGTTTTTTCGCCTCATCTTCACCCATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTg             >  3:323245/1‑149 (MQ=255)
                                                                                                                                     tcgcGGCGATGCCATCGCCGGGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTg             >  1:171518/4‑148 (MQ=255)
                                                                                                                                     tCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGc                            <  4:424009/133‑1 (MQ=255)
                                                                                                                                     tCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGc                            >  3:424009/1‑133 (MQ=255)
                                                                                                                                     tCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGc            >  3:152391/1‑149 (MQ=255)
                                                                                                                                      cccGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCGCCAATCAACCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCt           <  4:263967/149‑1 (MQ=255)
                                                                                                                                        cGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAAt                      >  1:171711/1‑136 (MQ=255)
                                                                                                                                        cGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAAt                      <  2:171711/136‑1 (MQ=255)
                                                                                                                                        cGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCtt          >  2:500782/1‑148 (MQ=255)
                                                                                                                                           cGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTAcc       >  2:460794/1‑148 (MQ=255)
                                                                                                                                               ggCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGc   >  2:473523/1‑148 (MQ=255)
                                                                                                                                               ggCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCa  >  3:842823/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCA  >  NC_000913/1964281‑1964572

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