Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 554,327 G→T 10.5% intergenic (+117/‑103) sfmA → / → sfmC FimA homolog/putative periplasmic pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273554,3270GT10.5% 1185.9 / 73.1 335intergenic (+117/‑103)sfmA/sfmCFimA homolog/putative periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base G (154/145);  new base T (17/19);  total (171/164)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 340 total aligned reads displayed.

CTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAAT  >  CP009273/554179‑554476
                                                                                                                                                    |                                                                                                                                                     
cTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGc                                                                                                                                                      <  2:5015182/150‑1 (MQ=255)
   aCGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCttt                                                                                                                                                   <  2:1852667/150‑1 (MQ=255)
   aCGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCttt                                                                                                                                                   >  1:2072204/1‑150 (MQ=255)
     gcgACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTttat                                                                                                                                                 >  1:5142298/1‑150 (MQ=255)
      cgACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCTCTttatt                                                                                                                                                >  1:1607247/1‑150 (MQ=255)
      cgACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTttatt                                                                                                                                                <  1:4638895/150‑1 (MQ=255)
          ttttATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAATGTAATGTAACGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATAGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAAt                                                                                                                                            >  1:990604/1‑150 (MQ=255)
          ttttATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCTCTTTATTAAAt                                                                                                                                            >  2:2308303/1‑150 (MQ=255)
              attatGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGaaa                                                                                                                                        <  2:2164438/150‑1 (MQ=255)
              attatGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGaaa                                                                                                                                        <  1:3167227/150‑1 (MQ=255)
               ttatGAGATATGAATAATCAAAACCACTTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCTCTTTATTAAATGaaaa                                                                                                                                       <  2:1963283/150‑1 (MQ=255)
               ttatGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGaaaa                                                                                                                                       >  1:3619897/1‑150 (MQ=255)
                  tGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCTCTTTATTAAATGAAAAggg                                                                                                                                    >  1:3071830/1‑150 (MQ=255)
                   gagaTATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGa                                                                                                                                   <  2:5060036/150‑1 (MQ=255)
                   gagaTATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGa                                                                                                                                   >  1:506626/1‑150 (MQ=255)
                   gagaTATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGa                                                                                                                                   >  2:4473919/1‑150 (MQ=255)
                      atatGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAatg                                                                                                                                <  2:3459024/150‑1 (MQ=255)
                         tGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAatgatg                                                                                                                             >  2:1670240/1‑150 (MQ=255)
                         tGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAatgatg                                                                                                                             >  2:2835537/1‑150 (MQ=255)
                           aataatCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGAtgtg                                                                                                                           >  1:691473/1‑150 (MQ=255)
                            ataatCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATgtgt                                                                                                                          <  1:2286266/150‑1 (MQ=255)
                             gaatCAAAACCACGTTGTTTTGAATTATATATCAAGTCTTATAACAAAGTAATGTACCGTTTGTCTGAAGCTTTATGGTGGCAATGTATATCGAAATCATTTTCACTTTGTATCATGCAGATTTATTAAATGAAAAGGGAATTATGtgtt                                                                                                                         <  1:1747491/149‑1 (MQ=255)
                                tCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTa                                                                                                                      >  1:4534730/1‑150 (MQ=255)
                                 cAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTaa                                                                                                                     <  2:2714496/150‑1 (MQ=255)
                                 cAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTaa                                                                                                                     <  1:1565016/150‑1 (MQ=255)
                                 cAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTaa                                                                                                                     >  2:1700226/1‑150 (MQ=255)
                                    aaCCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGaa                                                                                                                  <  1:3727725/150‑1 (MQ=255)
                                      ccACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAAc                                                                                                                >  1:2971214/1‑150 (MQ=255)
                                      ccACGTTGATTTGACTTACATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAAc                                                                                                                <  1:4120332/150‑1 (MQ=255)
                                       cACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAAcc                                                                                                               <  2:20451/150‑1 (MQ=255)
                                           ttgttTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAg                                                                                                           <  2:926828/150‑1 (MQ=255)
                                            tgttTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGc                                                                                                          <  2:2924548/150‑1 (MQ=255)
                                               tttGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAAt                                                                                                       <  2:3101241/150‑1 (MQ=255)
                                                ttGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATc                                                                                                      <  2:4056574/150‑1 (MQ=255)
                                                  gAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCAt                                                                                                    <  2:3285715/150‑1 (MQ=255)
                                                  gAATTATATATCACGTCGTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCAt                                                                                                    <  2:4086892/150‑1 (MQ=255)
                                                    aTTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCAttt                                                                                                  <  2:63253/150‑1 (MQ=255)
                                                    aTTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCAttt                                                                                                  <  1:4788429/150‑1 (MQ=255)
                                                     ttATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTc                                                                                                 <  2:4356195/150‑1 (MQ=255)
                                                      tatataTCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCt                                                                                                <  2:796906/150‑1 (MQ=255)
                                                         atatCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTta                                                                                             <  2:4779736/150‑1 (MQ=255)
                                                         atatCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTta                                                                                             <  1:3997302/150‑1 (MQ=255)
                                                           atCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCACTCATTTCTTtata                                                                                           >  2:3997030/1‑150 (MQ=255)
                                                           atCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTtata                                                                                           <  2:1593096/150‑1 (MQ=255)
                                                           atCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTtata                                                                                           <  1:5005189/150‑1 (MQ=255)
                                                            tCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTatat                                                                                          <  1:739604/150‑1 (MQ=255)
                                                              aCGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTc                                                                                        >  2:4807907/1‑150 (MQ=255)
                                                               cGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTcc                                                                                       >  1:4667952/1‑150 (MQ=255)
                                                               cGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTcc                                                                                       >  2:3802750/1‑150 (MQ=255)
                                                                   ttATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCTATCATTTCTTTATATTCCTTAt                                                                                   <  1:2915877/150‑1 (MQ=255)
                                                                   ttATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAt                                                                                   >  1:3532789/1‑150 (MQ=255)
                                                                    tataACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAtt                                                                                  <  1:1289888/150‑1 (MQ=255)
                                                                     ataACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAttt                                                                                 >  1:4085716/1‑150 (MQ=255)
                                                                     ataACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAttt                                                                                 <  1:2147541/150‑1 (MQ=255)
                                                                      taACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAtttt                                                                                >  1:1000812/1‑150 (MQ=255)
                                                                       aaCAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAttttt                                                                               >  2:820615/1‑150 (MQ=255)
                                                                            aGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGt                                                                          >  2:1248719/1‑150 (MQ=255)
                                                                                 tGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGa                                                                     <  2:1074612/150‑1 (MQ=255)
                                                                                 tGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGa                                                                     >  1:2608663/1‑150 (MQ=255)
                                                                                  gTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAAACAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGaa                                                                    >  1:2215334/1‑150 (MQ=255)
                                                                                   tACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAAt                                                                   >  1:1989393/1‑150 (MQ=255)
                                                                                   tACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAAt                                                                   >  1:2790328/1‑150 (MQ=255)
                                                                                       ggTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATacac                                                               <  2:4225341/150‑1 (MQ=255)
                                                                                        gTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATAcaca                                                              <  2:2599670/150‑1 (MQ=255)
                                                                                            tCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGc                                                          <  1:3344880/150‑1 (MQ=255)
                                                                                            tCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCATTATTTTTGACGTCAGGAATACACAAGGc                                                          >  2:4575046/1‑150 (MQ=255)
                                                                                             cTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCg                                                         <  1:115578/150‑1 (MQ=255)
                                                                                              tGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGt                                                        >  1:3815531/1‑150 (MQ=255)
                                                                                               gAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTa                                                       >  1:1376808/1‑150 (MQ=255)
                                                                                                       aTGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTat                                               >  2:3218787/1‑150 (MQ=255)
                                                                                                        tggtggCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTatg                                              <  1:4875755/150‑1 (MQ=255)
                                                                                                        tggtggCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTatg                                              <  1:762513/150‑1 (MQ=255)
                                                                                                         ggtggCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTatga                                             >  1:912857/1‑150 (MQ=255)
                                                                                                           tggCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTatgatg                                           >  2:2789718/1‑150 (MQ=255)
                                                                                                            ggCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTAtgatga                                          >  2:1349290/1‑150 (MQ=255)
                                                                                                             gCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGAc                                         <  1:26512/150‑1 (MQ=255)
                                                                                                                aTGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTaa                                      >  1:4573615/1‑150 (MQ=255)
                                                                                                                  gTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTaaaa                                    >  2:1782201/1‑150 (MQ=255)
                                                                                                                  gTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTaaaa                                    <  2:2455920/150‑1 (MQ=255)
                                                                                                                  gTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTaaaa                                    >  2:3591260/1‑150 (MQ=255)
                                                                                                                   tAAATCTATATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAATAGGGAATGATTTGTTTTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAAt                                   <  2:1994088/150‑1 (MQ=255)
                                                                                                                   tAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAAGCGTATTAACTATGATGACTAAAAt                                   >  2:2952884/1‑150 (MQ=255)
                                                                                                                    aaaTCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATa                                  >  2:498435/1‑150 (MQ=255)
                                                                                                                          aaaTCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTa                            <  1:2808678/150‑1 (MQ=255)
                                                                                                                           aaTCATGTTCACTTTGTATCATGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTAt                           >  2:1754786/1‑150 (MQ=255)
                                                                                                                             tCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTg                         >  2:3306297/1‑150 (MQ=255)
                                                                                                                             tCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTg                         <  1:1688024/150‑1 (MQ=255)
                                                                                                                               aTGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGAt                       >  1:3068357/1‑150 (MQ=255)
                                                                                                                               aTGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGAt                       <  1:2496153/150‑1 (MQ=255)
                                                                                                                                    cACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCa                  <  2:4701349/150‑1 (MQ=255)
                                                                                                                                       tttGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATta               >  2:5069948/1‑150 (MQ=255)
                                                                                                                                        ttGTATCATGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAATGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATtat              >  2:1192361/1‑150 (MQ=255)
                                                                                                                                        ttGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATtat              <  1:1847724/150‑1 (MQ=255)
                                                                                                                                           tATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTAtatt           <  2:1397815/150‑1 (MQ=255)
                                                                                                                                            aTCATTCCGATTTATTAAATGAAAAGGTAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCTTATTAACTATGATTACTAAAATAAAGTTATTGATGCTCATTAtattt          <  2:990604/150‑1 (MQ=255)
                                                                                                                                             tCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTAtatttt         <  1:1811987/150‑1 (MQ=255)
                                                                                                                                             tCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTAtatttt         <  1:4151880/150‑1 (MQ=255)
                                                                                                                                               aTGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTAtattttat       <  2:4886482/150‑1 (MQ=255)
                                                                                                                                                 gCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGCCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTAtattttattt     >  1:4225010/1‑150 (MQ=255)
                                                                                                                                                   cGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTaa   <  1:4298639/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAAT  >  CP009273/554179‑554476

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: