Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 119,536 G→T 100% P11P (CCG→CCT aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273119,5360GT97.9% 925.8 / 2.5 300P11P (CCG→CCTaceEpyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (3/7);  new base T (172/117);  total (176/124)
Fisher's exact test for biased strand distribution p-value = 9.90e-02
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 300 total aligned reads displayed.

GCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCAGCAACTA  >  CP009273/119389‑119682
                                                                                                                                                   |                                                                                                                                                  
  atatAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCg                                                                                                                                                >  1:4569127/1‑150 (MQ=255)
    ataAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGaa                                                                                                                                              >  2:4898772/1‑150 (MQ=255)
    ataAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGCTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGaa                                                                                                                                              >  2:164262/1‑150 (MQ=255)
     taAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGaaa                                                                                                                                             >  2:4485802/1‑150 (MQ=255)
         aGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTcg                                                                                                                                         >  1:4736512/1‑150 (MQ=255)
         aGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTcg                                                                                                                                         <  1:1805672/150‑1 (MQ=255)
           cgcgGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTcgcg                                                                                                                                       <  1:792922/150‑1 (MQ=255)
           cgcgGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTcgcg                                                                                                                                       <  1:2113310/150‑1 (MQ=255)
           cgcgGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTcgcg                                                                                                                                       <  1:3171612/150‑1 (MQ=255)
             cgGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGAc                                                                                                                                     >  2:2006986/1‑150 (MQ=255)
               gCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTg                                                                                                                                   >  2:1076992/1‑150 (MQ=255)
                 aaCTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGc                                                                                                                                 >  2:641081/1‑150 (MQ=255)
                     aaaCGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCa                                                                                                                             <  1:3318507/150‑1 (MQ=255)
                     aaaCGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGAAAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCa                                                                                                                             >  1:1068014/1‑150 (MQ=255)
                       aCGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAgg                                                                                                                           <  2:4970321/150‑1 (MQ=255)
                        cGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGc                                                                                                                          <  2:1397691/150‑1 (MQ=255)
                        cGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGc                                                                                                                          <  2:5066586/150‑1 (MQ=255)
                             aaCCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCg                                                                                                                     <  1:4200836/150‑1 (MQ=255)
                              aCCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGGCAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGa                                                                                                                    >  1:2316668/1‑150 (MQ=255)
                               ccTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGac                                                                                                                   >  2:4202699/1‑149 (MQ=255)
                                cTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAAt                                                                                                                  >  1:891788/1‑150 (MQ=255)
                                cTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAAt                                                                                                                  >  1:2564400/1‑150 (MQ=255)
                                 tGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATc                                                                                                                 <  1:2840421/150‑1 (MQ=255)
                                  gTCTGATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATca                                                                                                                <  2:4817891/150‑2 (MQ=255)
                                   tCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCgg                                                                                                               >  2:4435960/1‑150 (MQ=255)
                                       aTTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCAt                                                                                                           >  1:2918563/1‑150 (MQ=255)
                                         tGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATcc                                                                                                         <  1:1477983/150‑1 (MQ=255)
                                          gAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCg                                                                                                        <  1:3566651/150‑1 (MQ=255)
                                           aGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGt                                                                                                       <  2:1800575/150‑1 (MQ=255)
                                           aGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGt                                                                                                       >  1:4295573/1‑150 (MQ=255)
                                           aGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGt                                                                                                       >  2:727025/1‑150 (MQ=255)
                                               ttCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTgaag                                                                                                   >  2:62896/1‑150 (MQ=255)
                                               ttCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTgaag                                                                                                   >  1:3784907/1‑150 (MQ=255)
                                                tCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTgaata                                                                                                  >  1:5045526/1‑148 (MQ=255)
                                                    gCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGgt                                                                                              <  1:741512/150‑1 (MQ=255)
                                                     cGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACGGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGgtg                                                                                             <  1:2336333/150‑1 (MQ=255)
                                                         agTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGa                                                                                         >  2:1154616/1‑150 (MQ=255)
                                                          gTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAg                                                                                        >  1:911339/1‑150 (MQ=255)
                                                             cAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGt                                                                                     >  1:243120/1‑150 (MQ=255)
                                                                 gggACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTc                                                                                 >  2:4572564/1‑150 (MQ=255)
                                                                 gggACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTc                                                                                 >  2:3907255/1‑150 (MQ=255)
                                                                      tggTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTAt                                                                            >  1:1832538/2‑150 (MQ=255)
                                                                      aGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTAt                                                                            >  1:2051996/1‑150 (MQ=255)
                                                                       ggTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATc                                                                           >  1:3087939/1‑150 (MQ=255)
                                                                        gTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCt                                                                          >  2:3720302/1‑150 (MQ=255)
                                                                        gTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCt                                                                          >  1:192623/1‑150 (MQ=255)
                                                                           ccAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGAt                                                                       >  1:4965442/1‑150 (MQ=255)
                                                                            cAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATc                                                                      >  2:1250567/1‑150 (MQ=255)
                                                                              gAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGa                                                                    >  2:875189/1‑150 (MQ=255)
                                                                               aaaaCTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGAc                                                                   >  1:40421/1‑150 (MQ=255)
                                                                                 aaCTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCa                                                                 <  2:3427642/150‑1 (MQ=255)
                                                                                 aaCTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCa                                                                 <  1:3957588/150‑1 (MQ=255)
                                                                                     cAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACtg                                                             <  1:4685985/150‑1 (MQ=255)
                                                                                     cAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACtg                                                             >  1:3404643/1‑150 (MQ=255)
                                                                                     cAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACtg                                                             <  2:569483/150‑1 (MQ=255)
                                                                                      aaCGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACtgc                                                            <  2:76621/150‑1 (MQ=255)
                                                                                        cGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACtgctt                                                          >  2:4069343/1‑150 (MQ=255)
                                                                                        cGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACtgctt                                                          <  2:911286/150‑1 (MQ=255)
                                                                                           tattaGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACtgcttgct                                                       <  2:4622004/150‑1 (MQ=255)
                                                                                           tattaGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACtgcttgct                                                       <  1:4648006/150‑1 (MQ=255)
                                                                                              tagatagaTAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGaa                                                    <  2:3261149/150‑1 (MQ=255)
                                                                                               agatagatAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAg                                                   >  2:3768802/1‑150 (MQ=255)
                                                                                               agatagatAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAg                                                   >  1:1117653/1‑150 (MQ=255)
                                                                                               agatagatAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAg                                                   <  1:2720430/150‑1 (MQ=255)
                                                                                                      taAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCaa                                            <  1:4478097/150‑1 (MQ=255)
                                                                                                        aGGAATACCCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAg                                          <  2:2977097/150‑1 (MQ=255)
                                                                                                        aGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAg                                          <  2:4369854/150‑1 (MQ=255)
                                                                                                        aGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAg                                          >  1:588516/1‑150 (MQ=255)
                                                                                                         ggAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAgg                                         >  2:2556640/1‑150 (MQ=255)
                                                                                                         ggAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAgg                                         >  2:3328386/1‑150 (MQ=255)
                                                                                                          gAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGc                                        <  1:817163/150‑1 (MQ=255)
                                                                                                            aTAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCgg                                      >  1:966987/1‑150 (MQ=255)
                                                                                                             tAACCCATGTCAGAACGTTTCCCAACTGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGgt                                     <  2:1879291/150‑1 (MQ=255)
                                                                                                             tAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTACTGAAGCCCGCAAAGGCGgt                                     <  1:39295/150‑1 (MQ=255)
                                                                                                               aCCCATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGgtgt                                   <  1:4842313/150‑1 (MQ=255)
                                                                                                                cccATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTa                                  >  1:3521478/1‑150 (MQ=255)
                                                                                                                 ccATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTaa                                 <  1:4571002/150‑1 (MQ=255)
                                                                                                                  cATGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTaaa                                >  2:4553257/1‑150 (MQ=255)
                                                                                                                   aTGTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAAc                               <  2:93587/150‑1 (MQ=255)
                                                                                                                     gTCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGt                             >  2:1322804/1‑150 (MQ=255)
                                                                                                                      tCAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTa                            <  1:1239275/150‑1 (MQ=255)
                                                                                                                       cAGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAg                           <  1:162565/150‑1 (MQ=255)
                                                                                                                        aGAACGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGc                          >  1:3727577/1‑150 (MQ=255)
                                                                                                                          aaCGTTTCCCAAATGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGCTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCg                        >  1:1329863/1‑150 (MQ=255)
                                                                                                                                  ccAAATGACGTGGATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGcaca                <  1:219126/150‑1 (MQ=255)
                                                                                                                                    aaaaGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAgg              >  1:1735061/1‑150 (MQ=255)
                                                                                                                                    aaaTGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAgg              >  2:2471340/1‑150 (MQ=255)
                                                                                                                                     aaTGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGt             <  2:2816682/150‑1 (MQ=255)
                                                                                                                                     aaTGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGt             >  2:4866147/1‑150 (MQ=255)
                                                                                                                                     aaTGACGTGGATCCCATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGt             >  2:1078858/1‑150 (MQ=255)
                                                                                                                                      aTGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTCAACGTAGCCGCAGGCACAGGTa            >  2:2973172/1‑150 (MQ=255)
                                                                                                                                      aTGACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTa            >  1:4863712/1‑150 (MQ=255)
                                                                                                                                        gccGTGGATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATc          <  1:4361230/148‑1 (MQ=255)
                                                                                                                                        gACGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATc          <  1:3435984/150‑1 (MQ=255)
                                                                                                                                         aCGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCa         <  1:4082182/150‑1 (MQ=255)
                                                                                                                                         aCGTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCa         <  2:1611772/150‑1 (MQ=255)
                                                                                                                                           gTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCAGc       <  1:4051772/150‑1 (MQ=255)
                                                                                                                                           gTGGATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCAGc       >  1:1556770/1‑150 (MQ=255)
                                                                                                                                             ggATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCAGCaa     <  2:2316668/150‑1 (MQ=255)
                                                                                                                                              gATCCTATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCAGCAAc    >  1:3229001/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCAGCAACTA  >  CP009273/119389‑119682

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: