Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,088,850 C→A 100% R93S (CGT→AGT)  rhaR → transcriptional activator of rhaSR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,088,8500CA97.5% 625.8 / 6.0 198R93S (CGT→AGT) rhaRtranscriptional activator of rhaSR
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (110/83);  total (113/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 198 total aligned reads displayed.

CGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGC  >  CP009273/4088709‑4088991
                                                                                                                                             |                                                                                                                                             
cgcgGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGa                                                                                                                                         <  1:2711917/148‑1 (MQ=255)
  cgGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAg                                                                                                                                       <  1:1473564/148‑1 (MQ=255)
   ggTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGc                                                                                                                                      >  1:2924724/1‑148 (MQ=255)
      aaTGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGa                                                                                                                                   >  2:1205238/1‑148 (MQ=255)
      aaTGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGa                                                                                                                                   <  1:1659433/148‑1 (MQ=255)
       aTGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGaa                                                                                                                                  >  2:2866935/1‑148 (MQ=255)
        tGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAAt                                                                                                                                 >  1:1721169/1‑148 (MQ=255)
        tGGCCTGCATGTACTCAACGATCGCCCTTATCGCATAACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATGTATTGCCCGGAGAGTCTGAAGCTGAAt                                                                                                                                 >  1:1972152/1‑148 (MQ=255)
         ggCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATc                                                                                                                                >  1:869388/1‑148 (MQ=255)
         ggCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATc                                                                                                                                <  1:2050626/148‑1 (MQ=255)
         ggCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATc                                                                                                                                <  1:3240809/148‑1 (MQ=255)
               cATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACt                                                                                                                          >  2:3514114/1‑148 (MQ=255)
                 tGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTgg                                                                                                                        <  1:363133/148‑1 (MQ=255)
                 tGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTgg                                                                                                                        >  1:365808/1‑148 (MQ=255)
                          cGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCg                                                                                                               >  1:1088536/1‑148 (MQ=255)
                            aTCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGAt                                                                                                             >  1:2863177/1‑148 (MQ=255)
                             tCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGAtt                                                                                                            >  2:1922678/1‑148 (MQ=255)
                             tCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGAtt                                                                                                            >  2:3300926/1‑148 (MQ=255)
                              cGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTc                                                                                                           <  2:3201490/148‑1 (MQ=255)
                              cGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTc                                                                                                           >  1:1354139/1‑148 (MQ=255)
                                cccTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCg                                                                                                         <  1:3341863/148‑1 (MQ=255)
                                cccTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCg                                                                                                         <  2:1816433/148‑1 (MQ=255)
                                cccTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCg                                                                                                         >  2:956437/1‑148 (MQ=255)
                                 ccTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCgg                                                                                                        >  2:3516023/1‑148 (MQ=255)
                                  cTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCggg                                                                                                       >  1:2971315/1‑148 (MQ=255)
                                  cTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCggg                                                                                                       <  1:2539372/148‑1 (MQ=255)
                                  cTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCggg                                                                                                       <  2:497430/148‑1 (MQ=255)
                                    tATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGAt                                                                                                     >  1:1727703/1‑148 (MQ=255)
                                       cGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTa                                                                                                  <  1:550007/148‑1 (MQ=255)
                                       cGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTa                                                                                                  <  2:1898987/148‑1 (MQ=255)
                                       cGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTa                                                                                                  <  2:101357/148‑1 (MQ=255)
                                           ttACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGc                                                                                              <  2:2967779/148‑1 (MQ=255)
                                           ttACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGc                                                                                              <  1:472293/148‑1 (MQ=255)
                                                  tGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCa                                                                                       >  2:3242352/1‑148 (MQ=255)
                                                   ggCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAg                                                                                      >  2:2689635/1‑148 (MQ=255)
                                                     cGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAggg                                                                                    >  1:1000782/1‑148 (MQ=255)
                                                     cGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAggg                                                                                    >  1:1000627/1‑148 (MQ=255)
                                                     cGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAggg                                                                                    <  1:3422013/148‑1 (MQ=255)
                                                     cGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAggg                                                                                    >  1:2633044/1‑148 (MQ=255)
                                                        tctcTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCaa                                                                                 >  2:3114509/1‑148 (MQ=255)
                                                          tctTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAAcc                                                                               >  2:3177163/1‑148 (MQ=255)
                                                              ttACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCacac                                                                           >  2:1367409/1‑148 (MQ=255)
                                                                   ttCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgcg                                                                      <  2:881658/148‑1 (MQ=255)
                                                                   ttCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgcg                                                                      <  2:1464814/148‑1 (MQ=255)
                                                                     cATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCt                                                                    <  2:2066958/148‑1 (MQ=255)
                                                                     cATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCt                                                                    >  1:1920122/1‑148 (MQ=255)
                                                                        gCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAg                                                                 >  2:540264/1‑148 (MQ=255)
                                                                         cTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAgg                                                                <  1:3513352/148‑1 (MQ=255)
                                                                         cTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAgg                                                                <  2:617264/148‑1 (MQ=255)
                                                                            aCGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAg                                                             >  1:2899763/1‑148 (MQ=255)
                                                                              gATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCa                                                           <  1:313869/148‑1 (MQ=255)
                                                                                tAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGNATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATg                                                         >  1:1802415/1‑148 (MQ=255)
                                                                                  aaCACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATggg                                                       <  2:309865/148‑1 (MQ=255)
                                                                                  aaCACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATggg                                                       >  1:2201675/1‑148 (MQ=255)
                                                                                   acacTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCGTAGGTAGCATgggg                                                      >  2:905653/1‑148 (MQ=255)
                                                                                    cacTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGa                                                     <  2:3429223/148‑1 (MQ=255)
                                                                                    cacTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGa                                                     <  1:1877425/148‑1 (MQ=255)
                                                                                      cTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATg                                                   >  1:3173884/1‑148 (MQ=255)
                                                                                      cTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTAGTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATg                                                   <  1:942794/148‑1 (MQ=255)
                                                                                       tCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATgg                                                  >  2:1614692/1‑148 (MQ=255)
                                                                                          tACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGgcgc                                               <  1:849678/148‑1 (MQ=255)
                                                                                               ttCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGgcg                                          <  1:1127200/148‑1 (MQ=255)
                                                                                                  cGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCgg                                       <  2:2728705/148‑1 (MQ=255)
                                                                                                    ttAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCa                                     >  1:954219/1‑148 (MQ=255)
                                                                                                     tAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAg                                    >  2:15450/1‑148 (MQ=255)
                                                                                                       aCGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGt                                  <  2:2425730/148‑1 (MQ=255)
                                                                                                        cGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGtt                                 <  2:2578404/148‑1 (MQ=255)
                                                                                                        cGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGtt                                 <  1:2455252/148‑1 (MQ=255)
                                                                                                           tCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATc                              >  2:441166/1‑148 (MQ=255)
                                                                                                           tCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATc                              >  2:1651610/1‑148 (MQ=255)
                                                                                                           tCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATc                              >  1:3541002/1‑148 (MQ=255)
                                                                                                            cTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCg                             >  1:1680209/1‑148 (MQ=255)
                                                                                                              ggTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGACGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGt                           <  2:484380/148‑1 (MQ=255)
                                                                                                              ggTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGt                           <  1:2984459/148‑1 (MQ=255)
                                                                                                               gTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTc                          >  1:1702416/1‑148 (MQ=255)
                                                                                                                 tttGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAg                        <  2:2167158/148‑1 (MQ=255)
                                                                                                                   tGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCt                      <  1:585066/148‑1 (MQ=255)
                                                                                                                   tGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCt                      <  2:1353514/148‑1 (MQ=255)
                                                                                                                     cAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTg                    >  2:770010/1‑148 (MQ=255)
                                                                                                                     cAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTg                    <  1:531861/148‑1 (MQ=255)
                                                                                                                     cAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTg                    >  1:1400379/1‑148 (MQ=255)
                                                                                                                       gAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAg                  >  1:131725/1‑148 (MQ=255)
                                                                                                                        aaTATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGc                 >  2:2347116/1‑148 (MQ=255)
                                                                                                                         atatTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCa                <  1:2117299/148‑1 (MQ=255)
                                                                                                                           attattTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATg              >  2:2538423/1‑148 (MQ=255)
                                                                                                                           attattTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCAGGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATg              >  1:3408147/1‑148 (MQ=255)
                                                                                                                           attattTATTGCCAGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGACAGCGCAGCGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATg              >  2:2538441/1‑148 (MQ=255)
                                                                                                                            ttatttatTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGa             >  2:1750122/1‑148 (MQ=255)
                                                                                                                            ttatttatTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGa             >  2:3507978/1‑148 (MQ=255)
                                                                                                                              atttattGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGaaa           >  2:2080113/1‑148 (MQ=255)
                                                                                                                               tttattGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAag          <  2:1088536/148‑1 (MQ=255)
                                                                                                                                ttattGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagt         >  1:1702688/1‑148 (MQ=255)
                                                                                                                                ttattGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagt         >  2:1304847/1‑148 (MQ=255)
                                                                                                                                 tatgGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagta        <  2:2924724/148‑1 (MQ=255)
                                                                                                                                  attGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagtag       >  2:1121492/1‑148 (MQ=255)
                                                                                                                                   ttGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagtagt      >  1:2053458/1‑148 (MQ=255)
                                                                                                                                   ttGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagtagt      >  2:1033781/1‑148 (MQ=255)
                                                                                                                                    tGCCCGGAGAGTCTGAAGCTTAATCTTTACTGGCAGGGGGCGATTCCGGGATTTTACGCCAGCGCAGGGCAACCACACTGGCGCTTATGTTGCATGGGTATGGCGCAGGCGGGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTc     <  2:3022400/148‑1 (MQ=255)
                                                                                                                                      cccGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAg   >  1:1426218/1‑148 (MQ=255)
                                                                                                                                       ccGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGc  <  2:1831205/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
CGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGC  >  CP009273/4088709‑4088991

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: ‑