Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,997,314 G→T 93.1% intergenic (‑227/‑39) fliC ← / → fliD flagellar filament structural protein (flagellin)/flagellar filament capping protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,997,3140GT93.1% 862.6 / 43.2 289intergenic (‑227/‑39)fliC/fliDflagellar filament structural protein (flagellin)/flagellar filament capping protein
Reads supporting (aligned to +/- strand):  ref base G (14/7);  new base T (147/120);  total (162/127)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 289 total aligned reads displayed.

TCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAG  >  CP009273/1997170‑1997460
                                                                                                                                                |                                                                                                                                                  
  tGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCg                                                                                                                                             >  2:3006763/1‑150 (MQ=255)
   gTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGt                                                                                                                                            >  1:477734/1‑150 (MQ=255)
   gTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGt                                                                                                                                            >  1:668882/1‑150 (MQ=255)
     cGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGttt                                                                                                                                          <  2:3148931/150‑1 (MQ=255)
       tAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTa                                                                                                                                        >  1:3143280/1‑150 (MQ=255)
       tAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTAGCGCAATTCAGACCGATAACCCCAGTATTCGTTTTa                                                                                                                                        >  1:3143765/1‑150 (MQ=255)
        aaaaCAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTAc                                                                                                                                       >  1:1420766/1‑150 (MQ=255)
        aaaaCAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTAc                                                                                                                                       >  1:3686879/1‑150 (MQ=255)
            caacCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACgtgt                                                                                                                                   <  2:1588269/150‑1 (MQ=255)
              acCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCg                                                                                                                                 <  2:64850/150‑1 (MQ=255)
              acCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCg                                                                                                                                 >  1:746038/1‑150 (MQ=255)
              acCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCg                                                                                                                                 <  1:417571/150‑1 (MQ=255)
              acCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCg                                                                                                                                 <  2:2163521/150‑1 (MQ=255)
               ccTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGa                                                                                                                                >  1:4049758/1‑150 (MQ=255)
                cTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGaa                                                                                                                               <  1:3030312/150‑1 (MQ=255)
                cTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGaa                                                                                                                               >  2:3172486/1‑150 (MQ=255)
                 ttcaGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGaaa                                                                                                                              <  1:3124342/147‑1 (MQ=255)
                    aGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGAt                                                                                                                           <  2:3350818/150‑1 (MQ=255)
                    aGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGAt                                                                                                                           <  2:1837587/150‑1 (MQ=255)
                      aaTTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATaa                                                                                                                         <  2:2410359/150‑1 (MQ=255)
                      aaTTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATaa                                                                                                                         >  2:2220500/1‑150 (MQ=255)
                      aaTTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATaa                                                                                                                         >  2:11439/1‑150 (MQ=255)
                         ttttttCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAg                                                                                                                      <  2:4418911/150‑1 (MQ=255)
                         ttttttCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAg                                                                                                                      >  1:2319731/1‑150 (MQ=255)
                         ttttttCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGATTTACGTGTCGAAAGATAAAAg                                                                                                                      <  2:764608/150‑1 (MQ=255)
                         ttttttCAAAAACAGCCAATTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAg                                                                                                                      >  2:208527/1‑150 (MQ=255)
                          tttttCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAgg                                                                                                                     >  2:1466705/1‑150 (MQ=255)
                              tCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAAt                                                                                                                 <  2:2550934/150‑1 (MQ=255)
                              tCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTGTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAAt                                                                                                                 <  2:3813192/150‑1 (MQ=255)
                                aaaaaCAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCg                                                                                                               >  1:1789212/1‑150 (MQ=255)
                                  aaaCAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCa                                                                                                             >  1:4221562/1‑150 (MQ=255)
                                      aGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGc                                                                                                         <  1:3172486/150‑1 (MQ=255)
                                       gCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCa                                                                                                        >  2:1225915/1‑150 (MQ=255)
                                       gCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCaa                                                                                                       >  1:3023726/1‑150 (MQ=255)
                                        ccATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCaa                                                                                                       >  1:2224180/1‑150 (MQ=255)
                                         cATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAg                                                                                                      <  1:2595992/150‑1 (MQ=255)
                                           ttttttGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTa                                                                                                    <  1:2018533/150‑1 (MQ=255)
                                              tttGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTAttt                                                                                                 <  2:1222888/150‑1 (MQ=255)
                                                  ttAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATc                                                                                             >  1:1132715/1‑150 (MQ=255)
                                                   tAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCg                                                                                            <  2:2244123/150‑1 (MQ=255)
                                                    aGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGc                                                                                           <  1:4375114/150‑1 (MQ=255)
                                                     gTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCt                                                                                          <  1:552676/150‑1 (MQ=255)
                                                            aaaTACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTc                                                                                   >  1:536575/1‑150 (MQ=255)
                                                              aTACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCgg                                                                                 >  2:883929/1‑150 (MQ=255)
                                                                  tctTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCa                                                                             <  2:942853/150‑1 (MQ=255)
                                                                    ttttCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAgg                                                                           <  2:4117998/150‑1 (MQ=255)
                                                                           tGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGAt                                                                    <  2:3600991/150‑1 (MQ=255)
                                                                           tGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGAt                                                                    >  2:1751631/1‑150 (MQ=255)
                                                                           tGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGAt                                                                    >  2:1751131/1‑150 (MQ=255)
                                                                            gCCCCTTATTCCCGCTATTAAAAAAAACAATTCAACGTAAACTTTGCGCAATTCAGACCGATACCCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGAtt                                                                   >  2:50855/1‑150 (MQ=255)
                                                                            gCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGAtt                                                                   <  2:1053291/150‑1 (MQ=255)
                                                                            gCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGAtt                                                                   <  1:324334/150‑1 (MQ=255)
                                                                                  tATTCCCTCTATTAAAAAAAACAATTAAACTTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGtt                                                             <  2:3048503/150‑1 (MQ=255)
                                                                                  tATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGtt                                                             >  1:476987/1‑150 (MQ=255)
                                                                                  tATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGtt                                                             >  1:550513/1‑150 (MQ=255)
                                                                                   aTTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTc                                                            <  2:4311072/150‑1 (MQ=255)
                                                                                     tCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCa                                                          >  2:610637/1‑150 (MQ=255)
                                                                                        cGCTATTAAAAAAAACAATTACACGTAAACTGTGAGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAACGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATcc                                                       <  2:4348228/150‑1 (MQ=255)
                                                                                        cGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATcc                                                       <  2:197828/150‑1 (MQ=255)
                                                                                         gCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATcat                                                      >  2:2157202/1‑148 (MQ=255)
                                                                                           tATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTg                                                    >  2:3458852/1‑150 (MQ=255)
                                                                                            aTTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGa                                                   >  2:2187815/1‑150 (MQ=255)
                                                                                             ttAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGAt                                                  >  2:769725/1‑150 (MQ=255)
                                                                                              tAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATa                                                 <  2:3686868/150‑1 (MQ=255)
                                                                                              tAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATa                                                 <  1:3753552/150‑1 (MQ=255)
                                                                                                 aaaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGcc                                              >  1:1447319/1‑150 (MQ=255)
                                                                                                  aaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCt                                             <  2:260771/150‑1 (MQ=255)
                                                                                                   aaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTc                                            >  1:3296800/1‑150 (MQ=255)
                                                                                                     aaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCAc                                          <  1:1780533/150‑1 (MQ=255)
                                                                                                     aaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCAc                                          <  1:4322193/150‑1 (MQ=255)
                                                                                                        aaTTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCAccgc                                       >  2:3561376/1‑150 (MQ=255)
                                                                                                        aaTTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCAccgc                                       >  2:3900912/1‑150 (MQ=255)
                                                                                                         aTTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCAccgcc                                      >  2:3592184/1‑150 (MQ=255)
                                                                                                         aTTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCAccgcc                                      >  1:4398402/1‑150 (MQ=255)
                                                                                                            aaaCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCcgcg                                   >  2:3616793/1‑150 (MQ=255)
                                                                                                            aaaCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCcgcg                                   >  1:625484/1‑150 (MQ=255)
                                                                                                            aaaCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCcgcg                                   >  1:2909511/1‑150 (MQ=255)
                                                                                                            aaaCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCcgcg                                   >  2:1346782/1‑150 (MQ=255)
                                                                                                             aaCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCgcgc                                  >  1:1328758/1‑150 (MQ=255)
                                                                                                              aCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCa                                 <  1:4495598/150‑1 (MQ=255)
                                                                                                              aCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCa                                 <  1:3385635/150‑1 (MQ=255)
                                                                                                                gTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCaaa                               >  1:4532739/1‑150 (MQ=255)
                                                                                                                 tAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCaaaa                              >  1:3108267/1‑150 (MQ=255)
                                                                                                                  aaaCTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCaaaaa                             >  1:2482878/1‑150 (MQ=255)
                                                                                                                   aaCTTTGCGCAATTCAGACCGATAAACCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATATAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAg                            >  2:4510506/1‑150 (MQ=255)
                                                                                                                       ttGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGAc                        <  1:870800/150‑1 (MQ=255)
                                                                                                                          cgcAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCt                     >  1:1848337/1‑150 (MQ=255)
                                                                                                                            cAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTaa                   >  1:2142933/1‑150 (MQ=255)
                                                                                                                            cAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTaa                   >  1:2835708/1‑150 (MQ=255)
                                                                                                                              aTTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAAcc                 >  1:3271080/1‑150 (MQ=255)
                                                                                                                                 cAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAAccccc              <  2:2852908/150‑1 (MQ=255)
                                                                                                                                   gACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCAt            >  2:397443/1‑150 (MQ=255)
                                                                                                                                     ccGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCAttt          <  2:3023726/150‑1 (MQ=255)
                                                                                                                                       gATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCa        >  2:2923382/1‑150 (MQ=255)
                                                                                                                                       gATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCa        >  1:1377850/1‑150 (MQ=255)
                                                                                                                                        aTAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCaa       >  2:3806729/1‑150 (MQ=255)
                                                                                                                                           aCCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGTAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATc    <  2:476997/150‑1 (MQ=255)
                                                                                                                                            ccccTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCa   <  2:1675995/150‑1 (MQ=255)
                                                                                                                                            ccccTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGAAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCa   <  1:3006763/150‑1 (MQ=255)
                                                                                                                                             cccTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAg  >  1:1040336/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAG  >  CP009273/1997170‑1997460

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