Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,996,056 A→T 6.4% 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,0560AT 6.4% 571.2 / 23.8 188I344I (ATT→ATA) fliCflagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (85/91);  new base T (10/2);  total (95/93)
Fisher's exact test for biased strand distribution p-value = 3.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 190 total aligned reads displayed.

AGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGTTGGCAGT  >  CP009273/1995909‑1996202
                                                                                                                                                   |                                                                                                                                                  
 gTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTaa                                                                                                                                                   >  1:212449/1‑148 (MQ=255)
   cAAAGCTGTGTGCAGATTACCGCCTTTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATa                                                                                                                                                 <  2:2569326/148‑1 (MQ=255)
     aaCCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGACGAGTCAGTATAGGTAATAGt                                                                                                                                               >  2:3620243/1‑148 (MQ=255)
     aaCCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGt                                                                                                                                               >  2:1730531/1‑148 (MQ=255)
     aaCCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGt                                                                                                                                               <  2:3668828/148‑1 (MQ=255)
      aCCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGtt                                                                                                                                              >  2:2917810/1‑148 (MQ=255)
       ccTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGttt                                                                                                                                             >  1:119950/1‑148 (MQ=255)
        cTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGtttt                                                                                                                                            >  2:1372481/1‑148 (MQ=255)
         tGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTACTAGTTTTa                                                                                                                                           >  2:3782552/1‑148 (MQ=255)
         tGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTa                                                                                                                                           <  1:1650777/148‑1 (MQ=255)
           tttGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAAc                                                                                                                                         >  1:1675515/1‑148 (MQ=255)
               cAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGaa                                                                                                                                     >  2:1884306/1‑148 (MQ=255)
                aGATTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGaaa                                                                                                                                    <  2:2273967/148‑1 (MQ=255)
                  aTTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACa                                                                                                                                  <  1:3236645/148‑1 (MQ=255)
                  aTTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACa                                                                                                                                  >  2:2315732/1‑148 (MQ=255)
                  aTTACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACa                                                                                                                                  >  2:396749/1‑148 (MQ=255)
                    tACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGc                                                                                                                                >  1:838464/1‑148 (MQ=255)
                    tACCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGGAATAGTTTTAACAGAAACAGc                                                                                                                                >  2:2774482/1‑148 (MQ=255)
                     aCCGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCa                                                                                                                               <  2:1351379/148‑1 (MQ=255)
                      ccGCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCAc                                                                                                                              <  1:2930779/148‑1 (MQ=255)
                        gCCATTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTTCCATCATCTACGCCCAGTTTGACCGCGGTTGGAGAACTGTCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCa                                                                                                                            >  1:288957/1‑148 (MQ=255)
                           aTTTAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCagta                                                                                                                         >  2:2119246/1‑148 (MQ=255)
                            tttAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCagtag                                                                                                                        <  2:1188090/148‑1 (MQ=255)
                             ttAAATCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCagtagt                                                                                                                       >  2:2440702/1‑148 (MQ=255)
                               aaaTCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGttt                                                                                                                     <  2:1956875/148‑1 (MQ=255)
                                 aTCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTttca                                                                                                                   >  2:3179530/1‑148 (MQ=255)
                                 aTCGGCAGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTttca                                                                                                                   <  2:2793299/148‑1 (MQ=255)
                                      cAGAATCGTATGTTTAACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTcac                                                                                                              >  1:264412/1‑148 (MQ=255)
                                       aGAATCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTcaca                                                                                                             <  2:2462547/148‑1 (MQ=255)
                                         aaTCGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATc                                                                                                           >  2:26875/1‑148 (MQ=255)
                                            cGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGc                                                                                                        >  1:1900716/1‑148 (MQ=255)
                                            cGTATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGc                                                                                                        >  2:2334297/1‑148 (MQ=255)
                                              tATGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAg                                                                                                      >  2:3337634/1‑148 (MQ=255)
                                                tGTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCg                                                                                                    <  1:3266900/148‑1 (MQ=255)
                                                 gTTTTACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGt                                                                                                   >  2:3457170/1‑148 (MQ=255)
                                                   tttACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTaa                                                                                                 <  2:1336918/148‑1 (MQ=255)
                                                     tACCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAg                                                                                               <  2:348825/148‑1 (MQ=255)
                                                      aCCATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGa                                                                                              <  2:105606/148‑1 (MQ=255)
                                                       ccATCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGAt                                                                                             <  2:2442600/148‑1 (MQ=255)
                                                          tCAATATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGc                                                                                          >  1:1222977/1‑148 (MQ=255)
                                                            aaTATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCa                                                                                        <  2:2020361/148‑1 (MQ=255)
                                                            aaTATCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCa                                                                                        <  2:2706387/148‑1 (MQ=255)
                                                             atatCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCAt                                                                                       >  2:1913571/1‑148 (MQ=255)
                                                              tatCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCAtt                                                                                      >  1:1651189/1‑148 (MQ=255)
                                                              tatCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCAtt                                                                                      >  1:1356029/1‑148 (MQ=255)
                                                               atCGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCAttt                                                                                     >  2:172122/1‑148 (MQ=255)
                                                                 cGACCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGt                                                                                   <  2:26946/148‑1 (MQ=255)
                                                                   aCCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTAt                                                                                 >  2:2094525/1‑148 (MQ=255)
                                                                   aCCACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCAGTTGTAt                                                                                 >  2:469121/1‑148 (MQ=255)
                                                                    ccACTTCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATc                                                                                >  1:688247/1‑148 (MQ=255)
                                                                         tCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTaa                                                                           <  1:2521450/148‑1 (MQ=255)
                                                                         tCTGTTTTGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTaa                                                                           >  2:1829008/1‑148 (MQ=255)
                                                                             ttttGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAgcgc                                                                       >  2:717868/1‑148 (MQ=255)
                                                                               ttGCCATCATCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCat                                                                     <  1:2222296/148‑1 (MQ=255)
                                                                                   catcatCTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGt                                                                 <  1:893530/148‑1 (MQ=255)
                                                                                    atcatcTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTa                                                                <  2:3081269/148‑1 (MQ=255)
                                                                                     tcatcTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTAt                                                               <  2:2060395/148‑1 (MQ=255)
                                                                                      catcTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATc                                                              >  2:1237107/1‑148 (MQ=255)
                                                                                      catcTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATc                                                              >  2:2526588/1‑148 (MQ=255)
                                                                                      catcTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATc                                                              <  1:2675137/148‑1 (MQ=255)
                                                                                      catcTCCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATc                                                              <  1:2674779/148‑1 (MQ=255)
                                                                                        tctcCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCgg                                                            >  2:2031108/1‑148 (MQ=255)
                                                                                         ctcCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGt                                                           <  1:3716209/148‑1 (MQ=255)
                                                                                         ctcCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGt                                                           >  1:437184/1‑148 (MQ=255)
                                                                                         ctcCGCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGt                                                           <  2:74132/148‑1 (MQ=255)
                                                                                             gCCCAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCa                                                       <  1:2794782/148‑1 (MQ=255)
                                                                                               ccAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCAtt                                                     <  1:3032963/148‑1 (MQ=255)
                                                                                                cAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTa                                                    <  2:2250283/148‑1 (MQ=255)
                                                                                                cAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTa                                                    <  1:768988/148‑1 (MQ=255)
                                                                                                cAGTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTa                                                    >  1:1987934/1‑148 (MQ=255)
                                                                                                  gTTTGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTAcc                                                  <  2:911722/148‑1 (MQ=255)
                                                                                                   tttGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTAccc                                                 >  1:583789/1‑148 (MQ=255)
                                                                                                   tttGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTAccc                                                 >  2:611842/1‑148 (MQ=255)
                                                                                                   tttGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTAccc                                                 >  2:564106/1‑148 (MQ=255)
                                                                                                    ttGACCGCGGTTGGAGAACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCt                                                <  1:2046739/148‑1 (MQ=255)
                                                                                                                  agaACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGt                                  <  1:1263824/148‑1 (MQ=255)
                                                                                                                  agaACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGt                                  <  1:86244/148‑1 (MQ=255)
                                                                                                                  agaACTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGt                                  <  1:3405901/148‑1 (MQ=255)
                                                                                                                    aaCTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGttt                                <  1:1173845/148‑1 (MQ=255)
                                                                                                                     aCTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGtttt                               <  1:200441/148‑1 (MQ=255)
                                                                                                                     aCTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGtttt                               <  2:3662687/148‑1 (MQ=255)
                                                                                                                      cTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTa                              <  2:2562847/148‑1 (MQ=255)
                                                                                                                      cTGGCGGCACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTa                              <  1:2440702/148‑1 (MQ=255)
                                                                                                                            gcACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGc                        >  2:163044/1‑148 (MQ=255)
                                                                                                                            gcACCGGAAGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGc                        <  1:2531089/148‑1 (MQ=255)
                                                                                                                                   aaGAGTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGc                 <  2:838464/148‑1 (MQ=255)
                                                                                                                                     gagTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCAc               >  1:2118583/1‑148 (MQ=255)
                                                                                                                                       gTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCa             <  2:2239965/148‑1 (MQ=255)
                                                                                                                                       gTCAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCa             <  2:8307/148‑1 (MQ=255)
                                                                                                                                         cAGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAg           >  2:3260697/1‑148 (MQ=255)
                                                                                                                                          aGTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAgg          <  1:1040959/148‑1 (MQ=255)
                                                                                                                                           gTATAGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGt         <  1:38822/148‑1 (MQ=255)
                                                                                                                                            tataGGTTATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGtt        >  1:428409/1‑148 (MQ=255)
                                                                                                                                            tataGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGtt        >  1:2841931/1‑148 (MQ=255)
                                                                                                                                            tataGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGtt        <  1:3461619/148‑1 (MQ=255)
                                                                                                                                            tataGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCAGTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGtt        >  2:1411256/1‑148 (MQ=255)
                                                                                                                                              taGGTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGTTgg      <  1:1484975/148‑1 (MQ=255)
                                                                                                                                                ggTAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGTTGGCa    >  2:1801867/1‑148 (MQ=255)
                                                                                                                                                  tAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGTTGGCAGt  >  2:1614952/1‑148 (MQ=255)
                                                                                                                                                  tAATAGTTTTAACAGAAACAGCACCAGTAGTTTCATTCACATCCGCAGCGTAAAGATTGCCATTTGTATCTTTAAGCGCATATGTATCGGTATCATTACCCTTGGAATCCTGCAGTTTTACTAAGCTAACAGCACCAAGGTTGGCAGt  <  2:2953576/148‑1 (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: ‑