Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 207,624 G→T 5.2% D487Y (GAT→TAT)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273207,6240GT 5.2% 754.5 / 14.7 192D487Y (GAT→TAT) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (91/91);  new base T (5/5);  total (96/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 203 total aligned reads displayed.

AAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTC  >  CP009273/207475‑207773
                                                                                                                                                     |                                                                                                                                                     
aaGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCg                                                                                                                                                       >  1:4046946/1‑150 (MQ=255)
aaGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCg                                                                                                                                                       >  2:209687/1‑150 (MQ=255)
  gagGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGAt                                                                                                                                                     >  1:3534151/1‑150 (MQ=255)
    ggTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATcc                                                                                                                                                   <  2:1558051/150‑1 (MQ=255)
     gTCCAGCGGCTGCGGGAAGAGTCTGACTGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCTCGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCg                                                                                                                                                  <  2:1431959/150‑1 (MQ=255)
     gTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCg                                                                                                                                                  <  2:101110/150‑1 (MQ=255)
      tCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCgg                                                                                                                                                 <  2:3479960/150‑1 (MQ=255)
       ccAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGt                                                                                                                                                <  1:3017174/150‑1 (MQ=255)
        cAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGtt                                                                                                                                               >  2:1584571/1‑150 (MQ=255)
        cAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGtt                                                                                                                                               <  2:1916534/150‑1 (MQ=255)
         aGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTa                                                                                                                                              >  1:843727/1‑150 (MQ=255)
         aGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTGTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTa                                                                                                                                              >  1:2887129/1‑150 (MQ=255)
            ggCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAg                                                                                                                                           >  2:849916/1‑150 (MQ=255)
             gCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCTATCCGGTTAAAGt                                                                                                                                          <  1:859695/150‑1 (MQ=255)
             gCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGt                                                                                                                                          >  2:3559182/1‑150 (MQ=255)
             gCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGAGATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGt                                                                                                                                          <  2:4085231/150‑1 (MQ=255)
               tGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCa                                                                                                                                        <  2:127770/150‑1 (MQ=255)
                 cGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACt                                                                                                                                      <  1:2552173/150‑1 (MQ=255)
                   ggAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCTATCCGGTTAAAGTCACTAt                                                                                                                                    <  1:2062261/150‑1 (MQ=255)
                     aaGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCCCTAttt                                                                                                                                  >  2:1638385/1‑150 (MQ=255)
                        agTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGa                                                                                                                               >  2:1977528/1‑150 (MQ=255)
                        agTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGa                                                                                                                               >  1:3722330/1‑150 (MQ=255)
                        agTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGa                                                                                                                               >  1:3570664/1‑150 (MQ=255)
                        agTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGa                                                                                                                               <  1:3095480/150‑1 (MQ=255)
                          tCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGaca                                                                                                                             >  1:3070959/1‑150 (MQ=255)
                            tGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACAcc                                                                                                                           >  1:1542055/1‑150 (MQ=255)
                            tGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACAcc                                                                                                                           <  1:268368/150‑1 (MQ=255)
                            tGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACAcc                                                                                                                           <  2:2578235/150‑1 (MQ=255)
                             gACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCg                                                                                                                          >  2:3515914/1‑150 (MQ=255)
                                 gttggttTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGa                                                                                                                      >  2:2349684/1‑150 (MQ=255)
                                          tCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGc                                                                                                             >  2:1116/1‑150 (MQ=255)
                                           cGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCa                                                                                                            <  1:800660/150‑1 (MQ=255)
                                             atatCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAgg                                                                                                          <  2:1052842/150‑1 (MQ=255)
                                             atatCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAgg                                                                                                          <  2:3329446/150‑1 (MQ=255)
                                             atatCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAgg                                                                                                          >  1:1311880/1‑150 (MQ=255)
                                               atCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGc                                                                                                        <  1:118678/150‑1 (MQ=255)
                                                    gCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAatat                                                                                                   <  1:1852210/150‑1 (MQ=255)
                                                     cAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATg                                                                                                  <  2:877080/150‑1 (MQ=255)
                                                     cAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATg                                                                                                  <  1:2773094/150‑1 (MQ=255)
                                                       aCCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAg                                                                                                <  1:1957308/150‑1 (MQ=255)
                                                        ccgccgCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGc                                                                                               <  2:3421953/150‑1 (MQ=255)
                                                        ccgccgCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGc                                                                                               >  1:898731/1‑150 (MQ=255)
                                                        ccgccgCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGc                                                                                               >  1:898697/1‑150 (MQ=255)
                                                                 gTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAgggg                                                                                      >  2:1216104/1‑150 (MQ=255)
                                                                    gATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATc                                                                                   <  2:3828748/150‑1 (MQ=255)
                                                                     aTGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATcc                                                                                  >  1:2508463/1‑150 (MQ=255)
                                                                      tGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATccc                                                                                 <  1:2743792/150‑1 (MQ=255)
                                                                       gAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCcg                                                                                <  1:3880800/150‑1 (MQ=255)
                                                                        aaGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCcgg                                                                               >  1:2362847/1‑150 (MQ=255)
                                                                                 gCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTgg                                                                      <  1:2107437/150‑1 (MQ=255)
                                                                                              gcCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCt                                                         <  2:426324/150‑1 (MQ=255)
                                                                                                cTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcg                                                       >  1:1054339/1‑150 (MQ=255)
                                                                                                cTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcg                                                       >  2:705843/1‑150 (MQ=255)
                                                                                                   gCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcgacg                                                    <  1:1329887/150‑1 (MQ=255)
                                                                                                    cGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcgacga                                                   <  2:253539/150‑1 (MQ=255)
                                                                                                      aaCAGTGGCACGGCTTTATCGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACTCCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAAc                                                 <  1:1778944/150‑1 (MQ=255)
                                                                                                        cAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGt                                               >  1:333170/1‑150 (MQ=255)
                                                                                                         aGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTg                                              <  2:326050/150‑1 (MQ=255)
                                                                                                           tGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTggg                                            <  1:762734/150‑1 (MQ=255)
                                                                                                              cACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATc                                         >  2:659161/1‑150 (MQ=255)
                                                                                                               aCGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCg                                        <  2:3670366/150‑1 (MQ=255)
                                                                                                                    tttAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGAGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtagta                                   <  1:1179363/150‑1 (MQ=255)
                                                                                                                     ttAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGtagtag                                  >  2:2933587/1‑150 (MQ=255)
                                                                                                                      tAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTaga                                 >  1:1702630/1‑150 (MQ=255)
                                                                                                                      tAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTaga                                 >  1:117081/1‑150 (MQ=255)
                                                                                                                       aaCGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTagag                                >  2:2080604/1‑150 (MQ=255)
                                                                                                                        aCGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAgaga                               <  1:684839/150‑1 (MQ=255)
                                                                                                                         cGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaa                              <  1:469626/150‑1 (MQ=255)
                                                                                                                         cGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaa                              >  1:1735473/1‑150 (MQ=255)
                                                                                                                         cGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaa                              >  2:2278689/1‑150 (MQ=255)
                                                                                                                           aTGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaaa                            <  1:1275887/150‑1 (MQ=255)
                                                                                                                               ggATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCgg                        >  2:464639/1‑150 (MQ=255)
                                                                                                                               ggATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCgg                        >  1:3732624/1‑150 (MQ=255)
                                                                                                                                 aTGCCGATCATATTTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACAGGcc                      >  1:1788988/1‑150 (MQ=255)
                                                                                                                                 aTGCCGATCATATGTTTCTCTATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGcc                      <  2:1239044/150‑1 (MQ=255)
                                                                                                                                 aTGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGcc                      >  2:767824/1‑150 (MQ=255)
                                                                                                                                  tGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGccc                     <  2:3216899/150‑1 (MQ=255)
                                                                                                                                   gCCTATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCt                    <  1:3945418/150‑1 (MQ=255)
                                                                                                                                   gCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGTTGGACGAGCTGGGCAATATGAGCGAGGAGGGTATCCCGGCGGCGCTGGTAGAAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCt                    <  1:127259/150‑1 (MQ=255)
                                                                                                                                   gCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCt                    >  1:2654700/1‑150 (MQ=255)
                                                                                                                                    ccGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCtt                   >  2:156707/1‑150 (MQ=255)
                                                                                                                                     cGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTta                  >  1:1821630/1‑150 (MQ=255)
                                                                                                                                     cGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGAGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTta                  >  1:649570/1‑150 (MQ=255)
                                                                                                                                      gATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtat                 >  2:2827557/1‑150 (MQ=255)
                                                                                                                                           tatGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACct            >  1:3303033/1‑150 (MQ=255)
                                                                                                                                            atGTTTCTCTATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctg           >  1:3024992/1‑150 (MQ=255)
                                                                                                                                            atGTTTCTCTATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCGTATAACctg           >  1:4070248/1‑150 (MQ=255)
                                                                                                                                            atGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctg           <  1:2200398/150‑1 (MQ=255)
                                                                                                                                            atGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctg           >  1:3126880/1‑150 (MQ=255)
                                                                                                                                            atGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctg           >  1:3823510/1‑150 (MQ=255)
                                                                                                                                             tttttCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgc          >  1:2972418/3‑150 (MQ=255)
                                                                                                                                              gTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgct         <  2:3750448/150‑1 (MQ=255)
                                                                                                                                              gTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgct         <  1:1902123/150‑1 (MQ=255)
                                                                                                                                                ttCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGCATCCAGGCGGCGCTGGTCGCAAAAATCCTCGACTAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGt       >  1:1408881/1‑150 (MQ=255)
                                                                                                                                                  ctcGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGttt     >  1:342085/1‑150 (MQ=255)
                                                                                                                                                  ctcGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGAAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGttt     <  2:3156684/150‑1 (MQ=255)
                                                                                                                                                   tcGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtc    <  2:888091/150‑1 (MQ=255)
                                                                                                                                                     gATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc  <  2:3914190/150‑1 (MQ=255)
                                                                                                                                                     gATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc  <  1:2980960/150‑1 (MQ=255)
                                                                                                                                                     gATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc  <  2:2130345/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGCCGAATGCTGGCCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTC  >  CP009273/207475‑207773

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: