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.8% 872.6 / 7.0 271R93S (CGT→AGT) rhaRtranscriptional activator of rhaSR
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (140/125);  total (143/128)
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 271 total aligned reads displayed.

CTGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCATGT  >  CP009273/4088705‑4088995
                                                                                                                                                 |                                                                                                                                                 
cTGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCt                                                                                                                                               >  1:3477577/1‑150 (MQ=255)
 tGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTg                                                                                                                                              <  2:2049595/150‑1 (MQ=255)
 tGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTg                                                                                                                                              >  1:3387232/1‑150 (MQ=255)
 tGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTg                                                                                                                                              >  2:4643067/1‑150 (MQ=255)
  ggCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGa                                                                                                                                             <  1:192673/150‑1 (MQ=255)
   gcgcGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCATACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGaa                                                                                                                                            >  1:1001587/1‑150 (MQ=255)
     gcgGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGc                                                                                                                                          >  1:1638189/1‑150 (MQ=255)
       ggTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTg                                                                                                                                        >  2:4731309/1‑150 (MQ=255)
       ggTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTg                                                                                                                                        >  2:2163700/1‑150 (MQ=255)
       ggTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTg                                                                                                                                        <  2:2341791/150‑1 (MQ=255)
          aaTGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAAt                                                                                                                                     <  2:3647818/150‑1 (MQ=255)
          aaTGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAAt                                                                                                                                     <  1:3177787/150‑1 (MQ=255)
              gCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTg                                                                                                                                 <  1:1980488/150‑1 (MQ=255)
                cTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGAc                                                                                                                               >  1:555905/1‑150 (MQ=255)
                 tGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACt                                                                                                                              <  1:673701/150‑1 (MQ=255)
                  gCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTg                                                                                                                             <  2:103136/150‑1 (MQ=255)
                   cATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTgg                                                                                                                            >  2:6644/1‑150 (MQ=255)
                       tACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAgg                                                                                                                        <  2:4110618/150‑1 (MQ=255)
                       tACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAgg                                                                                                                        >  2:535844/1‑150 (MQ=255)
                       tACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAgg                                                                                                                        >  2:2244650/1‑150 (MQ=255)
                        aCTCACCGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACTCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAggg                                                                                                                       <  2:3096200/150‑1 (MQ=255)
                        aCTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAggg                                                                                                                       >  1:1328707/1‑150 (MQ=255)
                            aaCGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCg                                                                                                                   <  1:4888392/150‑1 (MQ=255)
                              cGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGAt                                                                                                                 <  1:963360/150‑1 (MQ=255)
                               gATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGAtt                                                                                                                >  2:2409986/1‑150 (MQ=255)
                                 tCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTcc                                                                                                              <  2:1557708/150‑1 (MQ=255)
                                 tCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTcc                                                                                                              >  2:3157820/1‑150 (MQ=255)
                                 tCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTcc                                                                                                              <  1:4624282/150‑1 (MQ=255)
                                 tCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTcc                                                                                                              <  1:4372631/150‑1 (MQ=255)
                                   gCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCgg                                                                                                            >  1:3721415/1‑150 (MQ=255)
                                         aTCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTa                                                                                                      >  1:3135113/1‑150 (MQ=255)
                                           cGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAAc                                                                                                    >  2:1774677/1‑150 (MQ=255)
                                             cATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGc                                                                                                  >  1:990668/1‑150 (MQ=255)
                                              aTTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGcc                                                                                                 <  1:2073494/150‑1 (MQ=255)
                                                tACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAg                                                                                               <  1:3281903/150‑1 (MQ=255)
                                                tACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAg                                                                                               <  1:697580/150‑1 (MQ=255)
                                                  cccGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAgcg                                                                                             >  2:202476/1‑150 (MQ=255)
                                                   ccGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAgcgc                                                                                            >  2:1100651/1‑150 (MQ=255)
                                                    cGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCa                                                                                           <  1:817192/150‑1 (MQ=255)
                                                     gTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAg                                                                                          <  1:529248/150‑1 (MQ=255)
                                                      tGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAgg                                                                                         <  2:4462985/150‑1 (MQ=255)
                                                       ggCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAggg                                                                                        >  1:4556832/1‑150 (MQ=255)
                                                       ggCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAggg                                                                                        >  1:3630867/1‑150 (MQ=255)
                                                       ggCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAggg                                                                                        >  1:2206490/1‑150 (MQ=255)
                                                         cGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCa                                                                                      >  2:2856882/1‑150 (MQ=255)
                                                            tctcTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAAcc                                                                                   <  2:4226798/150‑1 (MQ=255)
                                                            tctcTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAAcc                                                                                   >  1:2014018/1‑150 (MQ=255)
                                                             ctctTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACca                                                                                  >  1:181924/1‑150 (MQ=255)
                                                                 tttACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACt                                                                              <  1:3948527/150‑1 (MQ=255)
                                                                   tACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTgg                                                                            <  1:1078519/150‑1 (MQ=255)
                                                                   tACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTgg                                                                            <  1:608135/150‑1 (MQ=255)
                                                                     cATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgcg                                                                          <  2:3469010/150‑1 (MQ=255)
                                                                      aTTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgcgc                                                                         >  2:2695287/1‑150 (MQ=255)
                                                                       ttCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCt                                                                        <  2:9162/150‑1 (MQ=255)
                                                                          aTGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAg                                                                     <  2:3404539/150‑1 (MQ=255)
                                                                          aTGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAg                                                                     >  2:3584966/1‑150 (MQ=255)
                                                                             cTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTa                                                                  >  1:5021434/1‑150 (MQ=255)
                                                                               gACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGc                                                                >  2:4057611/1‑150 (MQ=255)
                                                                                aCGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCa                                                               <  2:2362997/150‑1 (MQ=255)
                                                                                 cGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTCGGTAGCAt                                                              >  1:3886659/1‑150 (MQ=255)
                                                                                 cGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCAt                                                              >  2:4590360/1‑150 (MQ=255)
                                                                                  gATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCCTg                                                             >  2:3809210/1‑150 (MQ=255)
                                                                                     aaaCACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCCTgggg                                                          >  2:3282467/1‑150 (MQ=255)
                                                                                      aaCACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGa                                                         >  2:1347883/1‑150 (MQ=255)
                                                                                      aaCACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGa                                                         >  2:3184725/1‑150 (MQ=255)
                                                                                       acacTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGAt                                                        <  2:785395/150‑1 (MQ=255)
                                                                                         acTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTCGCATGGGGATgg                                                      >  2:2926776/1‑150 (MQ=255)
                                                                                         acTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATgg                                                      >  1:2628784/1‑150 (MQ=255)
                                                                                         acTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATgg                                                      >  2:1457678/1‑150 (MQ=255)
                                                                                             cTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCa                                                  <  2:870420/150‑1 (MQ=255)
                                                                                              tACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAg                                                 <  2:990542/150‑1 (MQ=255)
                                                                                                cGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGgc                                               >  2:3396687/1‑150 (MQ=255)
                                                                                                cGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGgc                                               <  1:3226938/150‑1 (MQ=255)
                                                                                                 gCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGgcg                                              <  1:1515656/150‑1 (MQ=255)
                                                                                                     ccGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCCTGGGGATGGCGCTGGCGCGGc                                          >  2:4332569/1‑150 (MQ=255)
                                                                                                      cGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCa                                         >  2:4300159/1‑150 (MQ=255)
                                                                                                         tAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGt                                      >  1:4052248/1‑150 (MQ=255)
                                                                                                          aaCGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGtt                                     >  1:4212310/1‑150 (MQ=255)
                                                                                                          aaCGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGtt                                     <  1:2335764/150‑1 (MQ=255)
                                                                                                          aaCGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGAGTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGtt                                     <  1:5083330/150‑1 (MQ=255)
                                                                                                            cGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTAt                                   <  2:4967362/150‑1 (MQ=255)
                                                                                                             gATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATc                                  <  1:2774923/150‑1 (MQ=255)
                                                                                                               tCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCgg                                >  1:4204004/1‑150 (MQ=255)
                                                                                                                  ggTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCa                             >  1:4097969/1‑150 (MQ=255)
                                                                                                                    ttttGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGc                           <  1:4142021/150‑1 (MQ=255)
                                                                                                                        gCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGa                       <  2:1237129/150‑1 (MQ=255)
                                                                                                                             atatTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATg                  <  2:616546/150‑1 (MQ=255)
                                                                                                                                ttatttatTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGaaa               >  1:4267492/1‑150 (MQ=255)
                                                                                                                                ttatttatTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGaaa               <  2:3919002/150‑1 (MQ=255)
                                                                                                                                 tatttagtGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAag              <  2:1697187/150‑1 (MQ=255)
                                                                                                                                  atttattGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagt             <  2:75071/150‑1 (MQ=255)
                                                                                                                                     tattGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagtagt          >  1:4785297/1‑150 (MQ=255)
                                                                                                                                      attGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTc         >  1:1210010/1‑150 (MQ=255)
                                                                                                                                       ttGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTCTCGGTCAGCTTGAGCATGAAAGTAGTCa        >  2:2829724/1‑150 (MQ=255)
                                                                                                                                       ttGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCa        <  2:1950130/150‑1 (MQ=255)
                                                                                                                                        tGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAg       >  1:4756145/1‑150 (MQ=255)
                                                                                                                                        tGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTGATCGGTCAGCTTGAGCATGAAAGTAGTCAg       >  2:4462052/1‑150 (MQ=255)
                                                                                                                                          cccGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCa     >  1:3005863/1‑150 (MQ=255)
                                                                                                                                           ccGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCAt    >  1:1743194/1‑150 (MQ=255)
                                                                                                                                             ggAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCAtgt  >  2:421322/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCATGT  >  CP009273/4088705‑4088995

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: