Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,087,638 C→A 21.5% intergenic (‑288/‑26) rhaD ← / → rhaS pseudogene, truncated, rhamnulose‑phosphate aldolase/transcriptional activator of rhaBAD and rhaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,087,6380CA21.5% 494.8 / 88.5 194intergenic (‑288/‑26)rhaD/rhaSpseudogene, truncated, rhamnulose‑phosphate aldolase/transcriptional activator of rhaBAD and rhaT
Reads supporting (aligned to +/- strand):  ref base C (98/53);  new base A (26/16);  total (125/69)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

GATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGT  >  CP009273/4087491‑4087785
                                                                                                                                                   |                                                                                                                                                   
gatgTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGcc                                                                                                                                                     >  2:2479052/1‑148 (MQ=255)
    ttCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCgtgt                                                                                                                                                 >  1:2590885/1‑148 (MQ=255)
     tCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGtt                                                                                                                                                >  2:609163/1‑148 (MQ=255)
        caATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTgac                                                                                                                                             >  1:2359042/1‑148 (MQ=255)
        caATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTgac                                                                                                                                             >  1:869619/1‑148 (MQ=255)
         aaTTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTgacg                                                                                                                                            <  1:368393/148‑1 (MQ=255)
         aaTTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTgacg                                                                                                                                            >  1:944313/1‑148 (MQ=255)
           tttGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTgacgac                                                                                                                                          <  2:3311265/148‑1 (MQ=255)
            ttGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACa                                                                                                                                         >  2:847177/1‑148 (MQ=255)
              gCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATc                                                                                                                                       >  2:451942/1‑148 (MQ=255)
               cTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCa                                                                                                                                      <  1:2968656/148‑1 (MQ=255)
               cTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCa                                                                                                                                      <  1:2398885/148‑1 (MQ=255)
                  aaTTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCagga                                                                                                                                   >  1:685750/1‑148 (MQ=255)
                    ttGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCaggagg                                                                                                                                 >  2:55411/1‑148 (MQ=255)
                    ttGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCaggagg                                                                                                                                 >  2:2848241/1‑148 (MQ=255)
                     tgtgGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGc                                                                                                                                >  2:2292142/1‑148 (MQ=255)
                          tgatgtgatgCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTa                                                                                                                           <  1:2107784/148‑1 (MQ=255)
                             tgtgatgCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGa                                                                                                                        >  2:1397894/1‑148 (MQ=255)
                             tgtgatgCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGa                                                                                                                        >  2:2363009/1‑148 (MQ=255)
                              gtgatgCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGAc                                                                                                                       >  2:2403754/1‑148 (MQ=255)
                                gatgCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCg                                                                                                                     >  2:482491/1‑148 (MQ=255)
                                      cACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTAc                                                                                                               <  1:645455/148‑1 (MQ=255)
                                      cACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTAc                                                                                                               <  1:645554/148‑1 (MQ=255)
                                        ccGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACAt                                                                                                             <  2:1601325/148‑1 (MQ=255)
                                            aTTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAgtg                                                                                                         >  1:456071/1‑148 (MQ=255)
                                            aTTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAgtg                                                                                                         >  2:704815/1‑148 (MQ=255)
                                             tttCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAgtgt                                                                                                        >  2:484721/1‑148 (MQ=255)
                                                  tGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGAtt                                                                                                   >  1:1491980/1‑148 (MQ=255)
                                                   gAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGAttt                                                                                                  <  1:1757877/148‑1 (MQ=255)
                                                    aaaaTTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGAtttt                                                                                                 <  2:1226678/148‑1 (MQ=255)
                                                     aaaTTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGAttttt                                                                                                <  1:506844/148‑1 (MQ=255)
                                                      aaTTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGAtttttt                                                                                               <  2:1832946/148‑1 (MQ=255)
                                                        ttCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTc                                                                                             >  1:938707/1‑148 (MQ=255)
                                                         tCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTcc                                                                                            >  1:170786/1‑148 (MQ=255)
                                                         tCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTcc                                                                                            >  1:168740/1‑148 (MQ=255)
                                                          cACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCg                                                                                           <  2:876644/148‑1 (MQ=255)
                                                           aCGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGt                                                                                          >  1:2712799/1‑148 (MQ=255)
                                                            cGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTc                                                                                         >  1:2364424/1‑148 (MQ=255)
                                                              cTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTg                                                                                       >  1:988980/1‑148 (MQ=255)
                                                                gTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGt                                                                                     <  2:2176628/148‑1 (MQ=255)
                                                                 tATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTa                                                                                    <  1:2677392/148‑1 (MQ=255)
                                                                 tATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTa                                                                                    >  2:640694/1‑148 (MQ=255)
                                                                 tATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTa                                                                                    >  1:3262265/1‑148 (MQ=255)
                                                                  aTCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTaa                                                                                   <  2:2879512/148‑1 (MQ=255)
                                                                    cTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAAcg                                                                                 <  1:2926011/148‑1 (MQ=255)
                                                                    cTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAAcg                                                                                 >  2:2712168/1‑148 (MQ=255)
                                                                        aaaaaTCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTc                                                                             >  2:2449599/1‑148 (MQ=255)
                                                                         aaaaTCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTcc                                                                            <  1:3173564/148‑1 (MQ=255)
                                                                          aaaTCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCg                                                                           <  2:2103894/148‑1 (MQ=255)
                                                                             tCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTgg                                                                        >  1:875078/1‑148 (MQ=255)
                                                                                aCGTTTTTTACGGGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTACGTCTGGTTACGCGTCCGGGGCGa                                                                     >  1:1126099/1‑148 (MQ=255)
                                                                                   ttttttACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAg                                                                  >  1:2982881/1‑148 (MQ=255)
                                                                                            tGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGc                                                         >  2:106668/1‑148 (MQ=255)
                                                                                            tGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGc                                                         >  1:1876090/1‑148 (MQ=255)
                                                                                             ggTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGGGGCGCTAGAACCCCGGCt                                                        >  2:2293994/1‑148 (MQ=255)
                                                                                               ttttCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTcc                                                      <  2:3276330/148‑1 (MQ=255)
                                                                                               ttttCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTcc                                                      >  2:1522759/1‑148 (MQ=255)
                                                                                                tttCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTccc                                                     >  1:2645854/1‑148 (MQ=255)
                                                                                                    cGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAg                                                 >  1:2825216/1‑148 (MQ=255)
                                                                                                      tCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGc                                               >  2:2640201/1‑148 (MQ=255)
                                                                                                       cGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCg                                              >  2:682136/1‑148 (MQ=255)
                                                                                                       cGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCg                                              >  1:2114249/1‑148 (MQ=255)
                                                                                                       cGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCg                                              >  1:3213866/1‑148 (MQ=255)
                                                                                                          aaaTTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGAt                                           >  1:872316/1‑148 (MQ=255)
                                                                                                          aaaTTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGAt                                           <  1:3087418/148‑1 (MQ=255)
                                                                                                           aaTTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGAtt                                          >  1:2084474/1‑148 (MQ=255)
                                                                                                            aTTTAAGGTAAGAACCTTACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGAttt                                         <  2:45433/148‑1 (MQ=255)
                                                                                                             tttAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGAtttt                                        >  1:3043411/1‑148 (MQ=255)
                                                                                                              ttAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTc                                       >  1:605343/1‑148 (MQ=255)
                                                                                                              ttAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTc                                       <  2:1434332/148‑1 (MQ=255)
                                                                                                               tAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTcc                                      <  1:2894469/148‑1 (MQ=255)
                                                                                                                  ggTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGa                                   <  2:2126116/148‑1 (MQ=255)
                                                                                                                  ggTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGa                                   >  1:2320885/1‑148 (MQ=255)
                                                                                                                   gTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGaa                                  >  1:2188951/1‑148 (MQ=255)
                                                                                                                    tAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAc                                 >  2:1893088/1‑148 (MQ=255)
                                                                                                                    tAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAc                                 <  1:1968317/148‑1 (MQ=255)
                                                                                                                      aGAACCTGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcat                               >  2:859514/1‑148 (MQ=255)
                                                                                                                           cTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcatcatca                          >  2:3094316/1‑148 (MQ=255)
                                                                                                                           cTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcatcatca                          >  2:1204543/1‑148 (MQ=255)
                                                                                                                            tGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcatcatcat                         >  2:918996/1‑148 (MQ=255)
                                                                                                                            tGACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcatcatcat                         >  2:765028/1‑148 (MQ=255)
                                                                                                                             gACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATg                        >  1:2519145/1‑148 (MQ=255)
                                                                                                                             gACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATg                        >  1:2933173/1‑148 (MQ=255)
                                                                                                                             gACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATg                        >  1:1627256/1‑148 (MQ=255)
                                                                                                                             gACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATg                        <  1:451942/148‑1 (MQ=255)
                                                                                                                             gACCTCGTGATTACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATg                        <  1:1365821/148‑1 (MQ=255)
                                                                                                                                 tCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGAttt                    <  2:2982881/148‑1 (MQ=255)
                                                                                                                                        tACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGaa             >  2:121541/1‑148 (MQ=255)
                                                                                                                                        tACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGaa             <  1:1018255/148‑1 (MQ=255)
                                                                                                                                        tACTATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcat                               >  1:1392819/1‑130 (MQ=255)
                                                                                                                                         aCTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGaaa            <  1:466828/148‑1 (MQ=255)
                                                                                                                                           tATTTCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAAtt          >  1:1565696/1‑148 (MQ=255)
                                                                                                                                            aTTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtg         >  1:1897751/1‑148 (MQ=255)
                                                                                                                                             tttCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtgt        >  2:1648354/1‑148 (MQ=255)
                                                                                                                                             tttCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtg         <  2:481303/147‑1 (MQ=255)
                                                                                                                                              ttCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtgtg       >  2:2925419/1‑148 (MQ=255)
                                                                                                                                              ttCGCAGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtgtg       <  2:1342743/148‑1 (MQ=255)
                                                                                                                                                cGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGAt     <  1:847177/148‑1 (MQ=255)
                                                                                                                                                  ccGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTGTCCTGAACATCATCATGATTTTCATGAAATTGTGATTg   <  2:2253791/148‑1 (MQ=255)
                                                                                                                                                   cGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGt  <  1:55411/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGT  >  CP009273/4087491‑4087785

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: