Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,067,908 T→A 92.4% 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,9080TA92.4% 555.9 / 29.7 183V963V (GTT→GTAfluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (9/5);  new base A (77/91);  total (87/96)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 185 total aligned reads displayed.

TGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTG  >  CP009273/2067761‑2068052
                                                                                                                                                   |                                                                                                                                                
tGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCATAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGtt                                                                                                                                                  >  1:834679/1‑148 (MQ=255)
  cacaACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAAt                                                                                                                                                <  2:2235284/148‑1 (MQ=255)
   acaacaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATc                                                                                                                                               <  1:2011886/148‑1 (MQ=255)
    caacaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATcc                                                                                                                                              >  1:2476590/1‑148 (MQ=255)
     aacaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCg                                                                                                                                             >  1:1719708/1‑148 (MQ=255)
     aacaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCg                                                                                                                                             >  1:3458230/1‑148 (MQ=255)
     aacaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCg                                                                                                                                             <  1:3411418/148‑1 (MQ=255)
      acaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGc                                                                                                                                            <  2:2690316/148‑1 (MQ=255)
        aTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCAc                                                                                                                                          <  1:2277134/148‑1 (MQ=255)
         tgtgCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCAcc                                                                                                                                         >  2:1281938/1‑148 (MQ=255)
          gtgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCt                                                                                                                                        >  1:161388/1‑148 (MQ=255)
          gtgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCt                                                                                                                                        <  2:1962402/148‑1 (MQ=255)
          gtgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCt                                                                                                                                        <  2:187781/148‑1 (MQ=255)
           tgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGACCTGGTGGGTACAGCCTTCTGTAATCCGCACCtt                                                                                                                                       >  2:2328117/1‑148 (MQ=255)
           tgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCtt                                                                                                                                       <  1:66950/148‑1 (MQ=255)
              gtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAg                                                                                                                                    <  1:1742024/148‑1 (MQ=255)
                gcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCt                                                                                                                                  >  1:1172162/1‑148 (MQ=255)
                 ccGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTc                                                                                                                                 <  1:506917/148‑1 (MQ=255)
                 ccGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTc                                                                                                                                 <  1:1816286/148‑1 (MQ=255)
                    gTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCg                                                                                                                              <  2:1761300/148‑1 (MQ=255)
                        ccGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGga                                                                                                                          >  2:787539/1‑148 (MQ=255)
                           tCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGat                                                                                                                       <  1:1578972/148‑1 (MQ=255)
                            cTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCACAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGgagaga                                                                                                                      >  1:1985080/1‑146 (MQ=255)
                             tGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGatat                                                                                                                     <  2:3271786/148‑1 (MQ=255)
                              gggCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATg                                                                                                                    <  2:1352495/148‑1 (MQ=255)
                               ggCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGc                                                                                                                   >  2:1284318/1‑148 (MQ=255)
                               ggCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGc                                                                                                                   >  1:2973836/1‑148 (MQ=255)
                                gCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCg                                                                                                                  >  1:452910/1‑148 (MQ=255)
                                gCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCg                                                                                                                  >  1:452658/1‑148 (MQ=255)
                                gCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCg                                                                                                                  <  1:296898/148‑1 (MQ=255)
                                gCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAAAACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCg                                                                                                                  <  1:717716/148‑1 (MQ=255)
                                 cAGCCACAACGATATGACCTTTGGCGTAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCgt                                                                                                                 <  2:2044534/148‑1 (MQ=255)
                                 cAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTTAACTGGTGGGTACAGCCTTCTGTAATCCGCACATTCAGCTCCCGGGGAGATATGCgt                                                                                                                 <  2:675056/148‑1 (MQ=255)
                                 cAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCgt                                                                                                                 <  2:674975/148‑1 (MQ=255)
                                    ccACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGtgtg                                                                                                              >  2:95995/1‑148 (MQ=255)
                                       caACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTgggg                                                                                                           >  2:1907693/1‑148 (MQ=255)
                                        aaCGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGa                                                                                                          <  2:1361549/148‑1 (MQ=255)
                                        aaCGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGa                                                                                                          >  1:2846778/1‑148 (MQ=255)
                                         aCGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGAc                                                                                                         <  1:1337933/148‑1 (MQ=255)
                                          cGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACt                                                                                                        >  2:2458947/1‑148 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                    <  1:2927299/148‑1 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                    <  2:2435126/148‑1 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                    >  2:2856360/1‑148 (MQ=255)
                                                  ccTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgc                                                                                                <  1:505678/148‑1 (MQ=255)
                                                    tttGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcag                                                                                              <  1:3302351/148‑1 (MQ=255)
                                                    tttGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcag                                                                                              >  1:2095963/1‑148 (MQ=255)
                                                      tGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcaggc                                                                                            <  2:709532/148‑1 (MQ=255)
                                                            aGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCggg                                                                                      >  2:1915276/1‑148 (MQ=255)
                                                            aGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCggg                                                                                      <  2:2055894/148‑1 (MQ=255)
                                                             ggCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGa                                                                                     <  1:2635160/148‑1 (MQ=255)
                                                              gCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCGAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGAt                                                                                    >  2:3410940/1‑148 (MQ=255)
                                                              gCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGAt                                                                                    >  2:2558587/1‑148 (MQ=255)
                                                                aCCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGa                                                                                  >  2:937838/1‑148 (MQ=255)
                                                                 ccTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGAc                                                                                 >  1:349994/1‑148 (MQ=255)
                                                                   tCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGt                                                                               >  2:1513944/1‑148 (MQ=255)
                                                                     aTCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTtc                                                                             <  2:290540/148‑1 (MQ=255)
                                                                      tCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTtct                                                                            <  2:610918/148‑1 (MQ=255)
                                                                      tCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTtct                                                                            >  2:1292626/1‑148 (MQ=255)
                                                                      tCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTtct                                                                            <  1:1040396/148‑1 (MQ=255)
                                                                         cGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTtctctc                                                                         <  1:2724904/148‑1 (MQ=255)
                                                                           tGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTccc                                                                       <  1:108659/148‑1 (MQ=255)
                                                                                cccTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTcaca                                                                  >  2:2990133/1‑148 (MQ=255)
                                                                                cccTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTcaca                                                                  >  2:2990433/1‑148 (MQ=255)
                                                                                   tGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGaa                                                               <  2:3244337/148‑1 (MQ=255)
                                                                                    gCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAAt                                                              <  1:598459/148‑1 (MQ=255)
                                                                                       tGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGc                                                           >  1:2750903/1‑148 (MQ=255)
                                                                                       tGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGc                                                           >  1:2750999/1‑148 (MQ=255)
                                                                                       tGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGc                                                           >  2:1178534/1‑148 (MQ=255)
                                                                                         aCAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGcac                                                         >  2:3108445/1‑148 (MQ=255)
                                                                                            gTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATc                                                      <  2:3458230/148‑1 (MQ=255)
                                                                                             tGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCa                                                     <  1:400821/148‑1 (MQ=255)
                                                                                               cAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACt                                                   >  1:1327589/1‑148 (MQ=255)
                                                                                               cAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACt                                                   >  1:2744886/1‑148 (MQ=255)
                                                                                                       aGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGc                                           <  2:2415018/148‑1 (MQ=255)
                                                                                                        gtgtGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCa                                          <  2:1481192/148‑1 (MQ=255)
                                                                                                        gtgtGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCa                                          <  2:500897/148‑1 (MQ=255)
                                                                                                          gtgagtgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAgg                                        <  2:2098174/148‑1 (MQ=255)
                                                                                                          gtgagtgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAgg                                        >  1:167312/1‑148 (MQ=255)
                                                                                                          gtgagtgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAgg                                        >  1:2184605/1‑148 (MQ=255)
                                                                                                           tgagtgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGc                                       >  1:1115113/1‑148 (MQ=255)
                                                                                                              gtgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCgg                                    >  1:2292519/1‑148 (MQ=255)
                                                                                                                   ttACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTgg                               >  1:3077279/1‑148 (MQ=255)
                                                                                                                    tACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGa                              >  2:813978/1‑148 (MQ=255)
                                                                                                                     aCCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGaa                             >  2:2967297/1‑148 (MQ=255)
                                                                                                                        ggTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGcc                          >  2:1195018/1‑148 (MQ=255)
                                                                                                                         gTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGccc                         >  1:2538430/1‑148 (MQ=255)
                                                                                                                            aaCTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCgtg                      <  1:156929/148‑1 (MQ=255)
                                                                                                                             aCTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCgtgt                     <  1:767584/148‑1 (MQ=255)
                                                                                                                             aCTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTATCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCgtgt                     <  1:912251/148‑1 (MQ=255)
                                                                                                                              cTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTc                    <  1:3108445/148‑1 (MQ=255)
                                                                                                                              cTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTc                    >  1:2409761/1‑148 (MQ=255)
                                                                                                                                      tACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaaaa            >  2:3494542/1‑148 (MQ=255)
                                                                                                                                          gCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATc        <  2:2973836/148‑1 (MQ=255)
                                                                                                                                          gCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATc        >  2:1785621/1‑148 (MQ=255)
                                                                                                                                           ccTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCa       <  2:1702992/148‑1 (MQ=255)
                                                                                                                                               cTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCt   <  2:452658/148‑1 (MQ=255)
                                                                                                                                                tGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAcc     <  2:2290652/145‑1 (MQ=255)
                                                                                                                                                tGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTg  <  2:3506162/148‑1 (MQ=255)
                                                                                                                                                tGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTg  <  1:2558587/148‑1 (MQ=255)
                                                                                                                                                tGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTg  <  1:3410940/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTG  >  CP009273/2067761‑2068052

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: