Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,357,229 C→A 6.0% S355* (TCG→TAG)  yhcD → putative outer membrane fimbrial subunit usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,357,2290CA 6.0% 945.0 / 21.8 247S355* (TCG→TAG) yhcDputative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base C (126/106);  new base A (10/5);  total (136/111)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 253 total aligned reads displayed.

CCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCAGGATACGAA  >  CP009273/3357082‑3357376
                                                                                                                                                   |                                                                                                                                                   
ccGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATc                                                                                                                                                     >  2:794148/1‑148 (MQ=255)
 cGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCg                                                                                                                                                    >  1:2553613/1‑148 (MQ=255)
  gTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGa                                                                                                                                                   >  1:2311380/1‑148 (MQ=255)
   ttGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAg                                                                                                                                                  >  2:429388/1‑148 (MQ=255)
    tGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGc                                                                                                                                                 <  1:1116967/148‑1 (MQ=255)
    tGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGc                                                                                                                                                 >  1:450618/1‑148 (MQ=255)
          aaCAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATgg                                                                                                                                           >  1:1313925/1‑148 (MQ=255)
          aaCAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATgg                                                                                                                                           <  1:3080051/148‑1 (MQ=255)
          aaCAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATgg                                                                                                                                           <  1:3243181/148‑1 (MQ=255)
          aaCAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATgg                                                                                                                                           >  2:2243294/1‑148 (MQ=255)
           aCAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGa                                                                                                                                          >  1:401981/1‑148 (MQ=255)
            cAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGaa                                                                                                                                         <  1:2895405/148‑1 (MQ=255)
            cAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGaa                                                                                                                                         >  2:718939/1‑148 (MQ=255)
               aCGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCt                                                                                                                                      >  2:1597466/1‑148 (MQ=255)
               aCGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCt                                                                                                                                      >  2:2901499/1‑148 (MQ=255)
                cGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTg                                                                                                                                     <  1:3212395/148‑1 (MQ=255)
                   tAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGaca                                                                                                                                  <  1:1486078/148‑1 (MQ=255)
                              aCTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGc                                                                                                                       >  2:1049767/1‑148 (MQ=255)
                              aCTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGc                                                                                                                       >  2:2602546/1‑148 (MQ=255)
                              aCTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGc                                                                                                                       <  2:2908868/148‑1 (MQ=255)
                               cTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCa                                                                                                                      <  2:564261/148‑1 (MQ=255)
                               cTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCa                                                                                                                      >  1:1045488/1‑148 (MQ=255)
                                tataTTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCaa                                                                                                                     >  2:864030/1‑148 (MQ=255)
                                tataTTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCaa                                                                                                                     >  2:863940/1‑148 (MQ=255)
                                 atatTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAAt                                                                                                                    <  2:2073866/148‑1 (MQ=255)
                                  tatTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATg                                                                                                                   >  2:3130420/1‑148 (MQ=255)
                                    ttGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGaa                                                                                                                 >  1:2562244/1‑148 (MQ=255)
                                    ttGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGaa                                                                                                                 <  2:656654/148‑1 (MQ=255)
                                    ttGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGaa                                                                                                                 >  1:2263100/1‑148 (MQ=255)
                                    ttGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGaa                                                                                                                 >  2:1318356/1‑148 (MQ=255)
                                         cccTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAAcggcg                                                                                                            >  1:3149408/1‑148 (MQ=255)
                                             ccACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGt                                                                                                        >  1:878002/1‑148 (MQ=255)
                                               aCCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTa                                                                                                      <  2:376620/148‑1 (MQ=255)
                                               aCCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTa                                                                                                      <  1:2927985/148‑1 (MQ=255)
                                                cccTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTAt                                                                                                     <  1:650935/148‑1 (MQ=255)
                                                cccTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTAt                                                                                                     <  1:625923/148‑1 (MQ=255)
                                                 ccTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATc                                                                                                    >  2:1847608/1‑148 (MQ=255)
                                                     tttAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAg                                                                                                <  1:710167/148‑1 (MQ=255)
                                                      ttAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGc                                                                                               >  1:278642/1‑148 (MQ=255)
                                                       tAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCa                                                                                              >  1:3117563/1‑148 (MQ=255)
                                                          tCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGaa                                                                                           <  2:53107/148‑1 (MQ=255)
                                                           cAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGaaa                                                                                          >  1:333939/1‑148 (MQ=255)
                                                            aGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGaaaa                                                                                         <  2:3219560/148‑1 (MQ=255)
                                                             gCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAAt                                                                                        >  1:1902561/1‑148 (MQ=255)
                                                             gCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAAt                                                                                        >  2:54990/1‑148 (MQ=255)
                                                              cATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAAtt                                                                                       >  2:281123/1‑148 (MQ=255)
                                                              cATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAAtt                                                                                       <  2:2308219/148‑1 (MQ=255)
                                                              cATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAAtt                                                                                       <  1:2666696/148‑1 (MQ=255)
                                                              cATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAAtt                                                                                       >  1:2397891/1‑148 (MQ=255)
                                                               aTGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTa                                                                                      <  1:1427362/148‑1 (MQ=255)
                                                               aTGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTa                                                                                      <  1:2259781/148‑1 (MQ=255)
                                                               aTGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTa                                                                                      <  2:1342128/148‑1 (MQ=255)
                                                                 gTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATc                                                                                    >  2:636587/1‑148 (MQ=255)
                                                                 gTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATc                                                                                    <  2:1828157/148‑1 (MQ=255)
                                                                 gTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATc                                                                                    >  2:1256976/1‑148 (MQ=255)
                                                                 gTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATc                                                                                    >  1:968902/1‑148 (MQ=255)
                                                                  tAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCa                                                                                   <  1:2832312/148‑1 (MQ=255)
                                                                  tAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCa                                                                                   <  1:2277940/148‑1 (MQ=255)
                                                                   aGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAg                                                                                  >  1:1087738/1‑148 (MQ=255)
                                                                       aGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAggcgg                                                                              >  1:2196516/1‑148 (MQ=255)
                                                                       aGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAggcgg                                                                              >  2:2363507/1‑148 (MQ=255)
                                                                       aGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAggcgg                                                                              >  2:2363397/1‑148 (MQ=255)
                                                                       aGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAAAGGCGGCGTTATCGTAGCAGAAAATTATCAggcgg                                                                              >  2:1132585/1‑148 (MQ=255)
                                                                        gCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAggcggc                                                                             >  1:2353540/1‑148 (MQ=255)
                                                                           tCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATc                                                                          >  1:661675/1‑148 (MQ=255)
                                                                             cAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCaa                                                                        <  2:1975072/148‑1 (MQ=255)
                                                                              aGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAAt                                                                       >  2:262831/1‑148 (MQ=255)
                                                                                 gTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCg                                                                    <  1:1256976/148‑1 (MQ=255)
                                                                                     tCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCg                                                                >  2:1878822/1‑148 (MQ=255)
                                                                                     tCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCg                                                                <  2:525520/148‑1 (MQ=255)
                                                                                        cTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGt                                                             >  1:544425/1‑148 (MQ=255)
                                                                                             cgATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGt                                                        >  2:1131614/1‑148 (MQ=255)
                                                                                                       ggTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGccc                                              >  2:2834221/1‑148 (MQ=255)
                                                                                                              gATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTaa                                       <  1:1913494/148‑1 (MQ=255)
                                                                                                                tGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACa                                     <  1:1345974/148‑1 (MQ=255)
                                                                                                                  aCTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGt                                   <  2:1213528/148‑1 (MQ=255)
                                                                                                                              cACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAAt                       <  1:863940/148‑1 (MQ=255)
                                                                                                                               aCCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAAtt                      <  2:2062822/148‑1 (MQ=255)
                                                                                                                               aCCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAAtt                      <  2:399074/148‑1 (MQ=255)
                                                                                                                                ccgtGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTa                     <  1:244512/145‑1 (MQ=255)
                                                                                                                                  tgtggTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGt                   <  2:906505/145‑1 (MQ=255)
                                                                                                                                   tGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTa                  <  1:1041429/148‑1 (MQ=255)
                                                                                                                                     ggTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACa                >  2:1423137/1‑148 (MQ=255)
                                                                                                                                        aaTTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCg             <  2:1707752/148‑1 (MQ=255)
                                                                                                                                         aTTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGc            >  2:2464183/1‑148 (MQ=255)
                                                                                                                                         aTTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGc            <  2:333939/148‑1 (MQ=255)
                                                                                                                                          ttAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGTAGCCAAATTAGTACATCGCa           <  2:2559635/148‑1 (MQ=255)
                                                                                                                                          ttAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCa           <  2:2263172/148‑1 (MQ=255)
                                                                                                                                          ttAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCa           <  1:1856103/148‑1 (MQ=255)
                                                                                                                                          ttAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCa           <  1:2630929/148‑1 (MQ=255)
                                                                                                                                          ttAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCa           >  1:1998545/1‑148 (MQ=255)
                                                                                                                                           tAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCAg          >  1:3258121/1‑148 (MQ=255)
                                                                                                                                            aGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCAgg         <  2:1666193/148‑1 (MQ=255)
                                                                                                                                            aGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAAGCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCAgg         <  2:3302002/148‑1 (MQ=255)
                                                                                                                                             gTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCAGGa        >  1:232704/1‑148 (MQ=255)
                                                                                                                                             gTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCAGGa        >  1:1663766/1‑148 (MQ=255)
                                                                                                                                              tGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCAGGAt       >  1:1548092/1‑148 (MQ=255)
                                                                                                                                                  tCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCAGGATACGa   <  2:563725/148‑1 (MQ=255)
                                                                                                                                                  tCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAAACCGGACGAGACTGGTTACCCTGCCCGATTTTACAGTGAGCACCCAAATTAGTACAACGCAGGATACGa   >  1:181522/1‑148 (MQ=255)
                                                                                                                                                   cGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCAGGATACGaa  >  2:316223/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAAATTAGTACATCGCAGGATACGAA  >  CP009273/3357082‑3357376

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: