Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 554,327 G→T 7.0% 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,3270GT 7.0% 697.7 / 23.5 186intergenic (+117/‑103)sfmA/sfmCFimA homolog/putative periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base G (84/89);  new base T (3/10);  total (87/99)
Fisher's exact test for biased strand distribution p-value = 8.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 189 total aligned reads displayed.

GCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTA  >  CP009273/554178‑554474
                                                                                                                                                     |                                                                                                                                                   
gCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCg                                                                                                                                                     <  1:107379/150‑1 (MQ=255)
gCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCg                                                                                                                                                     <  1:107373/150‑1 (MQ=255)
      gcgACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTttat                                                                                                                                               <  2:2294520/150‑1 (MQ=255)
          cTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTaaa                                                                                                                                           <  2:1299344/150‑1 (MQ=255)
            tttATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATg                                                                                                                                         <  1:3073154/150‑1 (MQ=255)
              tattatGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGaa                                                                                                                                       >  2:2119916/1‑150 (MQ=255)
                ttatGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGaaaa                                                                                                                                     <  1:4683472/150‑1 (MQ=255)
                ttatGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGaaaa                                                                                                                                     <  1:1431691/150‑1 (MQ=255)
                   tGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAggg                                                                                                                                  <  2:1419975/150‑1 (MQ=255)
                    gagaTATGAATCATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGa                                                                                                                                 <  1:1680990/150‑1 (MQ=255)
                       atatGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAatg                                                                                                                              <  2:2608212/150‑1 (MQ=255)
                         atGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAatgat                                                                                                                            >  1:4510492/1‑150 (MQ=255)
                          tGACTAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAatgatg                                                                                                                           <  2:4019031/150‑1 (MQ=255)
                            aataatCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGAtgtg                                                                                                                         >  1:4045283/1‑150 (MQ=255)
                            aataatCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGGTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGAtgtg                                                                                                                         >  1:933727/1‑150 (MQ=255)
                             ataatCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATgtgt                                                                                                                        <  1:1896442/150‑1 (MQ=255)
                              taatCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGtgtt                                                                                                                       <  2:1252914/150‑1 (MQ=255)
                               aatCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGtgttg                                                                                                                      >  2:4744454/1‑150 (MQ=255)
                                 tCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACGAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTa                                                                                                                    <  1:854815/150‑1 (MQ=255)
                                 tCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTa                                                                                                                    <  2:3274214/150‑1 (MQ=255)
                                  cAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTaa                                                                                                                   >  2:3040890/1‑150 (MQ=255)
                                   aaaaCCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAg                                                                                                                  <  1:2024651/150‑1 (MQ=255)
                                    aaaCCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGa                                                                                                                 >  2:2947607/1‑150 (MQ=255)
                                     aaCCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGaa                                                                                                                <  2:1041932/150‑1 (MQ=255)
                                      aCCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGaaa                                                                                                               >  1:4225112/1‑150 (MQ=255)
                                          cGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCaa                                                                                                           <  1:4199606/150‑1 (MQ=255)
                                              gttTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCa                                                                                                       >  2:2737754/1‑150 (MQ=255)
                                               ttttGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCaa                                                                                                      >  1:2854287/1‑150 (MQ=255)
                                                tttGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAAt                                                                                                     >  1:2668578/1‑150 (MQ=255)
                                                 ttGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATc                                                                                                    <  1:3457488/150‑1 (MQ=255)
                                                 ttGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATc                                                                                                    <  1:1207164/150‑1 (MQ=255)
                                                   gAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCAt                                                                                                  >  1:4424623/1‑150 (MQ=255)
                                                    aaTTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCAtt                                                                                                 <  2:4370883/150‑1 (MQ=255)
                                                       tatataTCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCt                                                                                              >  1:2572690/1‑150 (MQ=255)
                                                         tatatCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCttt                                                                                            <  1:3133494/150‑1 (MQ=255)
                                                          atatCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTta                                                                                           <  2:2814179/150‑1 (MQ=255)
                                                          atatCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTta                                                                                           >  1:928504/1‑150 (MQ=255)
                                                          atatCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTta                                                                                           >  1:928711/1‑150 (MQ=255)
                                                          atatCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACAAAAGCAATCATTTCTTta                                                                                           >  2:3116002/1‑150 (MQ=255)
                                                           tatCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTtat                                                                                          <  2:4684170/150‑1 (MQ=255)
                                                             tCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTatat                                                                                        <  2:3562053/150‑1 (MQ=255)
                                                             tCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTATTTatat                                                                                        >  1:2941845/1‑150 (MQ=255)
                                                              cACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATAtt                                                                                       <  2:3223834/150‑1 (MQ=255)
                                                              cACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATAtt                                                                                       <  2:625482/150‑1 (MQ=255)
                                                              cACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATAtt                                                                                       <  2:1652939/150‑1 (MQ=255)
                                                                 gTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCt                                                                                    >  2:1814589/1‑150 (MQ=255)
                                                                   cTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTa                                                                                  <  1:3830371/150‑1 (MQ=255)
                                                                    ttATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGTCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATTTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAt                                                                                 <  1:3860586/150‑1 (MQ=255)
                                                                     tataACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAtt                                                                                <  2:4660058/150‑1 (MQ=255)
                                                                      ataACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAttt                                                                               >  2:4242764/1‑150 (MQ=255)
                                                                       taACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAtttt                                                                              >  1:4077844/1‑150 (MQ=255)
                                                                       taACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAtttt                                                                              <  1:219666/150‑1 (MQ=255)
                                                                        aaCAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTAttttt                                                                             >  2:3395257/1‑150 (MQ=255)
                                                                           aaaGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGcc                                                                          >  1:2121509/1‑150 (MQ=255)
                                                                            aaGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCg                                                                         <  1:4649038/150‑1 (MQ=255)
                                                                            aaGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGACg                                                                         >  1:2831995/1‑150 (MQ=255)
                                                                             aGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGt                                                                        >  2:1884863/1‑150 (MQ=255)
                                                                             aGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGt                                                                        >  1:177076/1‑150 (MQ=255)
                                                                              gTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATATAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTc                                                                       >  1:2411283/1‑150 (MQ=255)
                                                                                 aTGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAgg                                                                    <  2:3687345/150‑1 (MQ=255)
                                                                                  tGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGa                                                                   <  1:1772668/150‑1 (MQ=255)
                                                                                        ggTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATacac                                                             >  1:967500/1‑150 (MQ=255)
                                                                                           tGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAg                                                          <  1:604309/150‑1 (MQ=255)
                                                                                              cTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCg                                                       <  1:1079249/150‑1 (MQ=255)
                                                                                               tGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGt                                                      >  1:773136/1‑150 (MQ=255)
                                                                                               tGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGt                                                      <  2:2660726/150‑1 (MQ=255)
                                                                                               tGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGt                                                      >  1:245672/1‑150 (MQ=255)
                                                                                               tGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGt                                                      >  1:3160642/1‑150 (MQ=255)
                                                                                                  aGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTAtt                                                   >  1:2933650/1‑150 (MQ=255)
                                                                                                      gTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACt                                               <  1:3071596/150‑1 (MQ=255)
                                                                                                       tATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTa                                              >  1:2762595/1‑150 (MQ=255)
                                                                                                        aTGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTat                                             >  2:4250051/1‑150 (MQ=255)
                                                                                                        aTGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTat                                             >  2:2911821/1‑150 (MQ=255)
                                                                                                           gtggCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTatgat                                          <  1:2710646/150‑1 (MQ=255)
                                                                                                            tggCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTatgatg                                         >  1:1178933/1‑150 (MQ=255)
                                                                                                              gCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGAc                                       <  1:4775132/150‑1 (MQ=255)
                                                                                                               cAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACt                                      <  1:2288361/150‑1 (MQ=255)
                                                                                                                aaTGTAAATCGAAATCATGTTCACTTTGTATCATGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTa                                     <  1:1123423/150‑1 (MQ=255)
                                                                                                                aaTGTAAATCGAAATCATGTTCACTTTGTATCATGCCTATTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATATACAAGGCGTATTAACTATGATGACTa                                     <  1:227484/150‑1 (MQ=255)
                                                                                                                       aTCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATaaa                              <  2:334512/150‑1 (MQ=255)
                                                                                                                          gAAATCATGTTCACTTTGTATCATGCCGCTTAAATAAATTAAAACGGAATTATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCATTATTATTGCCGTCAGGAATACACAAGGCGTATTAACAATGATGACTAAAATaaatat                           >  1:384290/1‑147 (MQ=255)
                                                                                                                           aaaTCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTa                          <  1:503411/150‑1 (MQ=255)
                                                                                                                            aaTCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTAt                         >  1:1217173/1‑150 (MQ=255)
                                                                                                                             aTCATGTTCACTTTGTATCATGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTAtt                        <  2:3441343/150‑1 (MQ=255)
                                                                                                                             aTCATGTTCACTTTGTATCATGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTAtt                        <  2:3441783/150‑1 (MQ=255)
                                                                                                                             aTCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATAcaca                                                            >  1:4042812/1‑114 (MQ=255)
                                                                                                                             aTCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATAcaca                                                            <  2:4042812/114‑1 (MQ=255)
                                                                                                                                aTGTTCACTTTGTATCATGCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGAt                     >  1:472677/1‑150 (MQ=255)
                                                                                                                                   ttCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCt                  <  1:3006071/150‑1 (MQ=255)
                                                                                                                                   ttCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCt                  >  2:543053/1‑150 (MQ=255)
                                                                                                                                     cACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCa                <  2:933727/150‑1 (MQ=255)
                                                                                                                                       cTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCAtt              <  2:4347310/150‑1 (MQ=255)
                                                                                                                                       cTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCAtt              <  2:4347190/150‑1 (MQ=255)
                                                                                                                                       cTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCAtt              <  2:1613066/150‑1 (MQ=255)
                                                                                                                                         ttGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATtat            <  2:4424902/150‑1 (MQ=255)
                                                                                                                                         ttGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATtat            >  2:1771817/1‑150 (MQ=255)
                                                                                                                                           gTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAATCAATCATTTATTTATATTCCTTATTTTTGCCGTCATGAATAAACAAGGCGGATTAACTATGATTACTAAAATAAATTTATTGATGCTCATtattt          >  1:4448569/1‑148 (MQ=255)
                                                                                                                                             aTCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTAtattt        >  2:1188112/1‑150 (MQ=255)
                                                                                                                                                  gCCTCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTAtattttattt   <  2:3773703/150‑1 (MQ=255)
                                                                                                                                                   ccGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATtttattta  <  2:3441662/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTA  >  CP009273/554178‑554474

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: