Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,334,687 G→T 100% A21A (GCG→GCT kdsD → D‑arabinose 5‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,334,6870GT97.5% 681.3 / 6.4 215A21A (GCG→GCTkdsDD‑arabinose 5‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base T (91/118);  total (95/120)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 215 total aligned reads displayed.

ATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCGATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAATGGCG  >  CP009273/3334542‑3334828
                                                                                                                                                 |                                                                                                                                             
attattAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTAt                                                                                                                                             >  1:2186562/1‑148 (MQ=255)
 ttattaACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTAtt                                                                                                                                            >  2:2999160/1‑148 (MQ=255)
   attaACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGa                                                                                                                                          <  2:3769206/148‑1 (MQ=255)
    ttaACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGaa                                                                                                                                         <  2:2282256/148‑1 (MQ=255)
     taACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAAc                                                                                                                                        <  1:2560912/148‑1 (MQ=255)
       aCTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAATAAGTCCTGTCGATTGAACGt                                                                                                                                      >  2:853473/1‑148 (MQ=255)
         tGGCGGGTTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGa                                                                                                                                    <  2:2948308/148‑1 (MQ=255)
         tGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGa                                                                                                                                    <  1:295825/148‑1 (MQ=255)
         tGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGa                                                                                                                                    <  1:499146/148‑1 (MQ=255)
          ggCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGaa                                                                                                                                   <  1:1759819/148‑1 (MQ=255)
          ggCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGaa                                                                                                                                   <  1:343775/148‑1 (MQ=255)
           gCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAAt                                                                                                                                  >  2:2809527/1‑148 (MQ=255)
           gCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAAt                                                                                                                                  >  1:212225/1‑148 (MQ=255)
              gATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGcc                                                                                                                               <  2:1596324/148‑1 (MQ=255)
               aTTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCt                                                                                                                              >  1:1006052/1‑148 (MQ=255)
               aTTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCt                                                                                                                              >  1:1816183/1‑148 (MQ=255)
               aTTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCt                                                                                                                              <  2:1473881/148‑1 (MQ=255)
               aTTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCt                                                                                                                              >  2:504907/1‑148 (MQ=255)
                 ttATCGTATGGCTGGCGATGTTGTACTGTTTATCGCCAATACTCGTTGAATAAATGGAAACTCATTATGTCGCACGTAGAGTTACAACCTGGTTTTGAATTTCATCAAGCAGTTAAAGAAGTCCTTTCTATTTAAAGTGAATTCCTgt                                                                                                                            >  1:512111/1‑147 (MQ=255)
                      gTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAg                                                                                                                       <  2:2807210/148‑1 (MQ=255)
                       tATGGCTGGCGATGTTGTACTGGTTATCGCCAATCCTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGc                                                                                                                      <  2:2032645/148‑1 (MQ=255)
                       tATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGc                                                                                                                      <  1:23125/148‑1 (MQ=255)
                       tATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGc                                                                                                                      <  2:2103578/148‑1 (MQ=255)
                       tATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGc                                                                                                                      <  1:2545075/148‑1 (MQ=255)
                       tATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGc                                                                                                                      <  1:1790439/148‑1 (MQ=255)
                       tATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGc                                                                                                                      <  1:1105138/148‑1 (MQ=255)
                        aTGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCt                                                                                                                     <  2:168585/148‑1 (MQ=255)
                          ggctggcGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTg                                                                                                                   <  1:3633/148‑1 (MQ=255)
                           gctggcGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGa                                                                                                                  <  1:3700543/148‑1 (MQ=255)
                           gctggcGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGa                                                                                                                  >  1:1707960/1‑148 (MQ=255)
                            ctggcTATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGACACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGACGTCCTGGCTATTGACCGTGAATGCCTGGCGGAGCTTGAt                                                                                                                 <  1:1238703/148‑1 (MQ=255)
                              ggcGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCa                                                                                                               <  1:970452/148‑1 (MQ=255)
                              ggcGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCa                                                                                                               >  1:893953/1‑148 (MQ=255)
                              ggcGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCa                                                                                                               <  2:777136/148‑1 (MQ=255)
                              ggcGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCa                                                                                                               <  1:101276/148‑1 (MQ=255)
                                 gATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATa                                                                                                            <  1:1966652/148‑1 (MQ=255)
                                  aTGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATAc                                                                                                           <  2:1062463/148‑1 (MQ=255)
                                   tgttgtACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACa                                                                                                          >  2:2492231/1‑148 (MQ=255)
                                   tgttgtACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACa                                                                                                          >  2:2821642/1‑148 (MQ=255)
                                     ttgtACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACatc                                                                                                        >  1:864757/1‑148 (MQ=255)
                                                aTCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAAttt                                                                                             <  1:436671/148‑1 (MQ=255)
                                                 tCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTc                                                                                            <  1:3193125/148‑1 (MQ=255)
                                                  cGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCa                                                                                           <  2:3111277/148‑1 (MQ=255)
                                                    ccAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACg                                                                                         >  2:785630/1‑148 (MQ=255)
                                                     cAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGc                                                                                        <  1:1722235/148‑1 (MQ=255)
                                                     cAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGc                                                                                        >  2:151969/1‑148 (MQ=255)
                                                      aaTACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCt                                                                                       <  1:270742/148‑1 (MQ=255)
                                                       aTACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCtt                                                                                      >  2:356078/1‑148 (MQ=255)
                                                         aCTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGc                                                                                    >  2:3235376/1‑148 (MQ=255)
                                                         aCTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGc                                                                                    >  2:3215245/1‑148 (MQ=255)
                                                          cTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGcc                                                                                   >  2:115404/1‑148 (MQ=255)
                                                          cTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGcc                                                                                   >  1:1105044/1‑148 (MQ=255)
                                                          cTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCGATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGcc                                                                                   <  1:1921453/148‑1 (MQ=255)
                                                                      cTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATg                                                                       >  2:3719521/1‑148 (MQ=255)
                                                                      cTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATg                                                                       >  2:614682/1‑148 (MQ=255)
                                                                       tGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGt                                                                      <  2:2031046/148‑1 (MQ=255)
                                                                         gAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTc                                                                    >  2:1653612/1‑148 (MQ=255)
                                                                            aCGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCGATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTgg                                                                 >  2:751359/1‑148 (MQ=255)
                                                                                aTTATGTCGCACGTAGAGTTACACCCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTa                                                             >  1:2788902/1‑148 (MQ=255)
                                                                                 ttATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTaa                                                            >  1:1132704/1‑148 (MQ=255)
                                                                                  tATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTaaa                                                           <  1:3815431/148‑1 (MQ=255)
                                                                                  tATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGGTCTGGTGTaaa                                                           >  2:1799432/1‑148 (MQ=255)
                                                                                   aTGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAg                                                          <  2:3544595/148‑1 (MQ=255)
                                                                                   aTGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAg                                                          <  2:848512/148‑1 (MQ=255)
                                                                                     gTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAATggg                                                        >  2:2004091/1‑148 (MQ=255)
                                                                                     gTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAggg                                                        <  1:3592048/148‑1 (MQ=255)
                                                                                              agagTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTgtc                                               >  2:947134/1‑148 (MQ=255)
                                                                                                        ccGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGa                                     >  2:1233832/1‑148 (MQ=255)
                                                                                                        ccGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGa                                     >  2:2611455/1‑148 (MQ=255)
                                                                                                          gggTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATg                                   <  1:2091504/148‑1 (MQ=255)
                                                                                                          gggTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATg                                   <  1:2856748/148‑1 (MQ=255)
                                                                                                           ggTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATgg                                  >  2:3358485/1‑148 (MQ=255)
                                                                                                              tttGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGaa                               >  1:3043871/1‑148 (MQ=255)
                                                                                                               ttGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGaaa                              <  1:1331685/148‑1 (MQ=255)
                                                                                                               ttGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGaaa                              <  1:2821642/148‑1 (MQ=255)
                                                                                                                tGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGaaaa                             <  2:183198/148‑1 (MQ=255)
                                                                                                                tGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGaaaa                             <  1:2147511/148‑1 (MQ=255)
                                                                                                                   cTTTCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCg                          <  1:2999160/148‑1 (MQ=255)
                                                                                                                    tttCAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCgg                         <  1:3107549/148‑1 (MQ=255)
                                                                                                                       cAGCAAGCAGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGc                      <  1:3729975/148‑1 (MQ=255)
                                                                                                                        agcaagcaGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCa                     <  2:3789288/148‑1 (MQ=255)
                                                                                                                         gcaagcaGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCat                    >  2:2746586/1‑148 (MQ=255)
                                                                                                                            agcaGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTACAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATAtt                 >  2:1108678/1‑148 (MQ=255)
                                                                                                                            agcaGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATAtt                 >  1:3579053/1‑148 (MQ=255)
                                                                                                                            agcaGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATAtt                 <  2:1377442/148‑1 (MQ=255)
                                                                                                                              caGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTgg               >  2:1012339/1‑148 (MQ=255)
                                                                                                                              caGGTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTgg               >  2:3637783/1‑148 (MQ=255)
                                                                                                                                 gTAAAGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCg            >  2:18716/1‑148 (MQ=255)
                                                                                                                                   aaaGAAGTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGaa          >  1:2346961/1‑148 (MQ=255)
                                                                                                                                     agaagTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGaaaa        <  2:1195998/148‑1 (MQ=255)
                                                                                                                                     agaagTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGaaaa        <  2:1196295/148‑1 (MQ=255)
                                                                                                                                     agaagTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGaaaa        >  2:436501/1‑148 (MQ=255)
                                                                                                                                     agaagTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGaaaa        >  2:100258/1‑148 (MQ=255)
                                                                                                                                     agaagTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGaaaa        >  2:1976196/1‑148 (MQ=255)
                                                                                                                                      gaagTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGaaaaa       <  1:151969/148‑1 (MQ=255)
                                                                                                                                       aagTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAAt      <  1:2161804/148‑1 (MQ=255)
                                                                                                                                        agTCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAATg     >  1:1367341/1‑148 (MQ=255)
                                                                                                                                          tCCTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAATggc   >  1:2399509/1‑148 (MQ=255)
                                                                                                                                           ccTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAATggcg  <  1:616021/148‑1 (MQ=255)
                                                                                                                                           ccTGGCTATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAATggcg  <  1:614532/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
ATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATAACTGGAAACGCATTATGTCGCACGTAGAGTTACAACCGGGTTTTGACTTTCAGCAAGCAGGTAAAGAAGTCCTGGCGATTGAACGTGAATGCCTGGCGGAGCTTGATCAATACATCAATCAGAATTTCACGCTTGCCTGTGAAAAGATGTTCTGGTGTAAAGGGAAAGTTGTCGTCATGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAATGGCG  >  CP009273/3334542‑3334828

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: