Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,136,368 G→T 48.4% R386S (CGT→AGT)  hybB ← putative hydrogenase 2 cytochrome b type component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,136,3680GT48.4% 136.1 / inf 289R386S (CGT→AGT) hybBputative hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base G (57/92);  new base T (60/80);  total (117/172)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 290 total aligned reads displayed.

TCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATG  >  CP009273/3136221‑3136517
                                                                                                                                                   |                                                                                                                                                     
tctTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGAt                                                                                                                                                     <  2:4436164/150‑1 (MQ=255)
tctTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGAt                                                                                                                                                     >  2:3692416/1‑150 (MQ=255)
tctTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGAt                                                                                                                                                     <  1:1931822/150‑1 (MQ=255)
tctTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGAt                                                                                                                                                     <  1:1932277/150‑1 (MQ=255)
   tCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTAt                                                                                                                                                  <  2:3303264/150‑1 (MQ=255)
      aTGACGCGCCACATGGTACCGGAAGCCCATGCTTACGAAACGACGCCATATTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGAt                                                                                                                                               >  2:2260587/1‑150 (MQ=255)
         ccGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgt                                                                                                                                            <  1:2996795/150‑1 (MQ=255)
           gcgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgttt                                                                                                                                          <  2:3534246/150‑1 (MQ=255)
           gcgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCGTTGCTCGCCTCATGACTATTATGATCgttt                                                                                                                                          <  1:4461487/150‑1 (MQ=255)
            cgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCgtttt                                                                                                                                         >  2:3906421/1‑150 (MQ=255)
            cgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtttt                                                                                                                                         <  1:3379589/150‑1 (MQ=255)
            cgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtttt                                                                                                                                         <  2:3664541/150‑1 (MQ=255)
            cgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtttt                                                                                                                                         <  2:3479074/150‑1 (MQ=255)
            cgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtttt                                                                                                                                         <  1:4188637/150‑1 (MQ=255)
            cgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtttt                                                                                                                                         <  1:4256133/150‑1 (MQ=255)
                acaTGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCgttttgttt                                                                                                                                     <  1:3643355/150‑1 (MQ=255)
                  aTGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTa                                                                                                                                   >  1:633274/1‑150 (MQ=255)
                   tGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTaa                                                                                                                                  <  1:3288662/150‑1 (MQ=255)
                   tGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTaa                                                                                                                                  <  2:1044980/150‑1 (MQ=255)
                    ggTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTaaa                                                                                                                                 <  2:4169376/150‑1 (MQ=255)
                       aCCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAagga                                                                                                                              >  2:2176592/1‑150 (MQ=255)
                       aCCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAagga                                                                                                                              >  2:25189/1‑150 (MQ=255)
                         cGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAaggagg                                                                                                                            >  2:2856515/1‑150 (MQ=255)
                            aaGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAg                                                                                                                         <  1:157278/150‑1 (MQ=255)
                            aaGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAg                                                                                                                         <  2:3109715/150‑1 (MQ=255)
                             aGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGt                                                                                                                        <  2:3867707/150‑1 (MQ=255)
                               cccATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTAt                                                                                                                      <  2:874101/150‑1 (MQ=255)
                               cccATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTAt                                                                                                                      <  2:3314328/150‑1 (MQ=255)
                               cccATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTAt                                                                                                                      <  2:2959360/150‑1 (MQ=255)
                                 cATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCg                                                                                                                    <  1:3201000/150‑1 (MQ=255)
                                         cGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGa                                                                                                            <  1:1952018/150‑1 (MQ=255)
                                         cGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCAAAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATAATGATCTTTTTGTTTTAAAGGAGGAAGTATCGGCGGTAGa                                                                                                            >  2:3485184/1‑150 (MQ=255)
                                          gAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGAc                                                                                                           <  2:498465/150‑1 (MQ=255)
                                          gAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGAc                                                                                                           <  2:498527/150‑1 (MQ=255)
                                           aaaCGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACg                                                                                                          >  2:1072102/1‑150 (MQ=255)
                                               gacgCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATg                                                                                                      <  2:486141/150‑1 (MQ=255)
                                                   ccATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTa                                                                                                  <  1:3076037/150‑1 (MQ=255)
                                                     aTTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATGCGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACg                                                                                                >  1:2800862/1‑150 (MQ=255)
                                                         tCGATTTCGCAATCGATGCGTAAATGCCCCTCATTACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTATAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGt                                                                                            <  1:2260587/150‑1 (MQ=255)
                                                          cGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTa                                                                                           >  2:1549270/1‑150 (MQ=255)
                                                          cGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTa                                                                                           >  1:4139376/1‑150 (MQ=255)
                                                              ttCGCAATCGATGCGTAAATTCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAgcg                                                                                       <  2:2964127/150‑1 (MQ=255)
                                                                cGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCa                                                                                     <  1:1899827/150‑1 (MQ=255)
                                                                cGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCa                                                                                     <  1:906498/150‑1 (MQ=255)
                                                                cGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCa                                                                                     <  1:222394/150‑1 (MQ=255)
                                                                 gCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCaa                                                                                    <  1:3741941/150‑1 (MQ=255)
                                                                 gCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCaa                                                                                    >  1:3166783/1‑150 (MQ=255)
                                                                  cAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGATTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCaaa                                                                                   >  2:1785999/1‑150 (MQ=255)
                                                                   aaTCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAt                                                                                  <  1:683654/150‑1 (MQ=255)
                                                                    aTCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAtc                                                                                 <  2:977705/150‑1 (MQ=255)
                                                                    aTCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAtc                                                                                 <  2:833509/150‑1 (MQ=255)
                                                                      cGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAtctc                                                                               >  1:3052272/1‑150 (MQ=255)
                                                                        aTGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCaa                                                                             >  1:2022196/1‑150 (MQ=255)
                                                                        aTGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCaa                                                                             >  1:3271630/1‑150 (MQ=255)
                                                                         tGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAAt                                                                            >  1:514111/1‑150 (MQ=255)
                                                                           cGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAg                                                                          <  1:3260867/150‑1 (MQ=255)
                                                                             tAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGcc                                                                        <  1:1056658/150‑1 (MQ=255)
                                                                             tAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGcc                                                                        >  2:1494312/1‑150 (MQ=255)
                                                                               aaTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCcac                                                                      >  2:1855617/1‑150 (MQ=255)
                                                                                aTGCCCCTCTATACGGGTTACCGGATCAATATTAAGTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTTTGATCGTTTTGTTTTATAGGAGGTAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCcaca                                                                     <  2:661098/150‑1 (MQ=255)
                                                                                aTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCcaca                                                                     >  2:763338/1‑150 (MQ=255)
                                                                                 tGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACaa                                                                    >  2:214897/1‑150 (MQ=255)
                                                                                 tGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACaa                                                                    >  1:2810479/1‑150 (MQ=255)
                                                                                   ccccTCAATAATGGTTACCTGATCATTAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAATGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTATATAGCCACaaaa                                                                  <  2:2682120/150‑1 (MQ=255)
                                                                                    cccTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAAc                                                                 >  2:2869057/1‑150 (MQ=255)
                                                                                    cccTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAAc                                                                 >  2:585139/1‑150 (MQ=255)
                                                                                    cccTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAAc                                                                 <  1:1318977/150‑1 (MQ=255)
                                                                                      cTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCa                                                               <  1:2430786/150‑1 (MQ=255)
                                                                                                   aCCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACa                                                  >  2:3930151/1‑150 (MQ=255)
                                                                                                    ccGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAg                                                 >  1:1349877/1‑150 (MQ=255)
                                                                                                      ggATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGtt                                               <  2:1016266/150‑1 (MQ=255)
                                                                                                      ggATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGtt                                               >  2:1193769/1‑150 (MQ=255)
                                                                                                       gATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGttc                                              <  1:2841726/150‑1 (MQ=255)
                                                                                                        aTCTATAGTTATTCTCTGGCTCATGCTTTGCTCGCATCATGACTATTATGATAGTTTTTTTTTAAAGGAGGAATTATCGGCAGTAGACGAATGAGTACGATGTAAGAGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGttct                                             <  1:2512717/150‑1 (MQ=255)
                                                                                                         tCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGttctt                                            <  1:4344105/150‑1 (MQ=255)
                                                                                                         tCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGttctt                                            <  1:4532664/150‑1 (MQ=255)
                                                                                                             tAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCa                                        <  2:1369235/150‑1 (MQ=255)
                                                                                                             tAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCa                                        >  1:1727263/1‑150 (MQ=255)
                                                                                                               gTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAgg                                      >  1:2828307/1‑150 (MQ=255)
                                                                                                               gTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAgg                                      <  1:1843637/150‑1 (MQ=255)
                                                                                                                  aTTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCg                                   <  2:2800862/150‑1 (MQ=255)
                                                                                                                  aTTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCCCAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCg                                   >  2:4228684/1‑150 (MQ=255)
                                                                                                                   ttCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCTGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATATCCACAAAACCAATAGAAATCAACAGATCTACACAGGTCgg                                  >  2:3776452/1‑150 (MQ=255)
                                                                                                                   ttCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCgg                                  <  2:1407559/150‑1 (MQ=255)
                                                                                                                      tctGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGaa                               <  2:1520384/150‑1 (MQ=255)
                                                                                                                         ggCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTa                            <  2:2321328/150‑1 (MQ=255)
                                                                                                                          gCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAg                           >  1:3377747/1‑150 (MQ=255)
                                                                                                                             cATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCg                        >  2:1740877/1‑150 (MQ=255)
                                                                                                                               tGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTa                      >  2:1987376/1‑150 (MQ=255)
                                                                                                                                  tttGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAcc                   <  1:1193769/150‑1 (MQ=255)
                                                                                                                                    tGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAccgc                 >  1:1427076/1‑150 (MQ=255)
                                                                                                                                     gCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAccgcc                >  2:1298618/1‑150 (MQ=255)
                                                                                                                                          ccTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCg           >  1:1729703/1‑150 (MQ=255)
                                                                                                                                             cATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACTATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGATCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGt        <  1:3692419/150‑1 (MQ=255)
                                                                                                                                             cATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGt        <  1:3486169/150‑1 (MQ=255)
                                                                                                                                              aTGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGtt       <  2:4124398/150‑1 (MQ=255)
                                                                                                                                              aTGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGtt       <  2:3647752/150‑1 (MQ=255)
                                                                                                                                              aTGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGtt       <  1:2739382/150‑1 (MQ=255)
                                                                                                                                              aTGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGtt       <  1:1684994/150‑1 (MQ=255)
                                                                                                                                               tGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTg      <  1:2973105/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATG  >  CP009273/3136221‑3136517

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: