Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,055,901 A→T 15.0% Q397L (CAA→CTA)  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,9010AT15.0% 629.1 / 93.4 214Q397L (CAA→CTA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (79/103);  new base T (16/16);  total (95/119)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 220 total aligned reads displayed.

CCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCA  >  CP009273/1055754‑1056050
                                                                                                                                                   |                                                                                                                                                     
ccGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCtaa                                                                                                                                                     <  1:3135124/150‑1 (MQ=255)
   gCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTg                                                                                                                                                  >  1:4824493/1‑150 (MQ=255)
     aTTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGc                                                                                                                                                >  2:2677250/1‑150 (MQ=255)
      ttGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCTAATTGGcc                                                                                                                                               <  1:1403081/150‑1 (MQ=255)
        gATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTg                                                                                                                                             <  2:1758170/150‑1 (MQ=255)
         aTGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCACATTGGCCTGc                                                                                                                                            >  2:3168708/1‑150 (MQ=255)
           gCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCa                                                                                                                                          <  2:3521441/150‑1 (MQ=255)
               aaaCCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATGGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAggtg                                                                                                                                      <  2:2952597/150‑1 (MQ=255)
                 aCCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAggtggt                                                                                                                                    <  1:265044/150‑1 (MQ=255)
                    aTTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAggtggtggt                                                                                                                                 >  1:5116463/1‑150 (MQ=255)
                     ttCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGtggtt                                                                                                                                >  2:804960/1‑150 (MQ=255)
                     ttCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGtggtt                                                                                                                                <  1:4015715/150‑1 (MQ=255)
                      tCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGtggttt                                                                                                                               <  2:4007088/150‑1 (MQ=255)
                      tCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGtggttt                                                                                                                               >  1:1789915/1‑150 (MQ=255)
                       cGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGtggtttt                                                                                                                              <  2:1538361/150‑1 (MQ=255)
                        gTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGtggttttg                                                                                                                             <  2:2294739/150‑1 (MQ=255)
                        gTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGtggttttg                                                                                                                             <  2:4462741/150‑1 (MQ=255)
                            gCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGACCAGTGGGCGTGGAGGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGgggttttggttt                                                                                                                         >  2:75478/1‑150 (MQ=255)
                                  gcgGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTtggctg                                                                                                                   <  1:2980708/150‑1 (MQ=255)
                                   cgGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTTTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTtggctgg                                                                                                                  >  1:5168225/1‑150 (MQ=255)
                                     gCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCa                                                                                                                <  2:414847/150‑1 (MQ=255)
                                      cAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCTAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCAc                                                                                                               <  1:4086662/150‑1 (MQ=255)
                                            gcggcgAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCTAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAAc                                                                                                         <  2:3615716/150‑1 (MQ=255)
                                            gcggcgAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAAc                                                                                                         <  1:4247178/150‑1 (MQ=255)
                                            gcggcgAACAGAACGCAAATTAGTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAAc                                                                                                         <  2:4655987/150‑1 (MQ=255)
                                             cggcgAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACg                                                                                                        <  2:4704709/150‑1 (MQ=255)
                                                cgAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCTAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGgcg                                                                                                     >  2:1135851/1‑150 (MQ=255)
                                                 gAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGgcgc                                                                                                    >  1:2754008/1‑150 (MQ=255)
                                                     aGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGc                                                                                                >  2:1050158/1‑150 (MQ=255)
                                                       aaCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCTAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCAc                                                                                              <  1:3722669/150‑1 (MQ=255)
                                                        aCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACg                                                                                             >  2:1500656/1‑150 (MQ=255)
                                                        aCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACg                                                                                             >  1:869448/1‑150 (MQ=255)
                                                        aCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACg                                                                                             <  1:2870100/150‑1 (MQ=255)
                                                          gCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGcc                                                                                           <  1:1735722/150‑1 (MQ=255)
                                                           cAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCg                                                                                          <  2:4251042/150‑1 (MQ=255)
                                                           cAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCg                                                                                          >  1:4555031/1‑150 (MQ=255)
                                                             aaTTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCggg                                                                                        <  1:1566121/150‑1 (MQ=255)
                                                             aaTTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTGTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCggg                                                                                        >  2:3604743/1‑150 (MQ=255)
                                                               ttattGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGc                                                                                      <  1:2591301/150‑1 (MQ=255)
                                                                   tGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTaa                                                                                  <  2:622175/150‑1 (MQ=255)
                                                                      tggctggTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAgg                                                                               <  1:804960/150‑1 (MQ=255)
                                                                      tgcctggTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAgg                                                                               >  2:1829253/4‑150 (MQ=255)
                                                                       ggctggTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGc                                                                              >  1:2698178/1‑150 (MQ=255)
                                                                           ggTGCGTGCAGCGTATGCAGCACGGTGACCAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTa                                                                          >  1:2819719/1‑150 (MQ=255)
                                                                           ggTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTa                                                                          >  1:2330852/1‑150 (MQ=255)
                                                                                 tgcAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGa                                                                    <  2:2457541/150‑1 (MQ=255)
                                                                                 tgcAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAAAGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGa                                                                    >  1:2379871/1‑150 (MQ=255)
                                                                                    aGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTg                                                                 >  1:2421752/1‑150 (MQ=255)
                                                                                      cGTATGCAGCACGGTGAACAGTGGGCGGGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGGTATTCTGAGTGGt                                                               >  2:3419312/1‑150 (MQ=255)
                                                                                       gTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGtt                                                              <  2:2351876/150‑1 (MQ=255)
                                                                                       gTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGtt                                                              >  1:3603095/1‑150 (MQ=255)
                                                                                       gTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGtt                                                              >  1:3447543/1‑150 (MQ=255)
                                                                                        tATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGttt                                                             <  2:3399990/150‑1 (MQ=255)
                                                                                        tATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGttt                                                             <  2:4981760/150‑1 (MQ=255)
                                                                                        tATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGttt                                                             <  2:5189620/150‑1 (MQ=255)
                                                                                         aTGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTtc                                                            <  1:1703070/150‑1 (MQ=255)
                                                                                          tGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTtct                                                           >  1:3554247/1‑150 (MQ=255)
                                                                                            cagcaCGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCTAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTcc                                                         >  1:2405002/1‑150 (MQ=255)
                                                                                            cagcaCGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTcc                                                         <  1:3338172/150‑1 (MQ=255)
                                                                                            cagcaCGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTcc                                                         <  1:754639/150‑1 (MQ=255)
                                                                                             agcaCGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAAGTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACACCGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCg                                                        <  2:3612142/150‑1 (MQ=255)
                                                                                                  ggTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGctct                                                   <  2:274199/150‑1 (MQ=255)
                                                                                                   gTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTa                                                  >  1:2624574/1‑150 (MQ=255)
                                                                                                   gTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTa                                                  <  2:4777065/150‑1 (MQ=255)
                                                                                                       aCAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTAc                                                 <  1:4353552/147‑1 (MQ=255)
                                                                                                       aCAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTAc                                                 >  2:4353552/1‑147 (MQ=255)
                                                                                                       aCAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTc                                              <  1:4325685/150‑1 (MQ=255)
                                                                                                            gggCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCTAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTc                                         >  1:2263923/1‑150 (MQ=255)
                                                                                                            gggCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTc                                         <  2:200692/150‑1 (MQ=255)
                                                                                                             ggCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTcc                                        <  1:42042/150‑1 (MQ=255)
                                                                                                             ggCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTcc                                        <  2:1208793/150‑1 (MQ=255)
                                                                                                             ggCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTcc                                        >  2:5031104/1‑150 (MQ=255)
                                                                                                              gCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCg                                       >  2:5123173/1‑150 (MQ=255)
                                                                                                               cGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGc                                      >  2:3082234/1‑150 (MQ=255)
                                                                                                                gTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGcc                                     >  1:222915/1‑150 (MQ=255)
                                                                                                                 tGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCt                                    <  1:4253525/150‑1 (MQ=255)
                                                                                                                   gatgatTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGt                                  <  1:4742359/150‑1 (MQ=255)
                                                                                                                     tgatTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTc                                <  2:4885179/150‑1 (MQ=255)
                                                                                                                       atTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCAc                              <  1:4519094/150‑1 (MQ=255)
                                                                                                                        ttGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACg                             <  1:941637/150‑1 (MQ=255)
                                                                                                                           tgGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGaca                          >  2:3954114/1‑150 (MQ=255)
                                                                                                                              ttCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGacaaca                       <  2:3049096/150‑1 (MQ=255)
                                                                                                                                 tGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTg                    <  1:4451206/150‑1 (MQ=255)
                                                                                                                                  ggcggcGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGa                   <  1:1969143/150‑1 (MQ=255)
                                                                                                                                    cggcgATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACt                 >  2:3689988/1‑150 (MQ=255)
                                                                                                                                     ggcgATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTa                >  2:1737542/1‑150 (MQ=255)
                                                                                                                                       cgATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACa              <  1:4129182/150‑1 (MQ=255)
                                                                                                                                        gATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACaa             <  2:4777578/150‑1 (MQ=255)
                                                                                                                                        gATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACaa             <  2:2191814/150‑1 (MQ=255)
                                                                                                                                         aTGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACaaa            >  1:4546163/1‑150 (MQ=255)
                                                                                                                                          tGCTGGGGCTAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAg           >  1:4689574/1‑150 (MQ=255)
                                                                                                                                           gCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAgg          <  2:3325269/150‑1 (MQ=255)
                                                                                                                                            cTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGc         <  1:322502/150‑1 (MQ=255)
                                                                                                                                             tGGGGCTAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCt        >  2:2595330/1‑150 (MQ=255)
                                                                                                                                             tGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCt        <  2:460405/150‑1 (MQ=255)
                                                                                                                                               gggCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTAc      <  1:3574326/150‑1 (MQ=255)
                                                                                                                                                ggCAAATTGGCCTTCCATGTGGTGGTTTTTGTTTTGGCTGGCACTACCACGGCGCATGGACGCCGGGGGGTAAAGGCGGTTTTCTGAGTGGTTTCTCCGGCTCTACGTCGTTTCCGCCTGTTCAAGTCAACAGTGTATACAAAGGCTAca     <  2:4198055/150‑1 (MQ=255)
                                                                                                                                                ggCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTAca     <  1:4066554/150‑1 (MQ=255)
                                                                                                                                                   aaaTTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTAcagca  >  1:2126670/1‑150 (MQ=255)
                                                                                                                                                   aaaTTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTAcagca  >  1:3552147/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCA  >  CP009273/1055754‑1056050

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: