Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,996,056 A→T 18.5% I344I (ATT→ATA fliC ← flagellar filament structural protein (flagellin)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,996,0560AT18.5% 425.0 / 91.8 179I344I (ATT→ATAfliCflagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (76/70);  new base T (11/22);  total (87/92)
Fisher's exact test for biased strand distribution p-value = 5.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

Only 100 of 182 total aligned reads displayed.

AGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGTTGGCAGT  >  CP009273/1995909‑1996202
                                                                                                                                                   |                                                                                                                                                  
 gTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTaa                                                                                                                                                   <  2:896672/148‑1 (MQ=255)
   cAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATa                                                                                                                                                 >  1:148544/1‑148 (MQ=255)
   cAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATa                                                                                                                                                 >  1:148580/1‑148 (MQ=255)
    aaaCCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAg                                                                                                                                                <  1:81539/148‑1 (MQ=255)
         tGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTa                                                                                                                                           >  2:2097137/1‑148 (MQ=255)
         tGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTa                                                                                                                                           <  2:3278705/148‑1 (MQ=255)
         tGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTa                                                                                                                                           >  1:743976/1‑148 (MQ=255)
         tGTTTGCAGATTACCGCCATGTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTa                                                                                                                                           <  2:3161603/148‑1 (MQ=255)
           tttGCAGATTACCGCCATTTAAATCGGCAGAATAGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAAc                                                                                                                                         >  2:2561392/1‑148 (MQ=255)
            ttGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACa                                                                                                                                        >  2:2047282/1‑148 (MQ=255)
               cAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGaa                                                                                                                                     <  2:622742/148‑1 (MQ=255)
                aGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGaaa                                                                                                                                    <  2:349473/148‑1 (MQ=255)
                  aTTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACa                                                                                                                                  >  2:3071966/1‑148 (MQ=255)
                  aTTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAAGATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACa                                                                                                                                  <  2:1340751/148‑1 (MQ=255)
                    tACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGc                                                                                                                                <  2:1671967/148‑1 (MQ=255)
                    tACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGc                                                                                                                                <  2:443196/148‑1 (MQ=255)
                    tACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGc                                                                                                                                >  2:1705551/1‑148 (MQ=255)
                     aCCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCa                                                                                                                               <  1:3272877/148‑1 (MQ=255)
                         ccATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCag                                                                                                                           >  2:2037365/1‑148 (MQ=255)
                           aTTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCagta                                                                                                                         >  2:1989025/1‑148 (MQ=255)
                             ttAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCagtagt                                                                                                                       >  2:132466/1‑148 (MQ=255)
                                aaTCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTGTAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTttc                                                                                                                    >  1:2297301/1‑148 (MQ=255)
                                 aTCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTttca                                                                                                                   <  2:2634242/148‑1 (MQ=255)
                                 aTCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTttca                                                                                                                   >  2:843910/1‑148 (MQ=255)
                                  tCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTttcat                                                                                                                  <  1:1874948/148‑1 (MQ=255)
                                  tCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTttcat                                                                                                                  <  1:2648991/148‑1 (MQ=255)
                                  tCGGCAGAATAGTATGTTTTACCATCAATATCGACCCCTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTTGAGAACTGGCGGCACCTGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTttcat                                                                                                                  <  1:58305/148‑1 (MQ=255)
                                   cGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTttcatt                                                                                                                 <  1:1561031/148‑1 (MQ=255)
                                    ggCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTttcattc                                                                                                                >  1:1871556/1‑148 (MQ=255)
                                    ggCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTttcattc                                                                                                                >  1:310919/1‑148 (MQ=255)
                                      cAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTcac                                                                                                              <  2:976737/148‑1 (MQ=255)
                                      cAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTcac                                                                                                              <  2:3082378/148‑1 (MQ=255)
                                       aGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTcaca                                                                                                             >  1:1323903/1‑148 (MQ=255)
                                       aGAATCGTATGTTTTACCATCAAAATCGACCACTTCTGTTTTGCCATCATCTACGCCCAGTTTGAACGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTACAAGTTTTAACAGAAACAGCACCACTAGTTTCATTcaca                                                                                                             >  1:2972499/1‑148 (MQ=255)
                                        gAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACAt                                                                                                            <  1:400086/148‑1 (MQ=255)
                                         aaTCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATc                                                                                                           >  1:1336933/1‑148 (MQ=255)
                                         aaTCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATc                                                                                                           >  1:146403/1‑148 (MQ=255)
                                         aaTCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATc                                                                                                           >  2:691183/1‑148 (MQ=255)
                                         aaTCGTATGTTTTACCATCAATATAGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATc                                                                                                           >  1:20831/1‑148 (MQ=255)
                                          aTCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATcc                                                                                                          >  1:2415823/1‑148 (MQ=255)
                                            cGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTAACATCCGc                                                                                                        >  1:632694/1‑148 (MQ=255)
                                              tATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAg                                                                                                      <  2:1457492/148‑1 (MQ=255)
                                               aTGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGc                                                                                                     <  2:117832/148‑1 (MQ=255)
                                                  ttttACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTa                                                                                                  >  2:1589018/1‑148 (MQ=255)
                                                   tttACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTaa                                                                                                 <  1:672148/148‑1 (MQ=255)
                                                    ttACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTaaa                                                                                                <  2:1232871/148‑1 (MQ=255)
                                                     tACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAg                                                                                               <  2:2144986/148‑1 (MQ=255)
                                                     tACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAg                                                                                               <  2:1310307/148‑1 (MQ=255)
                                                       ccATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAATGAt                                                                                             >  1:2017096/1‑148 (MQ=255)
                                                         aTCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTg                                                                                           <  1:2044072/148‑1 (MQ=255)
                                                         aTCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTg                                                                                           <  1:269965/148‑1 (MQ=255)
                                                          tCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGc                                                                                          >  2:2606268/1‑148 (MQ=255)
                                                           cAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGcc                                                                                         <  2:182443/148‑1 (MQ=255)
                                                            aaTATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCa                                                                                        <  2:191263/148‑1 (MQ=255)
                                                             atatCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCAt                                                                                       >  1:650760/1‑148 (MQ=255)
                                                              tatCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCAtt                                                                                      >  1:2728380/1‑148 (MQ=255)
                                                               atCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCAttt                                                                                     <  1:2791214/148‑1 (MQ=255)
                                                                tCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTg                                                                                    >  2:77845/1‑148 (MQ=255)
                                                                 cGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGt                                                                                   >  2:2005138/1‑148 (MQ=255)
                                                                  gACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTa                                                                                  >  1:941311/1‑148 (MQ=255)
                                                                    ccACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATc                                                                                >  2:1047454/1‑148 (MQ=255)
                                                                     cACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCt                                                                               <  2:878463/148‑1 (MQ=255)
                                                                     cACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCt                                                                               <  2:38312/148‑1 (MQ=255)
                                                                        ttCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTa                                                                            >  1:1708384/1‑148 (MQ=255)
                                                                          cTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAg                                                                          >  1:1681367/1‑148 (MQ=255)
                                                                                tGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCata                                                                    >  2:2369765/1‑148 (MQ=255)
                                                                                  ccATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATg                                                                  <  1:2215182/148‑1 (MQ=255)
                                                                                  ccATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATg                                                                  >  1:905691/1‑148 (MQ=255)
                                                                                   catcatCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGt                                                                 <  1:2109793/148‑1 (MQ=255)
                                                                                     tcatcTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGGTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTAt                                                               >  2:139401/1‑148 (MQ=255)
                                                                                      catcTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATc                                                              <  2:310919/148‑1 (MQ=255)
                                                                                       atcTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCg                                                             <  1:2523846/148‑1 (MQ=255)
                                                                                         ctcCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGt                                                           <  1:2097137/148‑1 (MQ=255)
                                                                                            cGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATc                                                        <  1:2898921/148‑1 (MQ=255)
                                                                                             gCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCa                                                       <  1:2561392/148‑1 (MQ=255)
                                                                                              cccAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAgcgc                                                                       <  2:2246026/131‑1 (MQ=255)
                                                                                                   tttGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTAccc                                                 <  2:676258/148‑1 (MQ=255)
                                                                                                   tttGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTAccc                                                 <  2:675770/148‑1 (MQ=255)
                                                                                                       aCCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTgg                                             >  1:3123994/1‑148 (MQ=255)
                                                                                                        ccGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGa                                            >  1:1410859/1‑148 (MQ=255)
                                                                                                         cgcgGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGaa                                           <  1:341547/148‑1 (MQ=255)
                                                                                                           cgGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATc                                         <  1:90517/148‑1 (MQ=255)
                                                                                                           cgGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAAGAGATTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATc                                         <  2:2517389/148‑1 (MQ=255)
                                                                                                              ttGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTg                                      >  1:1422458/1‑148 (MQ=255)
                                                                                                                 gagaACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAg                                   >  2:2929264/1‑148 (MQ=255)
                                                                                                                   gaACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGtt                                 <  2:3210956/148‑1 (MQ=255)
                                                                                                                              aCCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTa                      <  2:2802035/148‑1 (MQ=255)
                                                                                                                                cGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAAc                    <  1:86495/148‑1 (MQ=255)
                                                                                                                                  gAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAg                  >  1:2850175/1‑148 (MQ=255)
                                                                                                                                  gAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAg                  >  1:1157110/1‑148 (MQ=255)
                                                                                                                                  gAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAg                  >  1:1157216/1‑148 (MQ=255)
                                                                                                                                  gAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTAACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCAATACCCTTGGAATCCAGCAGTTTTACTAAGCTAACAg                  >  1:2991/1‑148 (MQ=255)
                                                                                                                                   aaGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGc                 >  1:2794706/1‑148 (MQ=255)
                                                                                                                                         cAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAg           >  2:2963784/1‑148 (MQ=255)
                                                                                                                                           gTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGt         <  2:2311973/148‑1 (MQ=255)
                                                                                                                                           gTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGt         <  2:1336933/148‑1 (MQ=255)
                                                                                                                                            tataGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGtt        >  2:1491926/1‑148 (MQ=255)
                                                                                                                                             ataGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGTTg       >  2:107747/1‑148 (MQ=255)
                                                                                                                                                ggTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGTTGGCa    <  1:2606268/148‑1 (MQ=255)
                                                                                                                                                 gTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATACTGCAGTTTTAATAAGCTAACAGCACCAAGGTTGGCAg   >  1:2149910/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGTTGGCAGT  >  CP009273/1995909‑1996202

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