Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,777,451 C→A 40.3% A55E (GCG→GAG)  fadK → short chain acyl‑CoA synthetase, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,777,4510CA40.3% 267.3 / 296.7 243A55E (GCG→GAG) fadKshort chain acyl‑CoA synthetase, anaerobic
Reads supporting (aligned to +/- strand):  ref base C (60/85);  new base A (46/52);  total (106/137)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 250 total aligned reads displayed.

AACGTTTAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGTGCCGC  >  CP009273/1777302‑1777600
                                                                                                                                                     |                                                                                                                                                     
aaCGTTTAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAgcgc                                                                                                                                                       >  2:1136696/1‑150 (MQ=255)
     ttAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCg                                                                                                                                                  >  2:609156/1‑150 (MQ=255)
     ttAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCg                                                                                                                                                  >  2:607770/1‑150 (MQ=255)
        aCGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATc                                                                                                                                               >  2:702441/1‑150 (MQ=255)
          gAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCAc                                                                                                                                             >  1:3424158/1‑150 (MQ=255)
          gAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAc                                                                                                                                             <  1:3330706/150‑1 (MQ=255)
          gAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAc                                                                                                                                             <  2:4620304/150‑1 (MQ=255)
           aacaacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTCTAGCGAGCTCGATCAcg                                                                                                                                            >  1:2102175/1‑150 (MQ=255)
             caacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcgcc                                                                                                                                          >  1:2593987/1‑150 (MQ=255)
                 gtcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGa                                                                                                                                      >  1:4655152/1‑150 (MQ=255)
                  tcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAg                                                                                                                                     <  2:1924042/150‑1 (MQ=255)
                     tGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAACATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTg                                                                                                                                  <  2:2072376/150‑1 (MQ=255)
                     tGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTg                                                                                                                                  <  2:2311575/150‑1 (MQ=255)
                     tGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTg                                                                                                                                  <  2:2938645/150‑1 (MQ=255)
                       cggcgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTCCACCTATAGCGCGCTCGATCACGCCGCGAGCTGTc                                                                                                                                >  1:529178/1‑150 (MQ=255)
                       cggcgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTc                                                                                                                                <  2:2058841/150‑1 (MQ=255)
                         gcgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTg                                                                                                                              <  1:3721382/150‑1 (MQ=255)
                         gcgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTg                                                                                                                              <  2:3449132/150‑1 (MQ=255)
                          cgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTgg                                                                                                                             <  1:2174024/150‑1 (MQ=255)
                          cgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTgg                                                                                                                             <  1:2173894/150‑1 (MQ=255)
                          cgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTgg                                                                                                                             <  2:668051/150‑1 (MQ=255)
                           gTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGc                                                                                                                            >  2:3989055/1‑150 (MQ=255)
                                gTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACt                                                                                                                       >  2:3125872/1‑150 (MQ=255)
                                  cAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTgg                                                                                                                     >  2:1645223/1‑150 (MQ=255)
                                   aGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGa                                                                                                                    >  1:2459516/1‑150 (MQ=255)
                                    gCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGAt                                                                                                                   <  2:1966805/150‑1 (MQ=255)
                                     cAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATg                                                                                                                  >  1:2521569/1‑150 (MQ=255)
                                      aaGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGt                                                                                                                 <  2:1727901/150‑1 (MQ=255)
                                       aGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGtt                                                                                                                <  1:2773564/150‑1 (MQ=255)
                                         ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAg                                                                                                              >  2:268196/1‑150 (MQ=255)
                                         ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAg                                                                                                              <  2:4370582/150‑1 (MQ=255)
                                         ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAg                                                                                                              >  2:1475647/1‑150 (MQ=255)
                                           ttATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCg                                                                                                            <  2:1540227/150‑1 (MQ=255)
                                             aTGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGaa                                                                                                          <  1:4772612/150‑1 (MQ=255)
                                               ggggCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAgg                                                                                                        <  2:4753505/150‑1 (MQ=255)
                                                gggCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAggg                                                                                                       >  1:948221/1‑150 (MQ=255)
                                                 ggCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGt                                                                                                      <  2:1012826/150‑1 (MQ=255)
                                                  gCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTa                                                                                                     <  1:287369/150‑1 (MQ=255)
                                                  gCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTa                                                                                                     >  2:3163210/1‑150 (MQ=255)
                                                       gCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGaa                                                                                                <  1:4491101/150‑1 (MQ=255)
                                                          tCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCa                                                                                             >  1:1617419/1‑150 (MQ=255)
                                                          tCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCa                                                                                             >  1:3558431/1‑150 (MQ=255)
                                                          tCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCa                                                                                             >  2:2931920/1‑150 (MQ=255)
                                                          tCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCa                                                                                             >  2:1823472/1‑150 (MQ=255)
                                                           cGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAg                                                                                            <  2:3335394/150‑1 (MQ=255)
                                                              tGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGAGGTTAGCGAAGGGTATTGAATCAGGCg                                                                                         >  2:3906130/1‑150 (MQ=255)
                                                                 ccGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatc                                                                                      <  1:2050418/150‑1 (MQ=255)
                                                                  cGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcg                                                                                     <  1:670741/150‑1 (MQ=255)
                                                                  cGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcg                                                                                     <  1:569730/150‑1 (MQ=255)
                                                                  cGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcg                                                                                     <  1:2978703/150‑1 (MQ=255)
                                                                   gATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgc                                                                                    >  2:2747967/1‑150 (MQ=255)
                                                                     ttACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgcat                                                                                  <  1:1246994/150‑1 (MQ=255)
                                                                       aCTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgcatcg                                                                                <  2:2161771/150‑1 (MQ=255)
                                                                        cTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgcatcgc                                                                               >  2:1283322/1‑150 (MQ=255)
                                                                          ggCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATcgcatcgcat                                                                             >  2:795518/1‑150 (MQ=255)
                                                                            cagcagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCAttt                                                                           <  2:1954605/150‑1 (MQ=255)
                                                                              gcagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCa                                                                         <  2:1590131/150‑1 (MQ=255)
                                                                              gcagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCa                                                                         <  1:2693687/150‑1 (MQ=255)
                                                                                agACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAAc                                                                       >  2:4611138/1‑150 (MQ=255)
                                                                                agACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAAc                                                                       >  2:2780732/1‑150 (MQ=255)
                                                                                     gCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCct                                                                  <  2:1786429/150‑1 (MQ=255)
                                                                                        cGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCctggc                                                               >  1:3918016/1‑150 (MQ=255)
                                                                                             gATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtg                                                          <  1:149403/150‑1 (MQ=255)
                                                                                              aTGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtgt                                                         <  1:1901490/150‑1 (MQ=255)
                                                                                              aTGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtgt                                                         <  1:4569434/150‑1 (MQ=255)
                                                                                              aTGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtgt                                                         <  2:3582729/150‑1 (MQ=255)
                                                                                                 ccAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGaa                                                      >  2:1360960/1‑150 (MQ=255)
                                                                                                 ccAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGaa                                                      <  2:648741/150‑1 (MQ=255)
                                                                                                 ccAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGaa                                                      >  1:1913511/1‑150 (MQ=255)
                                                                                                   aGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAtt                                                    <  2:2341042/150‑1 (MQ=255)
                                                                                                   aGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAtt                                                    <  1:709820/150‑1 (MQ=255)
                                                                                                     aCAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTa                                                  >  1:3191560/1‑150 (MQ=255)
                                                                                                     aCAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTa                                                  >  1:2942415/1‑150 (MQ=255)
                                                                                                     aCAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTa                                                  <  2:1056522/150‑1 (MQ=255)
                                                                                                      cAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTAc                                                 <  1:3473779/150‑1 (MQ=255)
                                                                                                       aaaaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTAcc                                                <  2:3243000/150‑1 (MQ=255)
                                                                                                       aaaaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTAcc                                                >  1:2802119/1‑150 (MQ=255)
                                                                                                         aaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGt                                              >  1:3847669/1‑150 (MQ=255)
                                                                                                          agtGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGtt                                             <  2:3415500/148‑1 (MQ=255)
                                                                                                          aTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGtt                                             <  1:1360960/150‑1 (MQ=255)
                                                                                                           ttGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTta                                            <  1:3625854/150‑1 (MQ=255)
                                                                                                             gCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatc                                          >  1:1105392/1‑150 (MQ=255)
                                                                                                                 tGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatctatc                                      >  1:1948807/1‑150 (MQ=255)
                                                                                                                   gTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCtt                                    >  2:2205583/1‑150 (MQ=255)
                                                                                                                      gATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCTGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcc                                 >  1:2011251/1‑150 (MQ=255)
                                                                                                                      gATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcc                                 >  2:4167769/1‑150 (MQ=255)
                                                                                                                           tCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTgcctgcct                            <  1:4292563/150‑1 (MQ=255)
                                                                                                                           tCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTgcctgcct                            <  2:2459516/150‑1 (MQ=255)
                                                                                                                           tCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTgcctgcct                            <  2:2593987/150‑1 (MQ=255)
                                                                                                                              tGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGaa                         >  2:4648006/1‑150 (MQ=255)
                                                                                                                                  gCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATc                     <  2:4230974/150‑1 (MQ=255)
                                                                                                                                   cATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCg                    >  1:553314/1‑150 (MQ=255)
                                                                                                                                   cATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCg                    <  1:815462/150‑1 (MQ=255)
                                                                                                                                      cGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTg                 <  2:2521530/150‑1 (MQ=255)
                                                                                                                                      cGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTg                 >  2:4462303/1‑150 (MQ=255)
                                                                                                                                          cacCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGt             >  1:2439911/1‑150 (MQ=255)
                                                                                                                                           acCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGtt            <  1:1687908/150‑1 (MQ=255)
                                                                                                                                           acCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGtt            >  1:3717740/1‑150 (MQ=255)
                                                                                                                                              tataGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTcc         <  2:3632626/150‑1 (MQ=255)
                                                                                                                                                 aGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGTg      >  1:4390547/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACGTTTAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGTGCCGC  >  CP009273/1777302‑1777600

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: