Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,777,451 C→A 91.8% 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,4510CA91.8% 500.1 / 37.8 184A55E (GCG→GAG) fadKshort chain acyl‑CoA synthetase, anaerobic
Reads supporting (aligned to +/- strand):  ref base C (9/6);  new base A (72/97);  total (81/103)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 190 total aligned reads displayed.

TTTAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGTG  >  CP009273/1777306‑1777596
                                                                                                                                                 |                                                                                                                                                 
tttAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAgcgagc                                                                                                                                                 <  2:2500852/148‑1 (MQ=255)
tttAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAgcgagc                                                                                                                                                 <  1:1133563/148‑1 (MQ=255)
  tAACGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTc                                                                                                                                               >  1:96134/1‑148 (MQ=255)
     cGAACAACGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGAt                                                                                                                                            >  1:347848/1‑148 (MQ=255)
         caacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcg                                                                                                                                        <  1:1132240/148‑1 (MQ=255)
         caacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcg                                                                                                                                        <  2:1730231/148‑1 (MQ=255)
         caacGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcg                                                                                                                                        >  2:285024/1‑148 (MQ=255)
           acGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCAcgcc                                                                                                                                      >  2:1282095/1‑148 (MQ=255)
           acGTCGTGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcgcc                                                                                                                                      >  1:2547809/1‑148 (MQ=255)
            cgtcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCAcgccg                                                                                                                                     >  1:3287485/1‑148 (MQ=255)
              tcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCcgcg                                                                                                                                   <  2:1232999/148‑1 (MQ=255)
              tcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCcgcg                                                                                                                                   >  1:1242840/1‑148 (MQ=255)
              tcgtGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCcgcg                                                                                                                                   >  1:296932/1‑148 (MQ=255)
                 tGCGGCGTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGc                                                                                                                                >  2:1153474/1‑148 (MQ=255)
                   cggcgTATCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTg                                                                                                                              >  2:1469031/1‑148 (MQ=255)
                          tCGTCAGCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCa                                                                                                                       >  1:2668767/1‑148 (MQ=255)
                                gCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTgg                                                                                                                 <  2:434203/148‑1 (MQ=255)
                                gCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTgg                                                                                                                 >  2:2918187/1‑148 (MQ=255)
                                gCAAGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTgg                                                                                                                 <  2:1893458/148‑1 (MQ=255)
                                   aGGGTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATg                                                                                                              <  1:808715/148‑1 (MQ=255)
                                    gggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGt                                                                                                             <  1:235082/148‑1 (MQ=255)
                                    gggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGt                                                                                                             <  1:1073470/148‑1 (MQ=255)
                                     ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGtt                                                                                                            <  1:1377396/148‑1 (MQ=255)
                                     ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGtt                                                                                                            <  2:1947243/148‑1 (MQ=255)
                                     ggTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGtt                                                                                                            <  1:1377956/148‑1 (MQ=255)
                                      gTTATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTa                                                                                                           <  2:1119434/148‑1 (MQ=255)
                                       ttATGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAg                                                                                                          <  1:1933701/148‑1 (MQ=255)
                                          tGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGa                                                                                                       >  1:2781181/1‑148 (MQ=255)
                                          tGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGa                                                                                                       >  1:1005553/1‑148 (MQ=255)
                                          tGGGGCGATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGa                                                                                                       <  2:529916/148‑1 (MQ=255)
                                                gATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTa                                                                                                 <  2:1186491/148‑1 (MQ=255)
                                                gATGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTa                                                                                                 <  2:2837173/148‑1 (MQ=255)
                                                 aTGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTAt                                                                                                >  2:2716707/1‑148 (MQ=255)
                                                 aTGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTAt                                                                                                >  1:1725376/1‑148 (MQ=255)
                                                  tGCTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTAtt                                                                                               >  1:3039798/1‑148 (MQ=255)
                                                    cTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGa                                                                                             <  1:2342031/148‑1 (MQ=255)
                                                    cTTCGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGa                                                                                             <  1:2051422/148‑1 (MQ=255)
                                                       cGCTGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATc                                                                                          <  1:2918325/148‑1 (MQ=255)
                                                          tGGCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAgg                                                                                       <  2:243710/148‑1 (MQ=255)
                                                           ggCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTTCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGc                                                                                      >  1:1364582/1‑148 (MQ=255)
                                                            gCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCg                                                                                     >  2:2105529/1‑148 (MQ=255)
                                                            gCCGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCg                                                                                     >  2:15812/1‑148 (MQ=255)
                                                             ccGATTACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGa                                                                                    <  2:2305801/148‑1 (MQ=255)
                                                                  tACTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgca                                                                               <  2:2802111/148‑1 (MQ=255)
                                                                   aCTGGCAGCAGACCGCTCGTGCGATTACAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCTAACTGTCTGTCAAACTTGATGTTAGAGAAGGGTATTGAATCAGGCGatcgcat                                                                              <  2:2974045/148‑1 (MQ=255)
                                                                   aCTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgcat                                                                              <  2:911265/148‑1 (MQ=255)
                                                                   aCTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgcat                                                                              <  2:2370738/148‑1 (MQ=255)
                                                                   aCTGGCAGCAGACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGatcgcat                                                                              >  2:62849/1‑148 (MQ=255)
                                                                       gcagcaGACCGCTCGTGCGATGCCAGACTATATTGCAGTGGTCGATAATCTTGGTGCATCGTACACATATAGCGAGCTCGCTCACGCCGCGAGCTGTATGGCAAAGTGGATGTTAGCGAAGGGTATTGAATCAGGCGAtcgcatcgca                                                                          <  1:360027/148‑1 (MQ=255)
                                                                             gACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCaa                                                                    <  2:983244/148‑1 (MQ=255)
                                                                             gACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCaa                                                                    <  2:149586/148‑1 (MQ=255)
                                                                             gACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCaa                                                                    <  2:1820723/148‑1 (MQ=255)
                                                                             gACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCaa                                                                    <  1:449684/148‑1 (MQ=255)
                                                                             gACCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCaa                                                                    <  1:2906346/148‑1 (MQ=255)
                                                                              aCCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAAc                                                                   >  2:2646939/1‑148 (MQ=255)
                                                                              aCCGCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAAc                                                                   >  1:3131786/1‑148 (MQ=255)
                                                                                 gCTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGc                                                                <  2:3284956/148‑1 (MQ=255)
                                                                                  cTCGTGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGcc                                                               <  2:2519147/148‑1 (MQ=255)
                                                                                      tGCGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCctggc                                                           <  2:2739320/148‑1 (MQ=255)
                                                                                        cGATGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCctggctg                                                         >  1:2390970/1‑148 (MQ=255)
                                                                                           tGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtg                                                      >  2:1288345/1‑148 (MQ=255)
                                                                                           tGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtg                                                      <  1:2097292/148‑1 (MQ=255)
                                                                                           tGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtg                                                      >  2:3121940/1‑148 (MQ=255)
                                                                                           tGCCAGACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGgtg                                                      >  2:1381244/1‑148 (MQ=255)
                                                                                                gACAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAt                                                 <  1:1657838/148‑1 (MQ=255)
                                                                                                 aCAAAATTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAAtt                                                >  2:3248617/1‑148 (MQ=255)
                                                                                                   aaaaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTa                                              >  1:860946/1‑148 (MQ=255)
                                                                                                   aaaaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTa                                              <  2:575050/148‑1 (MQ=255)
                                                                                                    aaaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTAc                                             <  2:1233490/148‑1 (MQ=255)
                                                                                                     aaTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTAcc                                            >  1:1750338/1‑148 (MQ=255)
                                                                                                      aTTGCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCg                                           >  1:2552586/1‑148 (MQ=255)
                                                                                                         gCCGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTta                                        <  1:1479264/148‑1 (MQ=255)
                                                                                                          ccGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtat                                       <  1:2118808/148‑1 (MQ=255)
                                                                                                          ccGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtat                                       <  1:141444/148‑1 (MQ=255)
                                                                                                           cGTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatc                                      >  1:320130/1‑148 (MQ=255)
                                                                                                            gTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatct                                     <  1:1001522/148‑1 (MQ=255)
                                                                                                            gTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatct                                     >  2:1510431/1‑148 (MQ=255)
                                                                                                            gTGGTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatct                                     <  2:1414414/148‑1 (MQ=255)
                                                                                                            gTGGTCGAGAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatct                                     <  2:2708892/148‑1 (MQ=255)
                                                                                                              ggTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatctat                                   <  1:2186032/148‑1 (MQ=255)
                                                                                                              ggTCGATAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatctat                                   <  1:2817465/148‑1 (MQ=255)
                                                                                                               gTCGATAATCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTtatctatc                                  <  1:1282160/148‑1 (MQ=255)
                                                                                                                   aTAATCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgc                              >  1:2630090/1‑148 (MQ=255)
                                                                                                                     aaTCATGGTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcct                            >  1:895870/1‑148 (MQ=255)
                                                                                                                       tCATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcctgc                          >  1:1328843/1‑148 (MQ=255)
                                                                                                                        cATGGTGCATCGTACACCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTtgcctgcc                         >  1:3126952/1‑148 (MQ=255)
                                                                                                                            gTGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGaa                     >  1:3310603/1‑148 (MQ=255)
                                                                                                                             tGCATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGaaa                    >  1:1150419/1‑148 (MQ=255)
                                                                                                                               cATCGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAAt                  <  2:96134/148‑1 (MQ=255)
                                                                                                                                  cGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCgg               >  2:1317395/1‑148 (MQ=255)
                                                                                                                                  cGTACACCTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCgg               <  2:2547809/148‑1 (MQ=255)
                                                                                                                                       acCTATAGCGCGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAg          <  2:1161035/148‑1 (MQ=255)
                                                                                                                                        ccTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGt         <  1:62849/148‑1 (MQ=255)
                                                                                                                                         cTATAGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGtt        >  1:76941/1‑148 (MQ=255)
                                                                                                                                          tataGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGttt       <  1:2441014/148‑1 (MQ=255)
                                                                                                                                          tataGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGttt       >  1:2267557/1‑148 (MQ=255)
                                                                                                                                            taGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGttt       <  2:762654/146‑1 (MQ=255)
                                                                                                                                             aGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCg    >  1:2281703/1‑148 (MQ=255)
                                                                                                                                             aGCGAGCTCGATCACGCCGCGAGCTGTCTGGCAAACTGGATGTTAGCGAAGGGTATTGAATCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCg    <  1:3104877/148‑1 (MQ=255)
                                                                                                                                               cgcgcTCGATCACGCCGCGAGCTGTCTGGCAAACTGGCTGTTAGCGAAGGGTATTGACTCAGGCGATCGCATCGCATTTCAACTGCCTGGCTGGTGTGAATTTACCGTTATCTATCTTGCCTGCCTGAAAATCGGTGCAGTTTCCGTg  <  2:3030884/148‑1 (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: