Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,272,550 G→T 6.7% G93V (GGG→GTG)  kbaZ → tagatose 6‑phosphate aldolase 1, kbaZ subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,272,5500GT 6.7% 838.4 / 32.2 239G93V (GGG→GTG) kbaZtagatose 6‑phosphate aldolase 1, kbaZ subunit
Reads supporting (aligned to +/- strand):  ref base G (96/127);  new base T (11/5);  total (107/132)
Fisher's exact test for biased strand distribution p-value = 6.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 249 total aligned reads displayed.

GTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGAC  >  CP009273/3272401‑3272698
                                                                                                                                                     |                                                                                                                                                    
gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTggg                                                                                                                                                      <  1:2358048/150‑1 (MQ=255)
gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTggg                                                                                                                                                      <  2:520246/150‑1 (MQ=255)
    cTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCa                                                                                                                                                  >  1:101482/1‑150 (MQ=255)
    cTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCa                                                                                                                                                  >  1:100990/1‑150 (MQ=255)
       aTTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAAc                                                                                                                                               >  2:3841776/1‑150 (MQ=255)
         tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCTCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCg                                                                                                                                             >  2:2695380/1‑150 (MQ=255)
         tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCg                                                                                                                                             <  1:4145731/150‑1 (MQ=255)
          gAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGc                                                                                                                                            >  1:3920043/1‑150 (MQ=255)
          gAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGc                                                                                                                                            >  2:602428/1‑150 (MQ=255)
           aaGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCt                                                                                                                                           >  1:2870340/1‑150 (MQ=255)
              cAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGAGGCGTTGAGTCTGGGTGGTGACCATCTGGGGCCCAACCGCTGGc                                                                                                                                        >  2:2290465/1‑150 (MQ=255)
               aaCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCa                                                                                                                                       >  1:4500893/1‑150 (MQ=255)
                 ccTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCaaa                                                                                                                                     <  2:727326/150‑1 (MQ=255)
                  cTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCaaaa                                                                                                                                    >  1:1816097/1‑150 (MQ=255)
                    ccAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAAcc                                                                                                                                  >  2:2076541/1‑150 (MQ=255)
                     cAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCt                                                                                                                                 <  1:1326788/150‑1 (MQ=255)
                         cAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccg                                                                                                                             <  1:981057/150‑1 (MQ=255)
                         cAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGNGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccg                                                                                                                             <  1:4791844/150‑1 (MQ=255)
                          aGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccgg                                                                                                                            <  2:3443916/150‑1 (MQ=255)
                            gTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGTGCCAAACCGCTGGCAAAACCTgccggcc                                                                                                                          >  1:262757/1‑150 (MQ=255)
                            gTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccggcc                                                                                                                          <  2:1072178/150‑1 (MQ=255)
                             tAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccggccg                                                                                                                         >  2:1496338/1‑150 (MQ=255)
                              aGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGc                                                                                                                        >  1:3501086/1‑150 (MQ=255)
                                 ccAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCa                                                                                                                     >  2:368474/1‑150 (MQ=255)
                                   aGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAgg                                                                                                                   >  2:4237586/1‑150 (MQ=255)
                                    gTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGc                                                                                                                  >  2:376743/1‑150 (MQ=255)
                                       cggcggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAAt                                                                                                               <  1:92341/150‑1 (MQ=255)
                                       cggcggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAAt                                                                                                               <  1:2280993/150‑1 (MQ=255)
                                       cggcggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAAt                                                                                                               <  1:92328/150‑1 (MQ=255)
                                       cggcggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAAt                                                                                                               <  1:3978200/150‑1 (MQ=255)
                                       cggcggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAAt                                                                                                               >  2:4095694/1‑150 (MQ=255)
                                             ttATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCaa                                                                                                         <  1:721056/150‑1 (MQ=255)
                                              tataCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAAt                                                                                                        <  1:3026524/150‑1 (MQ=255)
                                                taCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGc                                                                                                      >  1:604087/1‑150 (MQ=255)
                                                   cGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCga                                                                                                   <  2:2783047/150‑1 (MQ=255)
                                                   cGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCga                                                                                                   <  1:301578/150‑1 (MQ=255)
                                                   cGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCga                                                                                                   <  2:2553666/150‑1 (MQ=255)
                                                    ggAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCgat                                                                                                  >  1:3552791/1‑150 (MQ=255)
                                                    ggAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCgat                                                                                                  >  2:1559393/1‑150 (MQ=255)
                                                      aaTGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCgatga                                                                                                >  2:2375538/1‑150 (MQ=255)
                                                      aaTGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCgatga                                                                                                <  2:1846950/150‑1 (MQ=255)
                                                          aCGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGTGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGAtgatttg                                                                                            >  2:4976572/1‑150 (MQ=255)
                                                           cgcccgccGATTTTCGCGGCTTTGTCTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGAtgatttga                                                                                           <  2:4644435/150‑1 (MQ=255)
                                                            gcccgccGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGTGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGAtgatttgat                                                                                          >  1:3062265/1‑150 (MQ=255)
                                                             cccgccGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGAtgatttgatt                                                                                         <  2:944635/150‑1 (MQ=255)
                                                             cccgccGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGAtgatttgatt                                                                                         >  1:2539771/1‑150 (MQ=255)
                                                               cgccgATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTaa                                                                                       <  1:1874049/150‑1 (MQ=255)
                                                                 ccgATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTaaaa                                                                                     >  2:3915535/1‑150 (MQ=255)
                                                                  cgATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAg                                                                                    <  2:2981188/150‑1 (MQ=255)
                                                                  cgATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAg                                                                                    <  2:3738708/150‑1 (MQ=255)
                                                                   gATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGc                                                                                   <  1:5144940/150‑1 (MQ=255)
                                                                     ttttCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTa                                                                                 <  2:5011119/150‑1 (MQ=255)
                                                                         cgcgGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGtt                                                                             >  2:2011876/1‑150 (MQ=255)
                                                                          gcgGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTg                                                                            <  1:1994467/150‑1 (MQ=255)
                                                                               tttgtttgTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCa                                                                       <  2:1753035/150‑1 (MQ=255)
                                                                                ttgtttgtCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAg                                                                      <  2:2527421/150‑1 (MQ=255)
                                                                                  gtttgtCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGa                                                                    <  1:4924082/150‑1 (MQ=255)
                                                                                    ttgtCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGAtt                                                                  >  1:2174032/1‑150 (MQ=255)
                                                                                         aGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGTGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGAGGATTTGATTAAAAGCTACGTTGCGGCAGGATTCaaaa                                                             >  2:2186888/1‑150 (MQ=255)
                                                                                         aGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCaaaa                                                             <  2:1217011/150‑1 (MQ=255)
                                                                                         aGCTCGCCGACTCGTTGAATGTCCCGCAGGATGCGTTGATTCTGGGGGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCaaaa                                                             <  2:4067277/150‑1 (MQ=255)
                                                                                               ccGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATcc                                                       <  2:3933100/150‑1 (MQ=255)
                                                                                                cGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCa                                                      >  2:1409211/1‑150 (MQ=255)
                                                                                                  aCTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCAcc                                                    >  1:3525357/1‑150 (MQ=255)
                                                                                                    tCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCtt                                                  >  2:2598330/1‑150 (MQ=255)
                                                                                                     cGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTg                                                 <  1:1541447/150‑1 (MQ=255)
                                                                                                      gTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCCAAAAAATCCACCTGGa                                                >  2:334949/1‑150 (MQ=255)
                                                                                                      gTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGa                                                >  2:2451871/1‑150 (MQ=255)
                                                                                                         gAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTg                                             >  2:3431243/1‑150 (MQ=255)
                                                                                                         gAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTg                                             <  1:374199/150‑1 (MQ=255)
                                                                                                          aaTTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCACTGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTgc                                            >  2:5075601/1‑150 (MQ=255)
                                                                                                           aTTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGTGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTgca                                           <  2:1116113/150‑1 (MQ=255)
                                                                                                             ttCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTgcagc                                         <  2:2870319/150‑1 (MQ=255)
                                                                                                                   cAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTcc                                   <  2:71439/150‑1 (MQ=255)
                                                                                                                   cAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTcc                                   >  2:3541439/1‑150 (MQ=255)
                                                                                                                   cAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTcc                                   <  2:908733/150‑1 (MQ=255)
                                                                                                                    aGGATGCGTTGATTCTGGGTGGTGACCATCTGGTGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCt                                  >  1:669005/1‑150 (MQ=255)
                                                                                                                    aGGATGCGTTGATTCTGGGTGGTGACCATCTGGTGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCt                                  >  1:1791518/1‑150 (MQ=255)
                                                                                                                    aGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCt                                  <  2:3501086/150‑1 (MQ=255)
                                                                                                                    aGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCt                                  >  2:4347134/1‑150 (MQ=255)
                                                                                                                     ggATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTg                                 <  1:2413425/150‑1 (MQ=255)
                                                                                                                        tGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCa                              <  2:4710920/150‑1 (MQ=255)
                                                                                                                           gTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGa                           >  2:2887666/1‑150 (MQ=255)
                                                                                                                            ttGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGAc                          <  1:1489897/150‑1 (MQ=255)
                                                                                                                                ttCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATc                      <  1:1655316/150‑1 (MQ=255)
                                                                                                                                ttCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATc                      <  1:2153888/150‑1 (MQ=255)
                                                                                                                                 tCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATcc                     <  1:4696922/150‑1 (MQ=255)
                                                                                                                                  cTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCg                    <  2:1823418/150‑1 (MQ=255)
                                                                                                                                  cTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCg                    <  1:2284775/150‑1 (MQ=255)
                                                                                                                                   tGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGa                   >  2:3069837/1‑150 (MQ=255)
                                                                                                                                     ggtggtGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGAtt                 <  2:4376912/150‑1 (MQ=255)
                                                                                                                                      gtggtgACCATCTGGTGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTc                <  1:3293803/150‑1 (MQ=255)
                                                                                                                                      gtggtgACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTc                <  1:4352726/150‑1 (MQ=255)
                                                                                                                                           gACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTa           <  2:3291857/150‑1 (MQ=255)
                                                                                                                                                tCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACtga      >  1:3221730/1‑150 (MQ=255)
                                                                                                                                                 cTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACtgat     >  2:653955/1‑150 (MQ=255)
                                                                                                                                                 cTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTGAACtgat     >  1:1100415/1‑150 (MQ=255)
                                                                                                                                                   ggggCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCACTGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGAGTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACtgatga   >  2:4255685/1‑150 (MQ=255)
                                                                                                                                                   ggggCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACtgatga   >  1:652192/1‑150 (MQ=255)
                                                                                                                                                   ggggCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACtgatga   >  1:4555755/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGAC  >  CP009273/3272401‑3272698

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