Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,604,731 C→T 5.0% R223W (CGG→TGG)  hyfI → hydrogenase 4, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,604,7310CT 5.0% 582.8 / 11.1 158R223W (CGG→TGG) hyfIhydrogenase 4, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (67/83);  new base T (4/4);  total (71/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

ACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGT  >  CP009273/2604587‑2604875
                                                                                                                                                |                                                                                                                                                
aCTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCggg                                                                                                                                               <  2:3192696/148‑1 (MQ=255)
aCTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCggg                                                                                                                                               >  1:2489832/1‑148 (MQ=255)
   aTGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTa                                                                                                                                            >  1:2890758/1‑148 (MQ=255)
   aTGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTa                                                                                                                                            >  1:2520393/1‑148 (MQ=255)
    tGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGCGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTaa                                                                                                                                           <  1:772288/148‑1 (MQ=255)
    tGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTaa                                                                                                                                           <  1:772294/148‑1 (MQ=255)
     gCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAac                                                                                                                                          >  1:1495716/1‑148 (MQ=255)
       aaCCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAacac                                                                                                                                        <  1:2898717/148‑1 (MQ=255)
         ccGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCt                                                                                                                                      >  1:2375701/1‑148 (MQ=255)
            ttGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGt                                                                                                                                   <  1:1842305/148‑1 (MQ=255)
              gtgGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTg                                                                                                                                 <  2:1197153/148‑1 (MQ=255)
              gtgGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTg                                                                                                                                 <  1:3100943/148‑1 (MQ=255)
              gtgGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTg                                                                                                                                 <  1:3100896/148‑1 (MQ=255)
                   cGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGa                                                                                                                            <  2:1202781/148‑1 (MQ=255)
                   cGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGa                                                                                                                            <  2:2344381/148‑1 (MQ=255)
                     cAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATg                                                                                                                          <  1:2073744/148‑1 (MQ=255)
                            cGCCATCACTGCGTAGCGCATTTGATCGAGACGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAATTTGGCGATGGGCTCATGCGCCATTTAAGCGCCGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCcgacga                                                                                                                   <  1:633892/148‑1 (MQ=255)
                            cGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCcgacga                                                                                                                   <  1:2296175/148‑1 (MQ=255)
                            cGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCcgacga                                                                                                                   <  1:3030306/148‑1 (MQ=255)
                            cGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCcgacga                                                                                                                   <  1:1788737/148‑1 (MQ=255)
                             gCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGAt                                                                                                                  <  2:3217205/148‑1 (MQ=255)
                               cATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATcc                                                                                                                >  2:2743250/1‑148 (MQ=255)
                                 tCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCAc                                                                                                              >  1:452207/1‑148 (MQ=255)
                                 tCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCAc                                                                                                              >  1:3137497/1‑148 (MQ=255)
                                       cGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCa                                                                                                        <  1:2245726/148‑1 (MQ=255)
                                        gTATCGCCATTGCGCGAGAAGCGCGGCGGCGGGCGGGATATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCaa                                                                                                       <  2:2244648/148‑1 (MQ=255)
                                         tATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCaag                                                                                                      >  2:2913764/1‑147 (MQ=255)
                                         tATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAAt                                                                                                      >  2:2913827/1‑148 (MQ=255)
                                         tATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAAt                                                                                                      >  1:1764369/1‑148 (MQ=255)
                                         tATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAAt                                                                                                      <  1:2846754/148‑1 (MQ=255)
                                                 ttGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGt                                                                                              >  1:1122770/1‑148 (MQ=255)
                                                 ttGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGt                                                                                              >  1:2696593/1‑148 (MQ=255)
                                                  tGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGtt                                                                                             <  2:2489832/148‑1 (MQ=255)
                                                  tGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGtt                                                                                             <  2:1334055/148‑1 (MQ=255)
                                                     gCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAg                                                                                          <  2:1374489/148‑1 (MQ=255)
                                                        agaAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCaa                                                                                       <  1:141592/148‑1 (MQ=255)
                                                        agaAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCaa                                                                                       >  2:2100334/1‑148 (MQ=255)
                                                          aaGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAAct                                                                                     >  1:798423/1‑148 (MQ=255)
                                                          aaGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAAct                                                                                     <  2:3099114/148‑1 (MQ=255)
                                                             cgcgGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCtt                                                                                  >  2:1351778/1‑148 (MQ=255)
                                                             cgcgGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCtt                                                                                  >  1:775942/1‑148 (MQ=255)
                                                              gcggcgGCTGGCGGGATATGGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCACCTCttt                                                                                 <  2:582075/148‑1 (MQ=255)
                                                                ggcggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTcg                                                                               <  1:1847157/148‑1 (MQ=255)
                                                                    gcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTa                                                                           <  2:2262707/148‑1 (MQ=255)
                                                                    gcTGGCGGGCTAGCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTa                                                                           <  2:1978207/148‑1 (MQ=255)
                                                                           ggCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGAGTACTCCggg                                                                    >  2:1877905/1‑148 (MQ=255)
                                                                             cTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCggggg                                                                  <  1:1015945/148‑1 (MQ=255)
                                                                             cTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCggggg                                                                  <  2:3137497/148‑1 (MQ=255)
                                                                                   tCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACAt                                                            >  1:2093959/1‑148 (MQ=255)
                                                                                    cAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATg                                                           >  2:1326116/1‑148 (MQ=255)
                                                                                     aGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGa                                                          >  1:961096/1‑148 (MQ=255)
                                                                                     aGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTGAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGa                                                          <  2:2133209/148‑1 (MQ=255)
                                                                                      gtggagagaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTAAAGGAAATCGGGTTAACACCTGGTTGCGCGATGCAGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGAc                                                         >  2:543338/6‑148 (MQ=255)
                                                                                       gggCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACt                                                        >  1:1955293/1‑148 (MQ=255)
                                                                                          cGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGaa                                                     >  1:3130703/1‑148 (MQ=255)
                                                                                            agaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAaga                                                   <  1:91132/148‑1 (MQ=255)
                                                                                            agaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAaga                                                   >  1:1501909/1‑148 (MQ=255)
                                                                                             gaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAagag                                                  <  1:1204492/148‑1 (MQ=255)
                                                                                             gaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAagag                                                  <  1:1223665/148‑1 (MQ=255)
                                                                                              aaaTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGt                                                 <  2:2834396/148‑1 (MQ=255)
                                                                                              aaaTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGt                                                 >  1:2084478/1‑148 (MQ=255)
                                                                                                aTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTcg                                                                               <  2:2322369/116‑1 (MQ=255)
                                                                                                aTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGc                                               <  2:1069458/148‑1 (MQ=255)
                                                                                                 tttGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCg                                              <  2:1505051/148‑1 (MQ=255)
                                                                                                  ttGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCgg                                             >  1:1009669/1‑148 (MQ=255)
                                                                                                  ttGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCgg                                             >  1:2364537/1‑148 (MQ=255)
                                                                                                    gCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCgggg                                           <  2:2468142/148‑1 (MQ=255)
                                                                                                    gCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCgggg                                           <  2:2467439/148‑1 (MQ=255)
                                                                                                    gCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCgggg                                           >  1:1859141/1‑148 (MQ=255)
                                                                                                     cGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGa                                          >  1:2999984/1‑148 (MQ=255)
                                                                                                      gATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGaa                                         >  2:2621423/1‑148 (MQ=255)
                                                                                                        tCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAt                                       <  1:404743/148‑1 (MQ=255)
                                                                                                           gCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGt                                    <  1:121559/148‑1 (MQ=255)
                                                                                                            cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtt                                   <  1:2696944/148‑1 (MQ=255)
                                                                                                            cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtt                                   <  2:2999984/148‑1 (MQ=255)
                                                                                                            cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtt                                   >  2:3019491/1‑148 (MQ=255)
                                                                                                              ccTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtttt                                 <  1:2073665/148‑1 (MQ=255)
                                                                                                               cTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTc                                >  1:885040/1‑148 (MQ=255)
                                                                                                                tGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCt                               <  2:2663161/148‑1 (MQ=255)
                                                                                                                tGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCt                               <  2:1008782/148‑1 (MQ=255)
                                                                                                                  cgcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTgg                             <  2:414245/148‑1 (MQ=255)
                                                                                                                  cgcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTgg                             <  2:414180/148‑1 (MQ=255)
                                                                                                                    ccATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCGGGTTGCGCGATGCCGAAGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGACGAGTGCGGGGAAATTGTTTTCTGGAc                           >  2:2213627/1‑148 (MQ=255)
                                                                                                                      aTTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGc                         <  2:2032235/148‑1 (MQ=255)
                                                                                                                      aTTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGc                         >  2:1907319/1‑148 (MQ=255)
                                                                                                                       tttAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCt                        >  1:1804172/1‑148 (MQ=255)
                                                                                                                            gCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGACATTGTTTTCTGGACGCTGCGaa                   >  1:2036356/1‑148 (MQ=255)
                                                                                                                            gCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGaa                   >  2:1523886/1‑148 (MQ=255)
                                                                                                                            gCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGaa                   >  1:1674832/1‑148 (MQ=255)
                                                                                                                                  aTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGt             <  2:1804172/148‑1 (MQ=255)
                                                                                                                                   tCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGtt            >  1:2711430/1‑148 (MQ=255)
                                                                                                                                    ccTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGttt           >  2:3089787/1‑148 (MQ=255)
                                                                                                                                    ccTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGttt           >  2:3229925/1‑148 (MQ=255)
                                                                                                                                          ggAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGcc     <  2:798423/148‑1 (MQ=255)
                                                                                                                                          ggAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGcc     >  1:2151450/1‑148 (MQ=255)
                                                                                                                                          ggAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGcc     >  1:1577907/1‑148 (MQ=255)
                                                                                                                                           gAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTTCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCa    >  2:3184021/1‑148 (MQ=255)
                                                                                                                                           gAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCa    >  2:2075188/1‑148 (MQ=255)
                                                                                                                                             aaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGt  <  2:1495716/148‑1 (MQ=255)
                                                                                                                                             aaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGt  <  2:2959554/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGT  >  CP009273/2604587‑2604875

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: