Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,777,451 C→A 100% 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,4510CA98.1% 844.7 / 5.3 259A55E (GCG→GAG) fadKshort chain acyl‑CoA synthetase, anaerobic
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (127/127);  total (130/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 263 total aligned reads displayed.

CGTTTAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGTGC  >  CP009273/1777304‑1777597
                                                                                                                                                   |                                                                                                                                                  
cGTTTAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAgcgagc                                                                                                                                                  >  2:4358825/1‑150 (MQ=255)
     aaCGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGAt                                                                                                                                             >  2:5129549/1‑150 (MQ=255)
        gAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAc                                                                                                                                          <  1:4466932/150‑1 (MQ=255)
          acaacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcgc                                                                                                                                        >  2:3283791/1‑150 (MQ=255)
           caacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcgcc                                                                                                                                       >  2:3566984/1‑150 (MQ=255)
           caacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcgcc                                                                                                                                       >  1:2468222/1‑150 (MQ=255)
             acGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCTCTGGCCGATTCCTGGCAGCAGACCGCTCGTGCGATGCCAGACCAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcgccgc                                                                                                                                     <  2:2954254/150‑1 (MQ=255)
             acGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcgccgc                                                                                                                                     <  1:3283791/150‑1 (MQ=255)
              cgtcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCcgcg                                                                                                                                    <  1:483969/150‑1 (MQ=255)
              cgtcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCcgcg                                                                                                                                    <  1:4698446/150‑1 (MQ=255)
               gtcgtGCGGCGTATCGTCAGCACGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGa                                                                                                                                   <  2:1449121/150‑1 (MQ=255)
               gtcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGa                                                                                                                                   <  2:2496482/150‑1 (MQ=255)
               gtcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGa                                                                                                                                   >  1:3012060/1‑150 (MQ=255)
               gtcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGa                                                                                                                                   <  1:4273881/150‑1 (MQ=255)
                 tgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGc                                                                                                                                 >  2:4306551/2‑150 (MQ=255)
                   tGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTg                                                                                                                               >  1:2546617/1‑150 (MQ=255)
                     tggcgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTc                                                                                                                             >  2:613619/2‑150 (MQ=255)
                      ggcgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCt                                                                                                                            >  2:4682079/1‑150 (MQ=255)
                       gcgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTg                                                                                                                           <  2:2697898/150‑1 (MQ=255)
                        cgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTgg                                                                                                                          >  2:1114681/1‑150 (MQ=255)
                        cgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTgg                                                                                                                          <  1:2005940/150‑1 (MQ=255)
                         gTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGc                                                                                                                         >  1:1603362/1‑150 (MQ=255)
                           aTCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCaa                                                                                                                       <  2:4685957/150‑1 (MQ=255)
                            tCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCaaa                                                                                                                      <  1:3179739/150‑1 (MQ=255)
                            tCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCaaa                                                                                                                      >  1:5128759/1‑150 (MQ=255)
                              gTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACt                                                                                                                    <  2:3511485/150‑1 (MQ=255)
                               tCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTg                                                                                                                   >  2:4101653/1‑150 (MQ=255)
                                 aGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGa                                                                                                                 >  2:4034050/1‑150 (MQ=255)
                                 aGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGa                                                                                                                 >  2:3781336/1‑150 (MQ=255)
                                      gggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTa                                                                                                            <  1:688829/150‑1 (MQ=255)
                                       ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAg                                                                                                           >  1:2747791/1‑150 (MQ=255)
                                       ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAg                                                                                                           <  1:2167610/150‑1 (MQ=255)
                                       ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAg                                                                                                           >  1:1271298/1‑150 (MQ=255)
                                       ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAg                                                                                                           >  1:3199680/1‑150 (MQ=255)
                                         ttATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCg                                                                                                         <  2:442215/150‑1 (MQ=255)
                                            tGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAg                                                                                                      >  2:1334079/1‑150 (MQ=255)
                                               ggCGATTTTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGt                                                                                                   >  2:4567345/1‑150 (MQ=255)
                                               ggCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTTTGGCAAACTGGATGTTAGCGAAGGGt                                                                                                   >  2:4039221/1‑150 (MQ=255)
                                                 cGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTAt                                                                                                 <  1:1005066/150‑1 (MQ=255)
                                                  gATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTCGCGAAGGGTAtt                                                                                                >  2:1751019/1‑150 (MQ=255)
                                                  gATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTAtt                                                                                                >  1:1900774/1‑150 (MQ=255)
                                                  gATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTAtt                                                                                                <  1:2234059/150‑1 (MQ=255)
                                                    tGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGa                                                                                              <  2:2697456/150‑1 (MQ=255)
                                                          gCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAgg                                                                                        <  2:3979782/150‑1 (MQ=255)
                                                          gCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAgg                                                                                        <  2:3804478/150‑1 (MQ=255)
                                                           cTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGc                                                                                       <  1:1298710/150‑1 (MQ=255)
                                                               ccGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatc                                                                                   <  2:2859629/150‑1 (MQ=255)
                                                               ccGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatc                                                                                   <  1:3717267/150‑1 (MQ=255)
                                                                 gATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgc                                                                                 >  2:3351666/1‑150 (MQ=255)
                                                                  aTTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgca                                                                                <  2:5047671/150‑1 (MQ=255)
                                                                      cTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgcatcgc                                                                            >  1:2892887/1‑150 (MQ=255)
                                                                       tGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGAtcgcatcgca                                                                           <  1:887302/150‑1 (MQ=255)
                                                                       tGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGAtcgcatcgca                                                                           <  1:4997031/150‑1 (MQ=255)
                                                                       tGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGAtcgcatcgca                                                                           <  1:4267052/150‑1 (MQ=255)
                                                                         gcagcaGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCAtt                                                                         <  1:544601/150‑1 (MQ=255)
                                                                           agcagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTc                                                                       >  1:2448474/1‑150 (MQ=255)
                                                                            gcagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCACACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCAGTTCa                                                                      >  2:4804733/1‑150 (MQ=255)
                                                                             cagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCaa                                                                     >  2:886406/1‑150 (MQ=255)
                                                                               gACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACt                                                                   <  2:3836786/150‑1 (MQ=255)
                                                                               gACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACt                                                                   >  1:843049/1‑150 (MQ=255)
                                                                                 ccGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGc                                                                 >  2:3197989/1‑150 (MQ=255)
                                                                                   gCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCct                                                               <  1:434943/150‑1 (MQ=255)
                                                                                   gCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCct                                                               >  2:1224358/1‑150 (MQ=255)
                                                                                    cTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTAGAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAATCTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCctg                                                              <  2:5072473/150‑1 (MQ=255)
                                                                                          cGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgt                                                        <  1:3071999/150‑1 (MQ=255)
                                                                                           gATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtg                                                       >  2:256874/1‑150 (MQ=255)
                                                                                           gATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtg                                                       <  2:3019368/150‑1 (MQ=255)
                                                                                             tGCCAGACAAACTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGtgtg                                                     <  2:2595685/150‑1 (MQ=255)
                                                                                               ccAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGaa                                                   <  1:3566408/150‑1 (MQ=255)
                                                                                               ccAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGaa                                                   >  1:3619725/1‑150 (MQ=255)
                                                                                               ccAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGaa                                                   >  2:4267321/1‑150 (MQ=255)
                                                                                               caaGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACAACTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGaa                                                   >  1:3818158/3‑150 (MQ=255)
                                                                                                cAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAt                                                  <  1:591039/150‑1 (MQ=255)
                                                                                                 aGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAtt                                                 >  2:1827990/1‑150 (MQ=255)
                                                                                                  gACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAttt                                                <  1:4753519/150‑1 (MQ=255)
                                                                                                  gACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAttt                                                >  2:944000/1‑150 (MQ=255)
                                                                                                   aCAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTa                                               <  2:2480977/150‑1 (MQ=255)
                                                                                                   aCAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTa                                               >  2:2153972/1‑150 (MQ=255)
                                                                                                    cAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTAc                                              >  2:1565305/1‑150 (MQ=255)
                                                                                                    cAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTAc                                              >  2:2512375/1‑150 (MQ=255)
                                                                                                     aaaaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTAcc                                             <  1:4725077/150‑1 (MQ=255)
                                                                                                       aaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGt                                           <  1:3934785/150‑1 (MQ=255)
                                                                                                        aTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGtt                                          <  2:4087322/150‑1 (MQ=255)
                                                                                                          tGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtat                                        >  2:4514549/1‑150 (MQ=255)
                                                                                                              gTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatctat                                    <  1:4233389/150‑1 (MQ=255)
                                                                                                                 gTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCtt                                 <  1:4596959/150‑1 (MQ=255)
                                                                                                                    gATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcc                              <  1:1558840/150‑1 (MQ=255)
                                                                                                                      tAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcctg                            <  1:2351428/150‑1 (MQ=255)
                                                                                                                        aTCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcctgcc                          >  2:1515972/1‑150 (MQ=255)
                                                                                                                         tCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcctgccc                         <  2:1932561/150‑2 (MQ=255)
                                                                                                                          cATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGcctgcctg                        >  1:3391452/1‑150 (MQ=255)
                                                                                                                             ggTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGaaa                     >  2:2167185/1‑150 (MQ=255)
                                                                                                                              gTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGaaaa                    >  2:2291332/1‑150 (MQ=255)
                                                                                                                                gCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATc                  <  1:1166917/150‑1 (MQ=255)
                                                                                                                                   tCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGt               >  2:357504/1‑150 (MQ=255)
                                                                                                                                    cGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTg              <  2:3279745/150‑1 (MQ=255)
                                                                                                                                    cGTACACCTAGAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTg              <  2:3704504/150‑1 (MQ=255)
                                                                                                                                             ataGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCg     >  1:4890372/1‑150 (MQ=255)
                                                                                                                                              taGCGAGCTCGATCACGCCGCGAGCTGTCTGTCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGt    <  2:180545/150‑1 (MQ=255)
                                                                                                                                                gcgagcTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGTGc  >  2:1902789/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTTTAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGTGC  >  CP009273/1777304‑1777597

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: