Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,067,205 G→A 5.4% G729D (GGC→GAC)  flu → CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,067,2050GA 5.4% 1063.6 / 31.8 295G729D (GGC→GAC) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base G (146/132);  new base A (10/6);  total (157/138)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 297 total aligned reads displayed.

AGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTG  >  CP009273/2067057‑2067353
                                                                                                                                                    |                                                                                                                                                    
aGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGAAGGc                                                                                                                                                     >  1:3445924/1‑150 (MQ=255)
   cTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTcc                                                                                                                                                  >  1:3491408/1‑150 (MQ=255)
    tGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTccc                                                                                                                                                 >  2:4192761/1‑150 (MQ=255)
      cAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGc                                                                                                                                               >  2:868090/1‑150 (MQ=255)
       aGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCa                                                                                                                                              >  2:883905/1‑150 (MQ=255)
             gTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATc                                                                                                                                        >  2:537863/1‑150 (MQ=255)
              tGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCa                                                                                                                                       >  2:4227340/1‑150 (MQ=255)
              tGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCa                                                                                                                                       <  1:1735672/150‑1 (MQ=255)
                   ttAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCg                                                                                                                                  <  2:1074962/150‑1 (MQ=255)
                    tAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCgg                                                                                                                                 <  1:935641/150‑1 (MQ=255)
                    tAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCgg                                                                                                                                 >  1:1918084/1‑150 (MQ=255)
                      actactCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgtg                                                                                                                               <  2:4397586/150‑1 (MQ=255)
                       ctactCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgtgt                                                                                                                              >  2:1610735/1‑150 (MQ=255)
                       ctactCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgtgt                                                                                                                              <  2:1911076/150‑1 (MQ=255)
                        tactCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTa                                                                                                                             >  2:488049/1‑150 (MQ=255)
                         actCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTaa                                                                                                                            >  2:1827633/1‑150 (MQ=255)
                            cccTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATg                                                                                                                         >  2:2643023/1‑150 (MQ=255)
                            cccTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATg                                                                                                                         <  2:3329708/150‑1 (MQ=255)
                             ccTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATgg                                                                                                                        <  1:1956302/150‑1 (MQ=255)
                              cTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGt                                                                                                                       >  1:2067523/1‑150 (MQ=255)
                                cAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGa                                                                                                                     <  1:493637/150‑1 (MQ=255)
                                 aaCCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGaa                                                                                                                    >  1:1009057/1‑150 (MQ=255)
                                 aaCCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGaa                                                                                                                    <  2:2284998/150‑1 (MQ=255)
                                    cGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAaac                                                                                                                 <  2:3617075/150‑1 (MQ=255)
                                     gggACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAaaca                                                                                                                <  1:1203306/150‑1 (MQ=255)
                                      ggACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAaacaa                                                                                                               >  1:12226/1‑150 (MQ=255)
                                        aCAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAacaaca                                                                                                             <  2:4023229/150‑1 (MQ=255)
                                         cAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCTGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAg                                                                                                            <  1:2618585/150‑1 (MQ=255)
                                         cAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAg                                                                                                            >  2:2423626/1‑150 (MQ=255)
                                          aGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGc                                                                                                           <  1:4082606/150‑1 (MQ=255)
                                          aGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGACAACAACAGc                                                                                                           >  1:257398/1‑150 (MQ=255)
                                          aGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGc                                                                                                           >  1:554914/1‑150 (MQ=255)
                                            tgatgaGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgt                                                                                                         <  2:1390559/150‑1 (MQ=255)
                                            tgatgaGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgt                                                                                                         <  1:4175458/150‑1 (MQ=255)
                                             gatgaGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgtc                                                                                                        <  2:188530/150‑1 (MQ=255)
                                              atgaGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgtcc                                                                                                       <  2:3441991/150‑1 (MQ=255)
                                               tgaGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgtccg                                                                                                      <  1:3047878/150‑1 (MQ=255)
                                                gagaGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgtccgt                                                                                                     >  1:3591025/1‑150 (MQ=255)
                                                gagaGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgtccgg                                                                                                     >  2:4089376/1‑149 (MQ=255)
                                                    gCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCa                                                                                                 <  1:3965624/150‑1 (MQ=255)
                                                     cTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAg                                                                                                >  2:574101/1‑150 (MQ=255)
                                                          aTCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTc                                                                                           >  1:464602/1‑150 (MQ=255)
                                                          aTCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTc                                                                                           <  1:3677196/150‑1 (MQ=255)
                                                          aTCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTc                                                                                           >  1:2726761/1‑150 (MQ=255)
                                                           tCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCa                                                                                          <  1:2495350/150‑1 (MQ=255)
                                                              gcgcAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAggg                                                                                       <  2:636057/150‑1 (MQ=255)
                                                                 cAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTACATGCTGACACCGGCCATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGACATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCAGTCTCAGCATTCAGTGCgg                                                                                    >  2:4276320/1‑150 (MQ=255)
                                                                   gTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTc                                                                                  <  2:3171709/150‑1 (MQ=255)
                                                                      aaaaTGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCAtc                                                                               <  2:3774920/150‑1 (MQ=255)
                                                                       aaaTGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCAtct                                                                              <  1:736904/150‑1 (MQ=255)
                                                                        aaTGCGTCTCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTCTGACCGGAGGGTGGCAGGCTCCCGCAGCCATCAGCCCGGTGTAAATGGTGACAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCAtctc                                                                             <  2:755648/150‑1 (MQ=255)
                                                                         aTGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCg                                                                            >  2:3250234/1‑150 (MQ=255)
                                                                             ttATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCa                                                                        <  2:351498/150‑1 (MQ=255)
                                                                                 cGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGAt                                                                    <  2:369129/150‑1 (MQ=255)
                                                                                 cGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGAt                                                                    >  2:3166322/1‑150 (MQ=255)
                                                                                  gTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATa                                                                   >  1:3423962/1‑150 (MQ=255)
                                                                                  gTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATa                                                                   >  1:2380554/1‑150 (MQ=255)
                                                                                  gTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATa                                                                   >  1:4461262/1‑150 (MQ=255)
                                                                                   tGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCATGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGCTaa                                                                  >  2:131448/1‑150 (MQ=255)
                                                                                   tGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATaa                                                                  >  2:3819942/1‑150 (MQ=255)
                                                                                    gCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAAc                                                                 <  2:2935935/150‑1 (MQ=255)
                                                                                      aGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCTGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACaa                                                               >  1:2859180/1‑150 (MQ=255)
                                                                                      aGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACaa                                                               >  1:130935/1‑150 (MQ=255)
                                                                                      aGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACaa                                                               >  1:1472572/1‑150 (MQ=255)
                                                                                      aGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACaa                                                               >  1:598390/1‑150 (MQ=255)
                                                                                        aaGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATg                                                             <  2:2893112/150‑1 (MQ=255)
                                                                                        aaGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGGGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATCACAATg                                                             >  2:3579228/1‑150 (MQ=255)
                                                                                          gTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGc                                                           >  1:4270265/1‑150 (MQ=255)
                                                                                              cccTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTa                                                       <  1:3004260/150‑1 (MQ=255)
                                                                                               ccTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTAt                                                      <  2:3034376/150‑1 (MQ=255)
                                                                                                   tATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGcc                                                  <  1:3751418/150‑1 (MQ=255)
                                                                                                       ccTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTg                                              >  2:4143816/1‑150 (MQ=255)
                                                                                                         tCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTggg                                            <  2:4186464/150‑1 (MQ=255)
                                                                                                          ccATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTgggg                                           >  2:301379/1‑150 (MQ=255)
                                                                                                                tGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGc                                     >  1:3227499/1‑150 (MQ=255)
                                                                                                                 gACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGcc                                    <  1:3911248/150‑1 (MQ=255)
                                                                                                                    acaGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGa                                 >  2:3218591/1‑150 (MQ=255)
                                                                                                                       ggCAATGGACTATGCCCGGATTGTGGCAGGCTCCCGCAGCAATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAag                              <  1:850881/150‑1 (MQ=255)
                                                                                                                       ggCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAag                              >  2:4356467/1‑150 (MQ=255)
                                                                                                                         cAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAagca                            >  2:2343019/1‑150 (MQ=255)
                                                                                                                            tGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCg                         <  2:1598598/150‑1 (MQ=255)
                                                                                                                                cTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAg                     >  1:1082717/1‑150 (MQ=255)
                                                                                                                                cTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTTTAAATGGTTAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAg                     <  1:4193003/150‑1 (MQ=255)
                                                                                                                                 tATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCAGTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGCTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGc                    >  2:319603/1‑150 (MQ=255)
                                                                                                                                  aTGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCt                   <  1:4227326/150‑1 (MQ=255)
                                                                                                                                     aCCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATg                <  2:555393/150‑1 (MQ=255)
                                                                                                                                        ggATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTCTGGAt             >  2:2075046/1‑150 (MQ=255)
                                                                                                                                        ggATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTCTGGAt             >  1:4288730/1‑150 (MQ=255)
                                                                                                                                        ggATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGAt             >  1:4070089/1‑150 (MQ=255)
                                                                                                                                        ggATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGAt             >  1:2013511/1‑150 (MQ=255)
                                                                                                                                         gtttGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGAGTCCGTCTCCGCATTCAGGGCGGTCATCTCGGTCAAGATAACAATGGCGTTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGTCAGCTATGGAtt            <  2:1734121/148‑1 (MQ=255)
                                                                                                                                         gATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGAtt            <  2:400085/150‑1 (MQ=255)
                                                                                                                                         gATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGAtt            <  2:1734470/150‑1 (MQ=255)
                                                                                                                                          aTTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTc           <  1:2213284/150‑1 (MQ=255)
                                                                                                                                          aTTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTc           <  1:1417753/150‑1 (MQ=255)
                                                                                                                                          aTTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTc           <  1:2541153/150‑1 (MQ=255)
                                                                                                                                            tgtgGCTGGCTCCCGCAGCCATCAGACCGGTGTAACTGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgt         <  1:2302220/150‑1 (MQ=255)
                                                                                                                                              tgGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtcc       <  2:1472378/150‑1 (MQ=255)
                                                                                                                                               ggCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtccg      <  1:3300894/150‑1 (MQ=255)
                                                                                                                                                   ggCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTg  <  2:1203809/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTG  >  CP009273/2067057‑2067353

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: