Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,055,617 A→T 8.3% A302A (GCA→GCT torA → trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,055,6170AT 8.3% 898.3 / 56.0 277A302A (GCA→GCTtorAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (117/137);  new base T (13/10);  total (130/147)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 285 total aligned reads displayed.

CAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGC  >  CP009273/1055468‑1055766
                                                                                                                                                     |                                                                                                                                                     
   cTAAAAGCGACAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACAt                                                                                                                                                    <  1:1566952/150‑1 (MQ=255)
    tAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATa                                                                                                                                                   >  2:2334763/1‑150 (MQ=255)
     aaaaGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATAc                                                                                                                                                  <  1:635714/150‑1 (MQ=255)
      aaaGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACg                                                                                                                                                 >  1:2822381/1‑150 (MQ=255)
        aGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCt                                                                                                                                               >  1:438796/1‑150 (MQ=255)
            aaaGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTAc                                                                                                                                           >  2:3233162/1‑150 (MQ=255)
             aaGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACa                                                                                                                                          >  1:2669975/1‑150 (MQ=255)
              aGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAg                                                                                                                                         >  1:1057195/1‑150 (MQ=255)
                  gccgccGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGaa                                                                                                                                     >  1:1417180/1‑150 (MQ=255)
                   ccgccgGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGaaa                                                                                                                                    <  2:2645596/150‑1 (MQ=255)
                   ccgccgGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGaaa                                                                                                                                    <  1:2347228/150‑1 (MQ=255)
                      ccgGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAAcc                                                                                                                                 <  1:2138304/150‑1 (MQ=255)
                       cgGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCt                                                                                                                                >  1:4551864/1‑150 (MQ=255)
                          tGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTa                                                                                                                             <  1:3875714/150‑1 (MQ=255)
                           gAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTAc                                                                                                                            <  2:2706858/150‑1 (MQ=255)
                             aaTTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGa                                                                                                                          >  1:1422463/1‑150 (MQ=255)
                                tGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACaa                                                                                                                       >  1:4007643/1‑150 (MQ=255)
                                tGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACaa                                                                                                                       <  2:3984524/150‑1 (MQ=255)
                                 gAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACaaa                                                                                                                      <  2:4507480/150‑1 (MQ=255)
                                      catcaGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACtt                                                                                                                 >  1:3071931/1‑150 (MQ=255)
                                       atcaGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTc                                                                                                                <  2:2017522/150‑1 (MQ=255)
                                         caGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCt                                                                                                              <  2:4288924/150‑1 (MQ=255)
                                         caGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCt                                                                                                              <  2:3271395/150‑1 (MQ=255)
                                          aGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCTCATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCtt                                                                                                             >  2:2811100/1‑150 (MQ=255)
                                          aGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCtt                                                                                                             <  1:3292071/150‑1 (MQ=255)
                                           gCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTg                                                                                                            <  2:2970114/150‑1 (MQ=255)
                                            cATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGc                                                                                                           <  2:4400156/150‑1 (MQ=255)
                                              tCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTa                                                                                                         <  2:411661/150‑1 (MQ=255)
                                                 aTCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAact                                                                                                      >  1:3044822/1‑150 (MQ=255)
                                                 aTCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAact                                                                                                      >  1:382062/1‑150 (MQ=255)
                                                  tCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCTCATACGCTGTACAGTGAAAACCTGTACGACAAAATCTTCCTTGCTAacta                                                                                                     >  2:3730128/1‑150 (MQ=255)
                                                    cGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAactact                                                                                                   <  2:3424335/150‑1 (MQ=255)
                                                    cGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAactact                                                                                                   >  1:1840676/1‑150 (MQ=255)
                                                    cGGTTGTCACATCCACCCATGAGTATCTGGGGCGAGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGACGCTGCAACTGGAGCTGGCACATACGCTGTACAGTGAAAACCTGTACGAAATTAACTTCCTTGATAactaat                                                                                                   <  2:1675708/150‑3 (MQ=255)
                                                     ggTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACtg                                                                                                  >  2:3148079/1‑150 (MQ=255)
                                                      gTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACtgt                                                                                                 >  1:3450573/1‑150 (MQ=255)
                                                        tGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACtgtgt                                                                                               <  1:1984294/150‑1 (MQ=255)
                                                        tGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACtgtgt                                                                                               <  2:2052346/150‑1 (MQ=255)
                                                                ccACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTg                                                                                       >  1:2770785/1‑150 (MQ=255)
                                                                  aCCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAg                                                                                     >  2:505453/1‑150 (MQ=255)
                                                                      aTGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCTCATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGt                                                                                 >  1:1076571/1‑150 (MQ=255)
                                                                       tGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGtt                                                                                <  1:3299767/150‑1 (MQ=255)
                                                                        gAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTc                                                                               <  2:1898735/150‑1 (MQ=255)
                                                                          gTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCACTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCt                                                                             <  1:3964302/150‑1 (MQ=255)
                                                                           tATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTg                                                                            <  2:2926563/150‑1 (MQ=255)
                                                                               tGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGt                                                                        >  1:1459404/1‑150 (MQ=255)
                                                                               tGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGt                                                                        >  2:154056/1‑150 (MQ=255)
                                                                                 gggCGCGAGCATGAGAAGCACATTGCGGTTAACCCGCAAACTGACGAGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTTCTGTGTGGGTTTTGAGCAGTTCCTGCCGTAt                                                                      <  2:2993610/150‑1 (MQ=255)
                                                                                   gcgcGAGCATGTGCAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGATGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTTCGACAAAAACTTCCTTGCTAACTAATGTGTGGGTTTTGAGCAGTTCCTGCCGTATct                                                                    <  2:2254079/150‑1 (MQ=255)
                                                                                   gcgcGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATct                                                                    <  2:1214049/150‑1 (MQ=255)
                                                                                      cgAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATctgct                                                                 <  1:2538127/150‑1 (MQ=255)
                                                                                      cgAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATctgct                                                                 <  2:3404655/150‑1 (MQ=255)
                                                                                       gAGCATGTGAAGCACATAGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGATGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTGTTGAGCAGTTCCTGCCGTATctgctg                                                                >  2:3879995/1‑150 (MQ=255)
                                                                                        aGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTgg                                                               >  1:5817/1‑150 (MQ=255)
                                                                                            tgtgAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCTCATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTga                                                           >  2:1559901/1‑150 (MQ=255)
                                                                                              tgAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTgaga                                                         >  2:3529649/1‑150 (MQ=255)
                                                                                              tgAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTgaga                                                         >  2:1721327/1‑150 (MQ=255)
                                                                                              tgAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTgaga                                                         >  2:3046935/1‑150 (MQ=255)
                                                                                               gAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGaa                                                        >  2:2108546/1‑150 (MQ=255)
                                                                                               gAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGaa                                                        >  1:4143274/1‑150 (MQ=255)
                                                                                                aaGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGaaa                                                       >  2:487761/1‑150 (MQ=255)
                                                                                                aaGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGGGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGaaa                                                       >  2:4452386/1‑150 (MQ=255)
                                                                                                 aGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAg                                                      >  1:1950595/1‑150 (MQ=255)
                                                                                                 aGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAg                                                      >  2:2896009/1‑150 (MQ=255)
                                                                                                  gCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGa                                                     <  1:1292577/150‑1 (MQ=255)
                                                                                                    acaTTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACg                                                   >  2:3675586/1‑150 (MQ=255)
                                                                                                    acaTTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACg                                                   >  1:3239119/1‑150 (MQ=255)
                                                                                                      aTTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCTCATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGt                                                 <  2:235041/150‑1 (MQ=255)
                                                                                                           ggTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGcc                                            <  2:1015559/150‑1 (MQ=255)
                                                                                                            gTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCg                                           >  2:758953/1‑150 (MQ=255)
                                                                                                             ttAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCTCATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGa                                          <  1:1007875/150‑1 (MQ=255)
                                                                                                               aaCCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGaaa                                        >  2:3195913/1‑150 (MQ=255)
                                                                                                                aCCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAg                                       <  2:2673589/150‑1 (MQ=255)
                                                                                                                aCCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAg                                       <  2:918814/150‑1 (MQ=255)
                                                                                                                   cGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCTCATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATg                                    >  2:1348647/1‑150 (MQ=255)
                                                                                                                          tGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATg                             <  2:3576550/150‑1 (MQ=255)
                                                                                                                          tGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATg                             >  2:961939/1‑150 (MQ=255)
                                                                                                                          tGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGAGGGTTTTGAGCAGTTCCTGCCTTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCCTg                             >  1:4080464/1‑150 (MQ=255)
                                                                                                                            aCGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATggg                           <  1:151832/150‑1 (MQ=255)
                                                                                                                            aCGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATggg                           >  1:367660/1‑150 (MQ=255)
                                                                                                                              gTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCt                         <  1:1273725/150‑1 (MQ=255)
                                                                                                                               tGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTg                        >  2:4057973/1‑150 (MQ=255)
                                                                                                                                gCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGa                       <  2:4595449/150‑1 (MQ=255)
                                                                                                                                   gCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGGTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGaaaa                    >  2:4222012/1‑150 (MQ=255)
                                                                                                                                     tGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGTAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAAc                  <  1:843482/150‑1 (MQ=255)
                                                                                                                                      gCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACt                 >  2:4196831/1‑150 (MQ=255)
                                                                                                                                       cAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTg                >  1:4173013/1‑150 (MQ=255)
                                                                                                                                         aCTGGCGCTGGCTCATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGAc              <  1:525676/150‑1 (MQ=255)
                                                                                                                                          cTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGAcc             <  2:596049/150‑1 (MQ=255)
                                                                                                                                          cTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGGTCCTGCCGTATCTGCTGGGTGAGACAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGAcc             <  2:1640243/150‑1 (MQ=255)
                                                                                                                                           tGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCg            <  1:1311350/150‑1 (MQ=255)
                                                                                                                                           tGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCg            <  2:2396908/150‑1 (MQ=255)
                                                                                                                                            ggCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCgg           <  1:2296521/150‑1 (MQ=255)
                                                                                                                                              cgcTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCa         <  2:2340935/150‑1 (MQ=255)
                                                                                                                                               gcTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCAt        <  1:171774/150‑1 (MQ=255)
                                                                                                                                               gcTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCAt        <  2:3588052/150‑1 (MQ=255)
                                                                                                                                                 tGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTg      <  2:506856/150‑1 (MQ=255)
                                                                                                                                                   gCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGAt    <  2:2822340/150‑1 (MQ=255)
                                                                                                                                                    cacaTACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATg   >  1:3115145/1‑150 (MQ=255)
                                                                                                                                                     acaTACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAAATGACCGGCATGTATGc  >  2:1706147/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGC  >  CP009273/1055468‑1055766

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: