Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,306,165 C→A 70.2% V125F (GTT→TTT)  rbfA ← 30s ribosome binding factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,306,1650CA70.2% 324.7 / inf 344V125F (GTT→TTT) rbfA30s ribosome binding factor
Reads supporting (aligned to +/- strand):  ref base C (51/51);  new base A (123/119);  total (174/170)
Fisher's exact test for biased strand distribution p-value = 9.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.46e-02

Only 100 of 350 total aligned reads displayed.

CACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTGCGGATGAAACCAG  >  CP009273/3306016‑3306311
                                                                                                                                                     |                                                                                                                                                  
    gTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACacga                                                                                                                                                >  1:437751/1‑150 (MQ=255)
     ttGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAacgac                                                                                                                                               >  1:3983579/1‑150 (MQ=255)
      tGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACacgacg                                                                                                                                              <  2:2798767/150‑1 (MQ=255)
        gCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGtt                                                                                                                                            >  1:3095605/1‑150 (MQ=255)
          gTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGttct                                                                                                                                          <  2:4653264/150‑1 (MQ=255)
           ttATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGttctt                                                                                                                                         >  1:3701/1‑150 (MQ=255)
            tatatataCGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGttcttc                                                                                                                                        >  2:2126261/1‑150 (MQ=255)
             atatataCGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCg                                                                                                                                       >  2:4544634/1‑150 (MQ=255)
               atataCGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTc                                                                                                                                     >  2:4334093/1‑150 (MQ=255)
                     gTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGttt                                                                                                                               >  1:4558687/1‑150 (MQ=255)
                     gTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGttt                                                                                                                               <  2:1118420/150‑1 (MQ=255)
                      tttCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTg                                                                                                                              >  2:4118849/1‑150 (MQ=255)
                       ttCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGa                                                                                                                             <  2:1983400/150‑1 (MQ=255)
                       ttCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGa                                                                                                                             >  1:3211163/1‑150 (MQ=255)
                         cACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGAcc                                                                                                                           >  2:364063/1‑150 (MQ=255)
                         cACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGAcc                                                                                                                           >  1:781040/1‑150 (MQ=255)
                          aCTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCa                                                                                                                          >  2:1669237/1‑150 (MQ=255)
                              tttGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCt                                                                                                                      <  1:958983/150‑1 (MQ=255)
                                  cAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTc                                                                                                                  >  1:264883/1‑150 (MQ=255)
                                  cAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTc                                                                                                                  <  2:1020563/150‑1 (MQ=255)
                                   aGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCa                                                                                                                 <  2:1257600/150‑1 (MQ=255)
                                    gcgcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCAc                                                                                                                <  1:3029497/150‑1 (MQ=255)
                                    gcgcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCAc                                                                                                                <  1:5053467/150‑1 (MQ=255)
                                     cgcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATGAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCAcc                                                                                                               >  1:4825070/1‑150 (MQ=255)
                                      gcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCa                                                                                                              <  1:1860621/150‑1 (MQ=255)
                                      gcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCa                                                                                                              >  2:1545019/1‑150 (MQ=255)
                                       cATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAg                                                                                                             >  1:5293090/1‑150 (MQ=255)
                                        aTCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAgg                                                                                                            >  1:2273061/1‑150 (MQ=255)
                                        aTCGGTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAgg                                                                                                            <  2:5238177/150‑1 (MQ=255)
                                          cGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGtt                                                                                                          <  1:3741732/150‑1 (MQ=255)
                                           gTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGttt                                                                                                         >  1:716701/1‑150 (MQ=255)
                                             tGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGa                                                                                                       >  1:2749272/1‑150 (MQ=255)
                                               cTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACa                                                                                                     >  1:1679168/1‑150 (MQ=255)
                                               cTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACa                                                                                                     >  1:3430889/1‑150 (MQ=255)
                                                tGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACAt                                                                                                    >  2:5315507/1‑150 (MQ=255)
                                                 ggACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATg                                                                                                   >  1:4047950/1‑150 (MQ=255)
                                                    cATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATgcgc                                                                                                >  2:1148366/1‑150 (MQ=255)
                                                     aTACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCa                                                                                               >  2:1043216/1‑150 (MQ=255)
                                                      tACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCAt                                                                                              >  2:3145833/1‑150 (MQ=255)
                                                       aCCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACTACTCATTATTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATTTTTGACCACGCTGGTCACCAGGTTTGACATTCGCATc                                                                                             >  1:4888169/1‑150 (MQ=255)
                                                       aCCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATc                                                                                             >  1:1154725/1‑150 (MQ=255)
                                                       aCCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATc                                                                                             >  1:519548/1‑150 (MQ=255)
                                                        cccTGAGGATTATCCAGCAACAAATCGACGTTAATGACCCGACAGCGACGACTAGGACGACTCATTAGTCCTCCTTGCTCTCGTCAGGGTTAAAACGACTTTCTACGACATGTTTAACAACTCTGGTCACAAGTTTAGACATACGCAtac                                                                                            >  2:3036210/1‑148 (MQ=255)
                                                           tGAGTTTTATCCAGCAACAAAACGCCGTTATTGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCtt                                                                                         <  2:1491773/150‑1 (MQ=255)
                                                           tGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCtt                                                                                         >  2:2339855/1‑150 (MQ=255)
                                                           tGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCtt                                                                                         <  2:1771952/150‑1 (MQ=255)
                                                            gAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTc                                                                                        >  1:4133823/1‑150 (MQ=255)
                                                             aGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCa                                                                                       <  2:3891701/150‑1 (MQ=255)
                                                             aGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCa                                                                                       >  1:564189/1‑150 (MQ=255)
                                                                tttATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAAcc                                                                                    <  1:1290526/150‑1 (MQ=255)
                                                                 ttATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCa                                                                                   >  2:1188498/1‑150 (MQ=255)
                                                                   aTCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCaga                                                                                 >  2:3781262/1‑150 (MQ=255)
                                                                   aTCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTACGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCaga                                                                                 >  2:3071876/1‑150 (MQ=255)
                                                                    tCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCagag                                                                                <  1:2926117/150‑1 (MQ=255)
                                                                       aGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGt                                                                             <  2:3208975/150‑1 (MQ=255)
                                                                         caacaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTg                                                                           <  1:5124405/150‑1 (MQ=255)
                                                                         caacaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTg                                                                           <  1:536902/150‑1 (MQ=255)
                                                                         caacaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTg                                                                           >  1:2087511/1‑150 (MQ=255)
                                                                          aacaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGt                                                                          >  2:3357550/1‑150 (MQ=255)
                                                                          aacaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGt                                                                          >  2:5149073/1‑150 (MQ=255)
                                                                           acaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTc                                                                         >  1:2386213/1‑150 (MQ=255)
                                                                           acaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTc                                                                         >  2:4111276/1‑150 (MQ=255)
                                                                            caaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCg                                                                        >  2:4272164/1‑150 (MQ=255)
                                                                              aaaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTa                                                                      >  1:4309464/1‑150 (MQ=255)
                                                                               aaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTag                                                                     >  1:3017521/1‑150 (MQ=255)
                                                                                aCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTaga                                                                    >  1:1419336/1‑150 (MQ=255)
                                                                                 cGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTagaa                                                                   <  2:1361858/150‑1 (MQ=255)
                                                                                  gCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTagaag                                                                  <  2:245097/150‑1 (MQ=255)
                                                                                  gCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTagaag                                                                  >  1:253022/1‑150 (MQ=255)
                                                                                     gTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGaagaag                                                               <  1:4255224/150‑1 (MQ=255)
                                                                                       tAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAAAGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGt                                                             >  2:1866976/1‑150 (MQ=255)
                                                                                        aaTGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTc                                                            <  1:1635363/150‑1 (MQ=255)
                                                                                           gTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGt                                                         <  2:4830633/150‑1 (MQ=255)
                                                                                            tCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGtt                                                        <  1:1560183/150‑1 (MQ=255)
                                                                                             cgcgaccgcgacGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTc                                                       >  2:4420950/1‑150 (MQ=255)
                                                                                              gcgaccgcgacGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTcc                                                      <  1:346157/150‑1 (MQ=255)
                                                                                              gcgaccgcgacGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTcc                                                      <  1:1806047/150‑1 (MQ=255)
                                                                                                 accgcgacGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGc                                                   <  2:4393410/150‑1 (MQ=255)
                                                                                                     cgacgacgaGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGa                                               >  2:3247958/1‑150 (MQ=255)
                                                                                                       acgacgaGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATa                                             <  2:863490/150‑1 (MQ=255)
                                                                                                         gacgaGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACg                                           >  1:3907457/1‑150 (MQ=255)
                                                                                                             aGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAgg                                       <  2:253040/150‑1 (MQ=255)
                                                                                                               gacgacTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGgcg                                     <  1:4490903/150‑1 (MQ=255)
                                                                                                                  gacTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACCTGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCAt                                  >  1:477018/1‑150 (MQ=255)
                                                                                                                  gacTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCAt                                  <  1:4918230/150‑1 (MQ=255)
                                                                                                                    cTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCg                                >  2:344228/1‑150 (MQ=255)
                                                                                                                          aGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTcc                          <  2:817522/150‑1 (MQ=255)
                                                                                                                             cctcctTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCca                       >  2:2750057/1‑150 (MQ=255)
                                                                                                                                 ctTGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCcagcag                   <  1:1507424/150‑1 (MQ=255)
                                                                                                                                  ttGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAgg                  <  1:381726/150‑1 (MQ=255)
                                                                                                                                  ttGCTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAgg                  >  2:2191357/1‑150 (MQ=255)
                                                                                                                                     cTGTCGTCCGGGTTAAAACGACGTTCTTCGTCATGTTAGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTg               >  2:5079876/1‑150 (MQ=255)
                                                                                                                                        tcgtcCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCTGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTGCgg            >  2:680318/1‑150 (MQ=255)
                                                                                                                                         cgtcCGGGTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTGCGGa           <  2:1799018/150‑1 (MQ=255)
                                                                                                                                                gTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTGCGGATGAAAcc    >  1:1966140/1‑150 (MQ=255)
                                                                                                                                                gTTAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTGCGGATGAAAcc    <  2:1406459/150‑1 (MQ=255)
                                                                                                                                                 ttAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTGCGGATGAAACCa   <  1:2138446/150‑1 (MQ=255)
                                                                                                                                                  tAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTGCGGATGAAACCAg  >  2:4048119/1‑150 (MQ=255)
                                                                                                                                                  tAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTGCGGATGAAACCAg  <  2:4750875/150‑1 (MQ=255)
                                                                                                                                                  tAAAACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTGCGGATGAAACCAg  >  2:4942723/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACCAGAGAGTTGTCGTAGAAGAAGGTCAGTTCCGGCACGATACGCAGGCGCATCGCTTTCCCCAGCAGGCTGCGGATGAAACCAG  >  CP009273/3306016‑3306311

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: