Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,777,451 C→A 74.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,4510CA74.3% 470.6 / 267.3 311A55E (GCG→GAG) fadKshort chain acyl‑CoA synthetase, anaerobic
Reads supporting (aligned to +/- strand):  ref base C (33/47);  new base A (76/155);  total (109/202)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 311 total aligned reads displayed.

TTTAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGTG  >  CP009273/1777306‑1777596
                                                                                                                                                 |                                                                                                                                                 
tttAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTc                                                                                                                                               <  2:583757/150‑1 (MQ=255)
  tAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGa                                                                                                                                             <  2:4025252/150‑1 (MQ=255)
   aaCGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGAt                                                                                                                                            >  1:2248044/1‑150 (MQ=255)
      gAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAc                                                                                                                                         <  1:3470988/150‑1 (MQ=255)
      gAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAc                                                                                                                                         <  1:3471072/150‑1 (MQ=255)
         caacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCAcgcc                                                                                                                                      >  1:3407052/1‑150 (MQ=255)
          aacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCAcgccg                                                                                                                                     >  2:3420586/1‑150 (MQ=255)
          aacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcgccg                                                                                                                                     <  2:632671/150‑1 (MQ=255)
            cgtcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCcgcg                                                                                                                                   <  1:2594064/150‑1 (MQ=255)
            cgtcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCcgcg                                                                                                                                   <  1:2593717/150‑1 (MQ=255)
                gtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCt                                                                                                                               <  1:1722283/150‑1 (MQ=255)
                    ggcgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCt                                                                                                                           >  2:150576/1‑150 (MQ=255)
                      cgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTgg                                                                                                                         >  1:893164/1‑150 (MQ=255)
                        tATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCa                                                                                                                       <  1:1239012/150‑1 (MQ=255)
                           cGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAAc                                                                                                                    >  1:2105756/1‑150 (MQ=255)
                            gTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACt                                                                                                                   <  2:3138359/150‑1 (MQ=255)
                            gTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACt                                                                                                                   <  2:856797/150‑1 (MQ=255)
                            gTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACt                                                                                                                   <  1:1296768/150‑1 (MQ=255)
                             tCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTg                                                                                                                  <  1:865636/150‑1 (MQ=255)
                              cAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTgg                                                                                                                 >  2:3737243/1‑150 (MQ=255)
                              cAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTgg                                                                                                                 >  2:1472096/1‑150 (MQ=255)
                               aGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGa                                                                                                                >  2:1363323/1‑150 (MQ=255)
                               aGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGa                                                                                                                <  1:1344061/150‑1 (MQ=255)
                                   aGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGtt                                                                                                            <  2:2208327/150‑1 (MQ=255)
                                     ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAg                                                                                                          >  1:1405512/1‑150 (MQ=255)
                                     ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAg                                                                                                          <  2:685495/150‑1 (MQ=255)
                                      gTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGc                                                                                                         <  2:1271234/150‑1 (MQ=255)
                                      gTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGc                                                                                                         <  1:3824015/150‑1 (MQ=255)
                                           ggggCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAgg                                                                                                    <  1:607189/150‑1 (MQ=255)
                                            gggCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGGTCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAggg                                                                                                   >  1:4398864/1‑150 (MQ=255)
                                             ggCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGt                                                                                                  <  1:670070/150‑1 (MQ=255)
                                               cGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTAt                                                                                                <  2:1768924/150‑1 (MQ=255)
                                                gATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTAtt                                                                                               >  2:3239223/1‑150 (MQ=255)
                                                gATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTAtt                                                                                               >  1:505770/1‑150 (MQ=255)
                                                 aTGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTg                                                                                              <  2:2576860/150‑1 (MQ=255)
                                                     ttCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATc                                                                                          <  1:4468752/150‑1 (MQ=255)
                                                      tCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCa                                                                                         <  1:4204907/150‑1 (MQ=255)
                                                      tCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCa                                                                                         <  2:4365894/150‑1 (MQ=255)
                                                      tCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCa                                                                                         <  1:1984921/150‑1 (MQ=255)
                                                       cGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAg                                                                                        >  1:4103822/1‑150 (MQ=255)
                                                            gCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGat                                                                                   <  2:4304151/150‑1 (MQ=255)
                                                             ccGATTACTGGCAGCATACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGAACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatc                                                                                  >  1:3207805/1‑150 (MQ=255)
                                                             ccGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatc                                                                                  <  2:239241/150‑1 (MQ=255)
                                                             ccGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatc                                                                                  <  1:289296/150‑1 (MQ=255)
                                                               gATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgc                                                                                >  1:577759/1‑150 (MQ=255)
                                                                 ttACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgcat                                                                              <  1:2409806/150‑1 (MQ=255)
                                                                  tACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgcatc                                                                             <  1:2705761/150‑1 (MQ=255)
                                                                   aCTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgcatcg                                                                            <  1:3287378/150‑1 (MQ=255)
                                                                      ggCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATcgcatcgcat                                                                         >  1:3763114/1‑150 (MQ=255)
                                                                        cagcagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCAttt                                                                       <  2:400280/150‑1 (MQ=255)
                                                                        cagcagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCAttt                                                                       <  2:400770/150‑1 (MQ=255)
                                                                        cagcagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCAttt                                                                       <  2:2688444/150‑1 (MQ=255)
                                                                        cagcagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGACTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCAttt                                                                       <  2:3648026/150‑1 (MQ=255)
                                                                           cagACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGACGGGTATTGAATCAGGCGATCGCATCGCATTTCaa                                                                    >  2:673336/1‑150 (MQ=255)
                                                                            agACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAAc                                                                   >  2:138706/1‑150 (MQ=255)
                                                                            agACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAAc                                                                   <  2:347520/150‑1 (MQ=255)
                                                                                 gCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCct                                                              <  2:4209850/150‑1 (MQ=255)
                                                                                      tGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCctggctg                                                         <  1:1142557/150‑1 (MQ=255)
                                                                                      tGCGAGGCCAGACAAAAGTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCctggctg                                                         <  1:1674017/150‑1 (MQ=255)
                                                                                       gCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCctggctgg                                                        >  2:1628350/1‑150 (MQ=255)
                                                                                        cGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgt                                                       <  2:581056/150‑1 (MQ=255)
                                                                                        cGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAGCTGCCTGGCTGgt                                                       <  2:4343608/150‑1 (MQ=255)
                                                                                         gATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtg                                                      >  2:424662/1‑150 (MQ=255)
                                                                                             ccAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGaa                                                  >  2:3790771/1‑150 (MQ=255)
                                                                                             ccAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGaa                                                  <  2:2666397/150‑1 (MQ=255)
                                                                                               aGACAAAATTGCCGTGGTCGATAATCATGGTGCATAGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAtt                                                <  1:1917401/150‑1 (MQ=255)
                                                                                                gACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAttt                                               <  1:2312085/150‑1 (MQ=255)
                                                                                                gACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAttt                                               <  1:2312228/150‑1 (MQ=255)
                                                                                                   aaaaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTAcc                                            >  2:2201758/1‑150 (MQ=255)
                                                                                                    aaaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCg                                           <  2:4222248/150‑1 (MQ=255)
                                                                                                      aTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGtt                                         <  2:2380583/150‑1 (MQ=255)
                                                                                                         gCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatc                                      <  2:4525169/150‑1 (MQ=255)
                                                                                                         gCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatc                                      <  1:2386063/150‑1 (MQ=255)
                                                                                                          ccGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatct                                     <  2:3591545/150‑1 (MQ=255)
                                                                                                           cGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatcta                                    <  1:697664/150‑1 (MQ=255)
                                                                                                            gTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatctat                                   >  2:2858876/1‑150 (MQ=255)
                                                                                                            gTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatctat                                   >  2:1005711/1‑150 (MQ=255)
                                                                                                              ggTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTatctatct                                 >  1:3993582/1‑150 (MQ=255)
                                                                                                                tCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtg                               <  2:2085923/150‑1 (MQ=255)
                                                                                                                  gATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcc                             <  2:3846011/150‑1 (MQ=255)
                                                                                                                  gATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcc                             >  2:1965210/1‑150 (MQ=255)
                                                                                                                         aTGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGa                      <  2:2286008/150‑1 (MQ=255)
                                                                                                                           ggTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGaaa                    >  1:3810682/1‑150 (MQ=255)
                                                                                                                            gTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGaaaa                   >  2:2091934/1‑150 (MQ=255)
                                                                                                                              gCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATc                 <  2:3648927/150‑1 (MQ=255)
                                                                                                                              gCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATc                 <  1:3664488/150‑1 (MQ=255)
                                                                                                                                 tCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGt              <  2:720103/150‑1 (MQ=255)
                                                                                                                                 tCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGt              <  2:3000987/150‑1 (MQ=255)
                                                                                                                                  cGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTg             >  2:526316/1‑150 (MQ=255)
                                                                                                                                  cGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTg             >  2:1850633/1‑150 (MQ=255)
                                                                                                                                    tACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCa           <  2:3957911/150‑1 (MQ=255)
                                                                                                                                     acacCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAg          >  1:3657590/1‑150 (MQ=255)
                                                                                                                                      cacCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGt         <  1:1549678/150‑1 (MQ=255)
                                                                                                                                      cacCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGt         <  1:1549632/150‑1 (MQ=255)
                                                                                                                                      cacCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGt         <  1:1094550/150‑1 (MQ=255)
                                                                                                                                       acCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGtt        <  2:4290248/150‑1 (MQ=255)
                                                                                                                                        ccTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGttt       <  1:3988978/150‑1 (MQ=255)
                                                                                                                                         cTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTc      >  2:1457027/1‑150 (MQ=255)
                                                                                                                                           ataGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCg    <  1:3233825/150‑1 (MQ=255)
                                                                                                                                           ataGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCg    >  2:4307133/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTTAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGTG  >  CP009273/1777306‑1777596

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: