Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,997,314 G→T 86.2% 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,3140GT86.2% 820.1 / 139.2 334intergenic (‑227/‑39)fliC/fliDflagellar filament structural protein (flagellin)/flagellar filament capping protein
Reads supporting (aligned to +/- strand):  ref base G (30/16);  new base T (179/109);  total (209/125)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 334 total aligned reads displayed.

TCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAG  >  CP009273/1997167‑1997460
                                                                                                                                                   |                                                                                                                                                  
   tCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTAtt                                                                                                                                               >  1:259525/1‑150 (MQ=255)
       tCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGtt                                                                                                                                           <  2:2396573/150‑1 (MQ=255)
         gTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGtttt                                                                                                                                         >  2:1748988/1‑150 (MQ=255)
         gTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGtttt                                                                                                                                         <  1:1414247/150‑1 (MQ=255)
          tAAAACAACCTTTATAATTTTTTTCAAAATCAGCCATTTTTTGTTATTAGCCGAAATACTCTTTTCTCTGCCCATTATTCCCTCTATTAAAAAAAACAATTAAACGTAAACTTTCAGAAAATCAGAACGATAAACCCTGTATTCGATTTa                                                                                                                                        >  2:3437724/1‑150 (MQ=255)
          tAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTa                                                                                                                                        >  1:4971910/1‑150 (MQ=255)
            aaaCAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACg                                                                                                                                      >  1:4489484/1‑150 (MQ=255)
            aaaCAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACg                                                                                                                                      <  1:3340000/150‑1 (MQ=255)
              acaacCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACgtg                                                                                                                                    <  2:503687/150‑1 (MQ=255)
              acaacCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACgtg                                                                                                                                    >  2:4487888/1‑150 (MQ=255)
               caacCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACgtgt                                                                                                                                   >  1:1617742/1‑150 (MQ=255)
               caacCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACgtgt                                                                                                                                   >  1:1479314/1‑150 (MQ=255)
                aacCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTc                                                                                                                                  <  2:1323619/150‑1 (MQ=255)
                  ccTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGa                                                                                                                                >  2:1507398/1‑150 (MQ=255)
                  ccTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGa                                                                                                                                <  1:2733523/150‑1 (MQ=255)
                  ccTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGa                                                                                                                                >  1:2790036/1‑150 (MQ=255)
                    tttAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGaaa                                                                                                                              >  1:4313704/1‑150 (MQ=255)
                    tttAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGaaa                                                                                                                              <  1:1172839/150‑1 (MQ=255)
                     ttAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAg                                                                                                                             >  1:3409971/1‑150 (MQ=255)
                     ttAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAg                                                                                                                             <  1:3095122/150‑1 (MQ=255)
                     ttAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAg                                                                                                                             >  1:3279065/1‑150 (MQ=255)
                        gAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATa                                                                                                                          <  1:3039660/150‑1 (MQ=255)
                         aaTTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATaa                                                                                                                         >  1:2831382/1‑150 (MQ=255)
                          aTTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATaaa                                                                                                                        <  2:3337809/150‑1 (MQ=255)
                             tttttCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAgg                                                                                                                     >  1:1756888/1‑150 (MQ=255)
                              ttttCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGa                                                                                                                    >  2:4929397/1‑150 (MQ=255)
                              ttttCAAAAACAGCCATTTTTTGTTAGTCCCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGa                                                                                                                    >  1:4403228/1‑150 (MQ=255)
                               tttCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGaa                                                                                                                   >  2:516582/1‑150 (MQ=255)
                                ttCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGaaa                                                                                                                  <  1:636209/150‑1 (MQ=255)
                                ttCAAAAAAAGCCATTTTTTTTTAGTCGCCGAAATACTATTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGaaa                                                                                                                  <  2:4091475/150‑1 (MQ=255)
                                   aaaaaCAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCg                                                                                                               >  2:1206304/1‑150 (MQ=255)
                                    aaaaCAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGc                                                                                                              >  1:4574195/1‑150 (MQ=255)
                                    aaaaCAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGc                                                                                                              >  1:2258546/1‑150 (MQ=255)
                                        cAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATgg                                                                                                          >  1:3188608/1‑150 (MQ=255)
                                        cAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATgg                                                                                                          <  2:2279391/150‑1 (MQ=255)
                                        cAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATgg                                                                                                          <  1:2928569/150‑1 (MQ=255)
                                             aTTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGt                                                                                                     >  2:4969261/1‑150 (MQ=255)
                                             aTTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGt                                                                                                     <  1:1000487/150‑1 (MQ=255)
                                              ttttttGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTa                                                                                                    >  2:4583836/1‑150 (MQ=255)
                                               tttttGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTAt                                                                                                   >  1:2336125/1‑150 (MQ=255)
                                               tttttGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTAt                                                                                                   >  2:4656910/1‑150 (MQ=255)
                                                ttttGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTAtt                                                                                                  >  1:5062236/1‑150 (MQ=255)
                                                 tttGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTAttt                                                                                                 >  2:177077/1‑150 (MQ=255)
                                                  ttGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTc                                                                                                >  1:1556042/1‑150 (MQ=255)
                                                   tGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCa                                                                                               <  2:1996914/150‑1 (MQ=255)
                                                    gTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCAt                                                                                              <  2:1413241/150‑1 (MQ=255)
                                                      tAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCg                                                                                            >  2:3041422/1‑150 (MQ=255)
                                                         tCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTg                                                                                         >  1:521363/1‑150 (MQ=255)
                                                         tCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAAGTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTg                                                                                         <  1:3591088/150‑1 (MQ=255)
                                                          cGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTgg                                                                                        >  2:251391/1‑150 (MQ=255)
                                                          cGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTgg                                                                                        >  1:1031767/1‑150 (MQ=255)
                                                          cGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTgg                                                                                        >  1:4480033/1‑150 (MQ=255)
                                                              gAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGt                                                                                    >  1:248249/1‑150 (MQ=255)
                                                               aaaTACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTc                                                                                   >  2:143736/1‑150 (MQ=255)
                                                                          tctcTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGGGTCGAAAGATAAAAGGAAATCGCATGGCACGTCTTTCATCGCTGGGAGTCGGGTCAGGTCt                                                                        >  2:3672284/1‑150 (MQ=255)
                                                                                 cccTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTa                                                                 >  1:3956557/1‑150 (MQ=255)
                                                                                  ccTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTaa                                                                <  1:3111421/150‑1 (MQ=255)
                                                                                     tATTCCCTCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGtt                                                             <  2:208639/150‑1 (MQ=255)
                                                                                        tCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCa                                                          <  2:812516/150‑1 (MQ=255)
                                                                                         cccGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCAt                                                         >  2:4795746/1‑150 (MQ=255)
                                                                                         cccGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCAt                                                         >  1:3471467/1‑150 (MQ=255)
                                                                                          ccGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATc                                                        >  1:1296115/1‑150 (MQ=255)
                                                                                           cGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAAGTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATcc                                                       <  2:273328/150‑1 (MQ=255)
                                                                                              tATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTg                                                    >  1:378491/1‑150 (MQ=255)
                                                                                               aTTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGa                                                   <  2:1670715/150‑1 (MQ=255)
                                                                                               aTTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGa                                                   <  2:4257743/150‑1 (MQ=255)
                                                                                                 tAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATa                                                 <  1:209641/150‑1 (MQ=255)
                                                                                                 tAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATa                                                 >  2:3796595/1‑150 (MQ=255)
                                                                                                 tAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATa                                                 >  1:5122267/1‑150 (MQ=255)
                                                                                                   aaaaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGc                                               >  2:4334150/1‑150 (MQ=255)
                                                                                                   aaaaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGc                                               <  1:4864202/150‑1 (MQ=255)
                                                                                                   aaaaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGc                                               >  1:2898268/1‑150 (MQ=255)
                                                                                                   aaaaaaaCAATTAAACGTAAACTTTGCGCAATACAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGc                                               >  2:4173922/1‑150 (MQ=255)
                                                                                                    aaaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGcc                                              <  2:265437/150‑1 (MQ=255)
                                                                                                     aaaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCt                                             >  1:3275039/1‑150 (MQ=255)
                                                                                                      aaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTc                                            <  2:1087126/150‑1 (MQ=255)
                                                                                                      aaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTc                                            >  1:4210439/1‑150 (MQ=255)
                                                                                                      aaaaCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTc                                            >  2:3086720/1‑150 (MQ=255)
                                                                                                         aCAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCAcc                                         >  1:806593/1‑150 (MQ=255)
                                                                                                             ttAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCAccgccg                                     <  2:4423868/150‑1 (MQ=255)
                                                                                                                aaCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCgcgc                                  >  2:3318737/1‑150 (MQ=255)
                                                                                                                aaCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCgcgc                                  >  2:3970978/1‑150 (MQ=255)
                                                                                                                aaCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCgcgc                                  >  1:2640946/1‑150 (MQ=255)
                                                                                                                 aCGTAAACTTTGCGCAATTCATACCGATAACCCCTGTATACGTTTTACGTGTCGAAAGATAAAAGGAAATCGAATGGCAAGTATTTCATCGCTGGGAGTCGGTTCAGGTCTGTATTTATGATCCATCCTTGATAGCCTCACAGCCGCGCa                                 >  2:1698992/1‑150 (MQ=255)
                                                                                                                 aCGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCa                                 >  2:147641/1‑150 (MQ=255)
                                                                                                                  cGTAAACTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCaa                                >  1:734159/1‑150 (MQ=255)
                                                                                                                    tAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCaaaa                              >  1:1860526/1‑150 (MQ=255)
                                                                                                                      aaCTTTGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAg                            >  1:2223709/1‑150 (MQ=255)
                                                                                                                          ttGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGAc                        <  2:1556146/150‑1 (MQ=255)
                                                                                                                           tGCGCAATTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTTATAGCCTCACCGCCGCGCAAAAAGCGACg                       <  2:498875/150‑1 (MQ=255)
                                                                                                                           tGCGCAATTCAGACCGATAACCCCTGTATTCGGTTTACGTGGCGAAAGATAAAAGGACATCGCATGGCAAGTAGTTCATCGCTGGGCGTCGGGTCAGGTCTGGAGGGAAGTTCCCTCCTGGATAGCCTCCCCGCCGCGCCCAAAGCGACg                       >  2:544198/1‑150 (MQ=255)
                                                                                                                                 aTTCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAAcc                 <  1:3722806/150‑1 (MQ=255)
                                                                                                                                  ttCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAAccc                >  1:1662481/1‑150 (MQ=255)
                                                                                                                                   tCAGACCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAAcccc               >  1:2812057/1‑150 (MQ=255)
                                                                                                                                       aCCGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCAtt           >  2:943393/1‑150 (MQ=255)
                                                                                                                                        ccGATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGAATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCAttt          <  2:3247577/150‑1 (MQ=255)
                                                                                                                                          taTAACCCCTGTTTTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGAAAGTATTTCAGCGCTGGGTGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCa        <  2:2406729/149‑1 (MQ=255)
                                                                                                                                          gATAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGAAAAAAGCGACGCTAACCCCCATTTCa        <  1:4074913/150‑1 (MQ=255)
                                                                                                                                            tAACCCCTGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCaaa      >  2:282046/1‑150 (MQ=255)
                                                                                                                                              aCCCCTGTATTCTTTTTAGGTGTCGAAAGATAAAAGTAAATCGCATGGAAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGTATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATc    <  1:4537193/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTATTCCCGCTATTAAAAAAAACAATTAAACGTAAACTTTGCGCAATTCAGACCGATAACCCCGGTATTCGTTTTACGTGTCGAAAGATAAAAGGAAATCGCATGGCAAGTATTTCATCGCTGGGAGTCGGGTCAGGTCTGGATTTAAGTTCCATCCTTGATAGCCTCACCGCCGCGCAAAAAGCGACGCTAACCCCCATTTCAAATCAG  >  CP009273/1997167‑1997460

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