Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,357,204 G→T 5.4% E347* (GAA→TAA)  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,2040GT 5.4% 1442.8 / 31.9 361E347* (GAA→TAA) yhcDputative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base G (169/172);  new base T (12/8);  total (181/180)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 365 total aligned reads displayed.

GCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAA  >  CP009273/3357055‑3357353
                                                                                                                                                     |                                                                                                                                                     
gCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACg                                                                                                                                                       >  2:3667560/1‑150 (MQ=255)
  aaTGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGaa                                                                                                                                                     <  1:468674/150‑1 (MQ=255)
  aaTGGCAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGaa                                                                                                                                                     >  1:1480402/1‑150 (MQ=255)
    tGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATc                                                                                                                                                   >  1:4464819/1‑150 (MQ=255)
    tGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATc                                                                                                                                                   <  1:3189326/150‑1 (MQ=255)
        aaTAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCt                                                                                                                                               >  2:3968052/1‑150 (MQ=255)
               aTCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTaa                                                                                                                                        <  1:2030444/150‑1 (MQ=255)
                   tAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAg                                                                                                                                    >  1:2781953/1‑150 (MQ=255)
                    aaTGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGt                                                                                                                                   <  1:4281517/150‑1 (MQ=255)
                      tGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACTAATCACCTTGGGTAATTAGTGc                                                                                                                                 >  1:547509/1‑150 (MQ=255)
                          aCCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCAtag                                                                                                                             <  2:2677971/150‑3 (MQ=255)
                          aCCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCg                                                                                                                             <  2:4167640/150‑1 (MQ=255)
                           ccGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACTAATCACCTTGGGTAATTAGTGCATCGa                                                                                                                            <  2:1001224/150‑1 (MQ=255)
                              ttGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCg                                                                                                                         <  2:4303732/150‑1 (MQ=255)
                              ttGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCg                                                                                                                         <  2:3542995/150‑1 (MQ=255)
                               tGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCgg                                                                                                                        <  2:2299724/150‑1 (MQ=255)
                                     aaCAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGaa                                                                                                                  >  2:1416993/1‑150 (MQ=255)
                                       cAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATc                                                                                                                <  2:5013341/150‑1 (MQ=255)
                                       cAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGACTc                                                                                                                >  2:4266244/1‑150 (MQ=255)
                                        aGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCt                                                                                                               >  2:4846076/1‑150 (MQ=255)
                                        aGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCt                                                                                                               >  2:103433/1‑150 (MQ=255)
                                         gACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATAACCTTGGGTAATTAGTGCATCTAGCGGATGGAATCTg                                                                                                              >  2:4118544/1‑150 (MQ=255)
                                          aCGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGa                                                                                                             <  2:3745100/150‑1 (MQ=255)
                                           cGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACTAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGac                                                                                                            <  1:2499866/150‑1 (MQ=255)
                                            ggTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGaca                                                                                                           >  2:5109553/1‑150 (MQ=255)
                                             gTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGacac                                                                                                          <  1:1218850/150‑1 (MQ=255)
                                             gTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGacac                                                                                                          <  2:2953905/150‑1 (MQ=255)
                                                gTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCt                                                                                                       <  2:4114433/150‑1 (MQ=255)
                                                    gCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAg                                                                                                   >  1:4362162/1‑150 (MQ=255)
                                                        aaCTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCa                                                                                               <  1:1015977/150‑1 (MQ=255)
                                                          cTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAAt                                                                                             <  1:1356587/150‑1 (MQ=255)
                                                            atatTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGa                                                                                           >  1:3740610/1‑150 (MQ=255)
                                                            atatTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGa                                                                                           >  1:3158177/1‑150 (MQ=255)
                                                             tatTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGaa                                                                                          <  2:2010788/150‑1 (MQ=255)
                                                              atTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAAc                                                                                         >  2:1330622/1‑150 (MQ=255)
                                                               ttGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAAcg                                                                                        <  1:3841819/150‑1 (MQ=255)
                                                                 gTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAAcggc                                                                                      <  2:4330791/150‑1 (MQ=255)
                                                                 gTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAAcggc                                                                                      <  1:4579791/150‑1 (MQ=255)
                                                                  ttCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAAcggcg                                                                                     >  1:4621211/1‑150 (MQ=255)
                                                                   tCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAAcggcgg                                                                                    >  2:864971/1‑150 (MQ=255)
                                                                    cccTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACggcggc                                                                                   <  2:1811262/150‑1 (MQ=255)
                                                                    cccTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACggcggc                                                                                   >  2:4792651/1‑150 (MQ=255)
                                                                       tCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGtt                                                                                <  2:1081570/150‑1 (MQ=255)
                                                                       tCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGtt                                                                                <  1:2547668/150‑1 (MQ=255)
                                                                         cACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTAt                                                                              <  1:4900654/150‑1 (MQ=255)
                                                                              tGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACTAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAg                                                                         >  1:1454353/1‑150 (MQ=255)
                                                                               gTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGc                                                                        <  2:4373792/150‑1 (MQ=255)
                                                                                tttAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCa                                                                       >  1:1331632/1‑150 (MQ=255)
                                                                                tttAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCa                                                                       <  2:1784025/150‑1 (MQ=255)
                                                                                  tAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGa                                                                     >  1:945013/1‑150 (MQ=255)
                                                                                  tAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGa                                                                     >  1:2114381/1‑150 (MQ=255)
                                                                                   aaTCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGaa                                                                    <  2:3145625/150‑1 (MQ=255)
                                                                                   aaTCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGaa                                                                    >  2:3207263/1‑150 (MQ=255)
                                                                                      cAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAAt                                                                 >  2:2861411/1‑150 (MQ=255)
                                                                                         cATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTAt                                                              >  1:264542/1‑150 (MQ=255)
                                                                                          aTGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATc                                                             >  1:4997310/1‑150 (MQ=255)
                                                                                            gTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAg                                                           <  2:3347164/150‑1 (MQ=255)
                                                                                              aGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAggc                                                         >  1:687531/1‑150 (MQ=255)
                                                                                               gAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACTAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAggcg                                                        >  1:3535888/1‑150 (MQ=255)
                                                                                                 aaGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAggcggc                                                      >  1:4833876/1‑150 (MQ=255)
                                                                                                    ctgccAGGGTTTTCGCTTCGCGTTAGGGGGGGTAGACGATGACGAGGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGTCACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGTAAATTATCAGGCGGCATc                                                   <  2:3171608/147‑1 (MQ=255)
                                                                                                    cTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATc                                                   >  1:1893298/1‑150 (MQ=255)
                                                                                                    cTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATc                                                   >  2:1237738/1‑150 (MQ=255)
                                                                                                    cTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATc                                                   <  1:3932541/150‑1 (MQ=255)
                                                                                                     ttCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCa                                                  >  1:5093777/1‑150 (MQ=255)
                                                                                                     ttCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCa                                                  <  2:4693635/150‑1 (MQ=255)
                                                                                                     ttCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCa                                                  <  2:4035630/150‑1 (MQ=255)
                                                                                                      tCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCaa                                                 <  1:3702482/150‑1 (MQ=255)
                                                                                                       ccAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAAt                                                >  2:2999589/1‑150 (MQ=255)
                                                                                                         aGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATcc                                              <  2:2140119/150‑1 (MQ=255)
                                                                                                         aGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATcc                                              <  1:2402378/150‑1 (MQ=255)
                                                                                                         aGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATcc                                              <  1:4084782/150‑1 (MQ=255)
                                                                                                           ggTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCgg                                            >  1:3375484/1‑150 (MQ=255)
                                                                                                            gTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGt                                           <  1:1193633/150‑1 (MQ=255)
                                                                                                               aTCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGa                                        >  1:754122/1‑150 (MQ=255)
                                                                                                               aTCGCTGCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAAGTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGa                                        >  2:2751894/1‑150 (MQ=255)
                                                                                                                  gCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTc                                     >  2:4934189/1‑150 (MQ=255)
                                                                                                                   cTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCt                                    <  1:4235279/150‑1 (MQ=255)
                                                                                                                   cTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCt                                    <  2:1857313/150‑1 (MQ=255)
                                                                                                                         gATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTcc                              <  2:689905/150‑1 (MQ=255)
                                                                                                                         gATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTcc                              >  2:1501729/1‑150 (MQ=255)
                                                                                                                          aTAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTccc                             >  1:4157861/1‑150 (MQ=255)
                                                                                                                                cGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGccc                       <  1:2131564/150‑1 (MQ=255)
                                                                                                                                 gggTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCg                      >  2:1678453/1‑150 (MQ=255)
                                                                                                                                   gTAGACGATGACTTTGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGAt                    <  1:157267/150‑1 (MQ=255)
                                                                                                                                   gTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGAt                    <  1:3614184/150‑1 (MQ=255)
                                                                                                                                   gTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGAt                    >  1:2353502/1‑150 (MQ=255)
                                                                                                                                    tAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGAtt                   >  2:2454913/1‑150 (MQ=255)
                                                                                                                                    tAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGAtt                   <  2:4495233/150‑1 (MQ=255)
                                                                                                                                     aGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGAttt                  >  2:3434699/1‑150 (MQ=255)
                                                                                                                                      gACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTa                 <  1:3326591/150‑1 (MQ=255)
                                                                                                                                       aCGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTaa                >  1:5047934/1‑150 (MQ=255)
                                                                                                                                        cGATGACTATGACTAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAAc               <  1:2670017/150‑1 (MQ=255)
                                                                                                                                        cGATGACTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAAc               <  2:2906421/150‑1 (MQ=255)
                                                                                                                                         gATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACa              >  2:3876248/1‑150 (MQ=255)
                                                                                                                                             aCTATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGa          >  2:4227421/1‑150 (MQ=255)
                                                                                                                                               tATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGagc        <  2:616325/150‑1 (MQ=255)
                                                                                                                                               tATGACGAATCACCTTGGGTAATTAGTGCATAGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGagc        >  1:1989098/1‑150 (MQ=255)
                                                                                                                                                aTGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGagca       <  2:2321748/150‑1 (MQ=255)
                                                                                                                                                     gAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCaa  <  2:1480402/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACACGCTGGAGTGCAATGAACGGCGGCGTTATCGTAGCAGAAAATTATCAGGCGGCATCAATCCGGTCGAGTCTGGTTCCCCTGCCCGATTTAACAGTGAGCAGCCAA  >  CP009273/3357055‑3357353

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: