Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,068,098 G→T 35.1% E1027* (GAA→TAA)  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,068,0980GT35.1% 340.6 / 205.3 224E1027* (GAA→TAA) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base G (75/70);  new base T (34/45);  total (109/115)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 229 total aligned reads displayed.

CTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCA  >  CP009273/2067949‑2068246
                                                                                                                                                     |                                                                                                                                                    
cTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTg                                                                                                                                                      >  2:3965537/1‑150 (MQ=255)
 ttCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGa                                                                                                                                                     <  1:2286897/150‑1 (MQ=255)
   ccACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAg                                                                                                                                                   <  1:3278203/150‑1 (MQ=255)
   ccACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAg                                                                                                                                                   <  2:3089929/150‑1 (MQ=255)
    cACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACCTCACTGGACCTGCAGGCCGGACTGGACGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGGTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAgg                                                                                                                                                  <  1:2384055/150‑1 (MQ=255)
       tGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGta                                                                                                                                               <  1:2879489/150‑1 (MQ=255)
           ggcagCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAAc                                                                                                                                           <  2:810043/150‑1 (MQ=255)
            gcagCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACg                                                                                                                                          >  2:1713825/1‑150 (MQ=255)
              agCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGt                                                                                                                                        >  1:354969/1‑150 (MQ=255)
              agCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGATATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGt                                                                                                                                        >  1:1933106/1‑150 (MQ=255)
                cGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCa                                                                                                                                      <  2:1311116/150‑1 (MQ=255)
                cGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCa                                                                                                                                      <  2:4142055/150‑1 (MQ=255)
                 gggATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAg                                                                                                                                     >  2:2382753/1‑150 (MQ=255)
                  ggATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAgg                                                                                                                                    >  1:1913940/1‑150 (MQ=255)
                   gATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGc                                                                                                                                   <  1:3788884/150‑1 (MQ=255)
                     tGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCca                                                                                                                                 <  1:2167636/150‑1 (MQ=255)
                     tGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCca                                                                                                                                 <  1:2167677/150‑1 (MQ=255)
                         gTTCCCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACt                                                                                                                             >  1:4031852/1‑150 (MQ=255)
                          ttCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTg                                                                                                                            <  1:152664/150‑1 (MQ=255)
                             tctcCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCGGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAAt                                                                                                                         >  2:1267940/1‑150 (MQ=255)
                               tcCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAAtgt                                                                                                                       <  2:1760656/150‑1 (MQ=255)
                               tcCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAAtgt                                                                                                                       >  1:3107689/1‑150 (MQ=255)
                               tcCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAAtgt                                                                                                                       >  1:1236654/1‑150 (MQ=255)
                                  cTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGAc                                                                                                                    >  1:444868/1‑150 (MQ=255)
                                     acaGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCtt                                                                                                                 >  2:2413049/1‑150 (MQ=255)
                                      caGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTc                                                                                                                <  1:71088/150‑1 (MQ=255)
                                         aaTGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGa                                                                                                             <  1:3055947/150‑1 (MQ=255)
                                          aTGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGAc                                                                                                            <  1:1342631/150‑1 (MQ=255)
                                          aTGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGAc                                                                                                            >  1:2433038/1‑150 (MQ=255)
                                           tGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACa                                                                                                           >  2:361349/1‑150 (MQ=255)
                                               acaTCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAAc                                                                                                       <  1:4105128/150‑1 (MQ=255)
                                               acaTCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAAc                                                                                                       >  1:3983483/1‑150 (MQ=255)
                                                  tCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCAt                                                                                                    <  2:1785554/150‑1 (MQ=255)
                                                  tCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCAt                                                                                                    <  2:2054390/150‑1 (MQ=255)
                                                  tCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCAt                                                                                                    >  1:690203/1‑150 (MQ=255)
                                                   cACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATc                                                                                                   >  1:1562854/1‑150 (MQ=255)
                                                        gACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATCACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCctc                                                                                              >  2:345000/1‑150 (MQ=255)
                                                        gACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCctc                                                                                              >  2:2125326/1‑150 (MQ=255)
                                                               aGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTgg                                                                                       <  2:3738684/150‑1 (MQ=255)
                                                               aGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTgg                                                                                       >  1:1953113/1‑150 (MQ=255)
                                                                 gCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTc                                                                                     <  1:2640738/150‑1 (MQ=255)
                                                                   cGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTccc                                                                                   <  2:4117018/150‑1 (MQ=255)
                                                                    ggACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCg                                                                                  <  2:1352342/150‑1 (MQ=255)
                                                                       cTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtc                                                                               <  1:3307063/150‑1 (MQ=255)
                                                                         ggAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcat                                                                             >  2:734553/1‑150 (MQ=255)
                                                                         ggAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcat                                                                             >  2:3334258/1‑150 (MQ=255)
                                                                         ggAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcat                                                                             <  2:2854853/150‑1 (MQ=255)
                                                                           aaGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatca                                                                           <  2:2007093/150‑1 (MQ=255)
                                                                           aaGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatca                                                                           >  1:1989405/1‑150 (MQ=255)
                                                                            aGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTcatcat                                                                          >  2:4012872/1‑150 (MQ=255)
                                                                             gCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATg                                                                         >  1:2928941/1‑150 (MQ=255)
                                                                             gCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATg                                                                         >  1:2555418/1‑150 (MQ=255)
                                                                              cccGTGTCCGGGAAAATATCACCCTGGGCTTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGa                                                                        >  2:3666388/1‑150 (MQ=255)
                                                                                     ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGAc                                                                 <  2:3186083/150‑1 (MQ=255)
                                                                                          aaaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGa                                                            <  2:2677358/150‑1 (MQ=255)
                                                                                            aaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACATTCTGACAGAACAATCGCCTCTGTGTGGTCCCGGTCATCATGACCGGGACCCGGAcc                                                          <  2:50394/150‑1 (MQ=255)
                                                                                             atatCACCCTGGGCGTTCAGGCCGGGTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTACTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCg                                                         <  2:1116666/150‑1 (MQ=255)
                                                                                              tatCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCgg                                                        >  2:4177266/1‑150 (MQ=255)
                                                                                                   cccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCaa                                                   >  1:3472117/1‑150 (MQ=255)
                                                                                                    ccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAAc                                                  >  2:1576742/1‑150 (MQ=255)
                                                                                                    ccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCAAACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAAc                                                  <  2:1867390/150‑1 (MQ=255)
                                                                                                      tGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACgg                                                >  1:964182/1‑150 (MQ=255)
                                                                                                       gggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGa                                               <  2:2510298/150‑1 (MQ=255)
                                                                                                         gCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATc                                             <  1:1676636/150‑1 (MQ=255)
                                                                                                          cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCt                                            >  2:3799135/1‑150 (MQ=255)
                                                                                                          cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCt                                            >  1:3479723/1‑150 (MQ=255)
                                                                                                           gTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCtt                                           >  2:1175682/1‑150 (MQ=255)
                                                                                                           gTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCtt                                           >  1:945157/1‑150 (MQ=255)
                                                                                                            ttCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTc                                          >  1:1258473/1‑150 (MQ=255)
                                                                                                             tCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCa                                         >  2:3479152/1‑150 (MQ=255)
                                                                                                             tCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCa                                         <  2:3434583/150‑1 (MQ=255)
                                                                                                             tCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCa                                         <  2:3937342/150‑1 (MQ=255)
                                                                                                                  ccGGTTAGGCCCACAGCGTCTGCGGCAGCTGCGCTGATGGTTATAACGGTCAGGCCACACTGACTTTGACCTTCTCACAGAACCATCGCCTCTCTGTGGCCCGGGAAATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGcc                                    <  1:4031314/150‑1 (MQ=255)
                                                                                                                     gTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCGTCAACGCcaca                                 >  1:2328297/1‑150 (MQ=255)
                                                                                                                        aTGCCCACAGCGTCAGCGGCAGAAGCGCTGATGGGTATAACGGTCAGGCCACACTGATTGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTc                              <  1:183791/150‑1 (MQ=255)
                                                                                                                         tGCCCACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCg                             <  2:2485180/150‑1 (MQ=255)
                                                                                                                           cccACAGCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCCCATTCGcg                           >  1:3347126/1‑149 (MQ=255)
                                                                                                                           cccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCt                           >  2:625336/1‑150 (MQ=255)
                                                                                                                             cacaGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTgg                         >  2:3252717/1‑150 (MQ=255)
                                                                                                                              accgcttCAGCGGCAGCCGGGTTTATGGTTATATCGTTCAGGCCACACTTAATGTGACCTTCTTACAGACCCATCGCCTCTCTGTGGCCCCGGTAATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGc                        <  1:752585/144‑1 (MQ=255)
                                                                                                                                   gTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTaa                   >  2:1926978/1‑150 (MQ=255)
                                                                                                                                    tCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAAc                  <  1:2996510/150‑1 (MQ=255)
                                                                                                                                     cAGCGGCAGCAGCTCTGAATGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTATGTGGTCCCGGTAATCATGAACGGTACCCGGACCGGCGCAACGGATCTTCAACGCCACATTAGCTGGCATTAACa                 >  2:3879820/1‑150 (MQ=255)
                                                                                                                                        cGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATa              <  2:852080/150‑1 (MQ=255)
                                                                                                                                        cGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATa              <  2:1002935/150‑1 (MQ=255)
                                                                                                                                         ggCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATaa             >  2:2997167/1‑150 (MQ=255)
                                                                                                                                           cagcagCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACa           >  2:1168326/1‑150 (MQ=255)
                                                                                                                                             gcagcGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATg         >  2:3003173/1‑150 (MQ=255)
                                                                                                                                             gcagcGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATg         >  2:913834/1‑150 (MQ=255)
                                                                                                                                              cagcGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGa        >  2:2605927/1‑150 (MQ=255)
                                                                                                                                              cagcGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGa        <  1:1362663/150‑1 (MQ=255)
                                                                                                                                              cagcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGa        <  1:1036651/150‑1 (MQ=255)
                                                                                                                                               agcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCCCCTTCGCTGGCATTAACAAAAACATGat       >  1:2473/1‑150 (MQ=255)
                                                                                                                                               agcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGat       >  2:276460/1‑150 (MQ=255)
                                                                                                                                               agcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGat       >  2:267934/1‑150 (MQ=255)
                                                                                                                                                 cgcTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGatat     >  1:578964/1‑150 (MQ=255)
                                                                                                                                                 cgcTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTATCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGatat     <  1:2245409/150‑1 (MQ=255)
                                                                                                                                                  gcTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGACCCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATAtt    <  1:2041278/150‑1 (MQ=255)
                                                                                                                                                  gcTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATAtt    <  2:2229153/150‑1 (MQ=255)
                                                                                                                                                    tGAAGTGTATACCGGACTGGCCACCCGGAAGTTGAGCTTGTTACTGATCCATCGCCTCTATGTGGTCCCGTTCTTCTTGACGGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtca  <  1:1887858/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCA  >  CP009273/2067949‑2068246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: