Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,067,908 T→A 65.5% V963V (GTT→GTA 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,9080TA65.5% 334.9 / 291.6 286V963V (GTT→GTAfluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (54/45);  new base A (95/92);  total (149/137)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

AGTGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCG  >  CP009273/2067759‑2068056
                                                                                                                                                     |                                                                                                                                                    
aGTGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGtt                                                                                                                                                      >  1:915906/1‑150 (MQ=255)
   gCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATc                                                                                                                                                   >  1:522471/1‑150 (MQ=255)
     acaacaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCg                                                                                                                                                 <  2:1336725/150‑1 (MQ=255)
     acaacaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCg                                                                                                                                                 <  2:4001812/150‑1 (MQ=255)
       aacaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCa                                                                                                                                               >  1:3036891/1‑150 (MQ=255)
       aacaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCa                                                                                                                                               >  2:3509323/1‑150 (MQ=255)
         caTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCAcc                                                                                                                                             <  1:3426877/150‑1 (MQ=255)
         caTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCAcc                                                                                                                                             >  1:1245210/1‑150 (MQ=255)
         caTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCAcc                                                                                                                                             >  1:273200/1‑150 (MQ=255)
          aTGTTCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACAAAGTGTGAGTGACTTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCt                                                                                                                                            >  2:1160962/1‑150 (MQ=255)
          aTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGGGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCt                                                                                                                                            <  2:2421926/150‑1 (MQ=255)
              gcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAg                                                                                                                                        <  1:2473803/150‑1 (MQ=255)
              gcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAg                                                                                                                                        <  2:514157/150‑1 (MQ=255)
               cgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGc                                                                                                                                       >  1:1156085/1‑150 (MQ=255)
                   ccGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTccc                                                                                                                                   <  1:2208179/150‑1 (MQ=255)
                   ccGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTccc                                                                                                                                   <  2:3067093/150‑1 (MQ=255)
                    cGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCg                                                                                                                                  <  2:2218821/150‑1 (MQ=255)
                        ttCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGga                                                                                                                              >  2:5005050/1‑150 (MQ=255)
                         tCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGgag                                                                                                                             >  1:2534814/1‑150 (MQ=255)
                         tCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGgag                                                                                                                             <  2:2138994/150‑1 (MQ=255)
                           cGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGat                                                                                                                           <  2:2830110/150‑1 (MQ=255)
                             tCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGatat                                                                                                                         <  2:2914977/150‑1 (MQ=255)
                             tCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGatat                                                                                                                         >  1:1271555/1‑150 (MQ=255)
                              cTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATg                                                                                                                        <  1:1052483/150‑1 (MQ=255)
                                  gCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCgtg                                                                                                                    >  2:3594944/1‑150 (MQ=255)
                                  gCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCgtg                                                                                                                    >  2:1087634/1‑150 (MQ=255)
                                    aGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGtgtg                                                                                                                  <  2:2048863/150‑1 (MQ=255)
                                        acaACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGa                                                                                                              >  2:957834/1‑150 (MQ=255)
                                         caACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGAc                                                                                                             <  2:3665716/150‑1 (MQ=255)
                                           aCGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACtt                                                                                                           <  2:5329791/150‑1 (MQ=255)
                                           aCGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACtt                                                                                                           <  2:1801212/150‑1 (MQ=255)
                                               tatGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCAc                                                                                                       <  1:1506616/150‑1 (MQ=255)
                                                   aCCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgca                                                                                                   <  1:82872/150‑1 (MQ=255)
                                                    ccTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcag                                                                                                  >  1:4660155/1‑150 (MQ=255)
                                                    ccTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcag                                                                                                  >  1:162911/1‑150 (MQ=255)
                                                      tttGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcaggc                                                                                                <  1:1380790/150‑1 (MQ=255)
                                                       ttGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcaggca                                                                                               <  2:2866672/150‑1 (MQ=255)
                                                        tGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGTGAGATATGCGTGTGGGGACTTCCACTgcaggcag                                                                                              <  2:4470378/150‑1 (MQ=255)
                                                          gCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCg                                                                                            <  2:2487331/150‑1 (MQ=255)
                                                                gCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGa                                                                                      >  2:4161156/1‑150 (MQ=255)
                                                                 cACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGAc                                                                                     <  2:2352024/150‑1 (MQ=255)
                                                                   ccTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGt                                                                                   >  2:3902233/1‑150 (MQ=255)
                                                                   ccTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGt                                                                                   >  2:1978591/1‑150 (MQ=255)
                                                                       aTCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTtctc                                                                               >  2:4252657/1‑150 (MQ=255)
                                                                           cGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCGGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTccc                                                                           >  1:3784749/1‑150 (MQ=255)
                                                                           cGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTccc                                                                           <  2:3109767/150‑1 (MQ=255)
                                                                            gTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCt                                                                          <  2:3075813/150‑1 (MQ=255)
                                                                             tGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTc                                                                         <  1:994444/150‑1 (MQ=255)
                                                                                cccccTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTcaca                                                                      >  1:3483261/1‑150 (MQ=255)
                                                                                 ccccTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAg                                                                     >  1:3111511/1‑150 (MQ=255)
                                                                                     tGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCTTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATg                                                                 >  2:881983/1‑150 (MQ=255)
                                                                                      gCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATgg                                                                <  1:15469/150‑1 (MQ=255)
                                                                                      gCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATgg                                                                <  1:3724681/150‑1 (MQ=255)
                                                                                        gTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGca                                                              >  2:3858590/1‑150 (MQ=255)
                                                                                         tGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGACTGGcac                                                             >  1:3235179/1‑150 (MQ=255)
                                                                                         atacacagCAAAACACAGTGTGAGTGAAATACCGGTGAACTGCAGGGAACAGCCTTCTGTAATCCGCACCATCAGCTCCCCGGGAGATCAGCGTGAGGGGACTTCCACTCCAGGCAGCCCGATGACGATCACTCCCGCACAGAATGGcac                                                             >  2:817013/8‑150 (MQ=255)
                                                                                          gACAGTGCAAAACACGGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGcaca                                                            <  2:4305954/150‑1 (MQ=255)
                                                                                          gACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATCTGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGcaca                                                            >  2:3204882/1‑150 (MQ=255)
                                                                                           aCAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACAt                                                           >  2:1516835/1‑150 (MQ=255)
                                                                                           aCAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACAt                                                           >  2:932221/1‑150 (MQ=255)
                                                                                             aGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCa                                                         <  1:4579655/150‑1 (MQ=255)
                                                                                              gTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCAc                                                        <  1:282117/150‑1 (MQ=255)
                                                                                              gTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCAc                                                        <  1:1771979/150‑1 (MQ=255)
                                                                                               tGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACt                                                       >  2:521512/1‑150 (MQ=255)
                                                                                                 cAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTgg                                                     <  1:362407/150‑1 (MQ=255)
                                                                                                  aaaaCACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGa                                                    >  1:753696/1‑150 (MQ=255)
                                                                                                  aaaaCACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGa                                                    >  2:2593874/1‑150 (MQ=255)
                                                                                                             tgagtgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCg                                         >  2:1315453/1‑150 (MQ=255)
                                                                                                             tgagtgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCg                                         >  2:1700673/1‑150 (MQ=255)
                                                                                                              gagtgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCgg                                        <  1:5005050/150‑1 (MQ=255)
                                                                                                                 tgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACt                                     <  2:5193242/150‑1 (MQ=255)
                                                                                                                 tgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACt                                     <  2:3427012/150‑1 (MQ=255)
                                                                                                                           gTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCgt                           <  1:4258950/150‑1 (MQ=255)
                                                                                                                           gTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCCCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCgt                           <  1:4026223/150‑1 (MQ=255)
                                                                                                                           gTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCgt                           >  2:4651913/1‑150 (MQ=255)
                                                                                                                            tGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCgtg                          <  2:3724281/150‑1 (MQ=255)
                                                                                                                              aaCTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTc                        >  2:3763700/1‑150 (MQ=255)
                                                                                                                              aaCTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTc                        <  2:4121096/150‑1 (MQ=255)
                                                                                                                               aCTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTcc                       <  1:3864284/150‑1 (MQ=255)
                                                                                                                                cTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCg                      <  2:4583497/150‑1 (MQ=255)
                                                                                                                                 tggtggGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCgg                     <  1:510746/150‑1 (MQ=255)
                                                                                                                                 tggtggGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCgg                     >  1:1992558/1‑150 (MQ=255)
                                                                                                                                  ggtggGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCggg                    <  1:855831/150‑1 (MQ=255)
                                                                                                                                    tggGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaa                  >  1:2912589/1‑150 (MQ=255)
                                                                                                                                    tggGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaa                  >  2:2989728/1‑150 (MQ=255)
                                                                                                                                       gTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAat               <  2:4228663/150‑1 (MQ=255)
                                                                                                                                         aCAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAatat             <  1:2397673/150‑1 (MQ=255)
                                                                                                                                         aCAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAatat             >  1:3544603/1‑150 (MQ=255)
                                                                                                                                         aCAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCACCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAatat             >  2:2278200/1‑150 (MQ=255)
                                                                                                                                          cAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATc            >  1:1018230/1‑150 (MQ=255)
                                                                                                                                          cAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATc            >  1:3691489/1‑150 (MQ=255)
                                                                                                                                           aGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCa           <  2:3377011/150‑1 (MQ=255)
                                                                                                                                           aGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCa           >  2:2777178/1‑150 (MQ=255)
                                                                                                                                             ccTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAcc         <  2:1296777/150‑1 (MQ=255)
                                                                                                                                              cTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGACAATATCAccc        >  2:5275653/1‑150 (MQ=255)
                                                                                                                                                tCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTg      <  1:1414/150‑1 (MQ=255)
                                                                                                                                                tCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTg      >  2:548821/1‑150 (MQ=255)
                                                                                                                                                 cTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTgg     >  2:2507151/1‑150 (MQ=255)
                                                                                                                                                  tGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTggg    >  1:27111/1‑150 (MQ=255)
                                                                                                                                                    ttATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCg  >  2:1616730/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGTGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCG  >  CP009273/2067759‑2068056

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: