Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,330,101 G→T 8.7% H80N (CAT→AAT)  adiA ← arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,330,1010GT 8.7% 735.4 / 29.6 207H80N (CAT→AAT) adiAarginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (81/108);  new base T (9/9);  total (90/117)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 208 total aligned reads displayed.

CAACGGCGCGTCCGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACGGTAACATTT  >  CP009273/4329954‑4330246
                                                                                                                                                   |                                                                                                                                                 
 aaCGGCGCGTCCGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGa                                                                                                                                                  <  2:2668014/148‑1 (MQ=255)
     gcgcGTCCGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCt                                                                                                                                              >  1:2949710/1‑148 (MQ=255)
            cGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGAt                                                                                                                                       <  1:1052336/148‑1 (MQ=255)
            cGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGAt                                                                                                                                       <  1:2133211/148‑1 (MQ=255)
            cGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGAt                                                                                                                                       <  2:101269/148‑1 (MQ=255)
              gCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCa                                                                                                                                     <  1:2355308/148‑1 (MQ=255)
                 aTAAAGTCGGCGGTATCGTCCAGCATCCAGGCGACGGCATCGACAAGCTCCAGCAGGTCGCGATCCATAGCGGCGCGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCAGGATTTGTTGGCGCTCATGATGCTTACCGATCAAtt                                                                                                                                  <  2:1820922/148‑1 (MQ=255)
                  tAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTg                                                                                                                                 >  1:1436789/1‑148 (MQ=255)
                  tAAAGTCGACGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTg                                                                                                                                 >  1:2111660/1‑148 (MQ=255)
                   aaaGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGt                                                                                                                                <  2:2735826/148‑1 (MQ=255)
                    aaGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATTAAGCTTACCGATCAATTGTc                                                                                                                               <  1:903560/148‑1 (MQ=255)
                     aGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGTCCGGCACGTTTTGTTGGCGCTCATTAAGCTTACCGATCAATTGTCt                                                                                                                              >  2:3027207/1‑148 (MQ=255)
                           cggTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGtt                                                                                                                        <  2:411881/148‑1 (MQ=255)
                           cggTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGtt                                                                                                                        <  2:2822385/148‑1 (MQ=255)
                              tATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTg                                                                                                                     <  1:3280229/148‑1 (MQ=255)
                              tATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTg                                                                                                                     <  1:257653/148‑1 (MQ=255)
                               aTCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGa                                                                                                                    <  2:164565/148‑1 (MQ=255)
                               aTCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGa                                                                                                                    <  1:1498984/148‑1 (MQ=255)
                                tCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGAt                                                                                                                   >  1:2806778/1‑148 (MQ=255)
                                  ttCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATTAAGCTTACCGATCAATTGTCTGACGTTTTGATGt                                                                                                                 <  1:2836233/148‑1 (MQ=255)
                                  ttCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGt                                                                                                                 <  1:333039/148‑1 (MQ=255)
                                     cAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCg                                                                                                              <  2:693613/148‑1 (MQ=255)
                                        aaTCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTcc                                                                                                           >  1:1403497/1‑148 (MQ=255)
                                        aaTCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTcc                                                                                                           >  2:917721/1‑148 (MQ=255)
                                           ccAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGa                                                                                                        >  1:616586/1‑148 (MQ=255)
                                            cAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGAt                                                                                                       <  2:1453873/148‑1 (MQ=255)
                                             aGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACCGGAAGACCGGCACGTTTTGTTGGCGCGCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATg                                                                                                      <  2:3063171/148‑1 (MQ=255)
                                             aGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATg                                                                                                      <  1:1308245/148‑1 (MQ=255)
                                             aGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATg                                                                                                      <  1:88104/148‑1 (MQ=255)
                                              ggCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGt                                                                                                     <  1:2296784/148‑1 (MQ=255)
                                               gCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATTAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGtt                                                                                                    >  2:2596741/1‑148 (MQ=255)
                                               gCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGtt                                                                                                    <  2:1703780/148‑1 (MQ=255)
                                               gCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGtt                                                                                                    >  1:1766425/1‑148 (MQ=255)
                                               gCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGtt                                                                                                    <  2:2710314/148‑1 (MQ=255)
                                                cGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTc                                                                                                   >  2:1067530/1‑148 (MQ=255)
                                                  aaTTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATTAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCa                                                                                                 <  2:2286134/148‑1 (MQ=255)
                                                  aaTTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCa                                                                                                 <  2:1970388/148‑1 (MQ=255)
                                                    ttCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCAtt                                                                                               <  2:524230/148‑1 (MQ=255)
                                                     tCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCAttt                                                                                              >  2:719006/1‑148 (MQ=255)
                                                         cGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATa                                                                                          <  2:300698/148‑1 (MQ=255)
                                                         cGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATa                                                                                          >  2:2683503/1‑148 (MQ=255)
                                                         cGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATa                                                                                          <  1:54118/148‑1 (MQ=255)
                                                           aCAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGc                                                                                        >  2:250625/1‑148 (MQ=255)
                                                            cAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCt                                                                                       >  1:2077022/1‑148 (MQ=255)
                                                                cTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATTAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAAc                                                                                   >  2:2660873/1‑148 (MQ=255)
                                                                cTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAAc                                                                                   <  1:1919002/148‑1 (MQ=255)
                                                                  ccAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTGTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACAt                                                                                 >  2:2083128/1‑148 (MQ=255)
                                                                   cagcagGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATc                                                                                <  1:1661122/148‑1 (MQ=255)
                                                                    agcagGTCGCGTTCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCAGGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCa                                                                               <  2:1486985/148‑1 (MQ=255)
                                                                      cagGTCGCGATCCCTTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAgg                                                                             <  1:3295929/148‑1 (MQ=255)
                                                                               aTCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATg                                                                    >  1:602860/1‑148 (MQ=255)
                                                                                tCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATgg                                                                   >  2:1949139/1‑148 (MQ=255)
                                                                                  cATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCt                                                                 <  1:2746243/148‑1 (MQ=255)
                                                                                    ttGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTc                                                               <  1:965150/148‑1 (MQ=255)
                                                                                    ttGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTc                                                               >  1:2368704/1‑148 (MQ=255)
                                                                                    ttGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTc                                                               >  1:2908649/1‑148 (MQ=255)
                                                                                      gcggcgAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGt                                                             <  1:1635751/148‑1 (MQ=255)
                                                                                        ggcgAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGttt                                                           <  1:277078/148‑1 (MQ=255)
                                                                                        ggcgAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGttt                                                           <  2:2106104/148‑1 (MQ=255)
                                                                                         gcgAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTg                                                          >  2:174553/1‑148 (MQ=255)
                                                                                             gggCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAaga                                                      >  2:1669134/1‑148 (MQ=255)
                                                                                             gggCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAaga                                                      >  2:1242631/1‑148 (MQ=255)
                                                                                             gggCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAaga                                                      >  2:284483/1‑148 (MQ=255)
                                                                                                cTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAagagag                                                   >  2:144163/1‑148 (MQ=255)
                                                                                                  ttttCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGaa                                                 <  1:741825/148‑1 (MQ=255)
                                                                                                   tttCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAAt                                                >  1:316356/1‑148 (MQ=255)
                                                                                                     tCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATgg                                              <  2:905904/148‑1 (MQ=255)
                                                                                                      cccGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGc                                             >  2:2308307/1‑148 (MQ=255)
                                                                                                       ccGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCa                                            <  1:1152627/148‑1 (MQ=255)
                                                                                                        cGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCaa                                           >  2:2335055/1‑148 (MQ=255)
                                                                                                         gATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATTAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCaaa                                          <  2:1301866/148‑1 (MQ=255)
                                                                                                          aTCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCaaaa                                         <  2:2908170/148‑1 (MQ=255)
                                                                                                            cGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATgg                                                                   <  2:748902/120‑1 (MQ=255)
                                                                                                            cGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAAcc                                       <  1:2174072/148‑1 (MQ=255)
                                                                                                            agCCCAACAGGTATTCAGGCTCGATTTGTTTGCGCTCATGATGCTTACCGATCAATTGTCTGTCGTTTTGATGTTCGTCAGGATGTTCCATTTGATATCTTTACATCAGGCATTCTTTGGCTTCGTTTGAAGAGAGAATGGCAAAAcc                                       <  1:2174306/147‑1 (MQ=255)
                                                                                                              cccAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACcat                                     >  2:3219370/1‑148 (MQ=255)
                                                                                                               ccAACAGGTAGACCGGCTCGTTTTGTTGGCGCTCATGAATCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTGCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACcatc                                    <  2:484934/148‑1 (MQ=255)
                                                                                                               ccAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATTAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACcatc                                    <  1:2623639/148‑1 (MQ=255)
                                                                                                               ccAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACcatc                                    <  2:2256275/148‑1 (MQ=255)
                                                                                                               ccAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACcatc                                    <  1:1243915/148‑1 (MQ=255)
                                                                                                                 aaCAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACcatcat                                  >  2:2376253/1‑148 (MQ=255)
                                                                                                                 aaCAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACcatcat                                  >  2:280569/1‑148 (MQ=255)
                                                                                                                 aaCAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACcatcat                                  <  2:2028089/148‑1 (MQ=255)
                                                                                                                   cAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCAtcatca                                >  1:226234/1‑148 (MQ=255)
                                                                                                                    aGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCaa                               <  2:3014427/148‑1 (MQ=255)
                                                                                                                     ggAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCaaa                              >  2:948272/1‑148 (MQ=255)
                                                                                                                            cGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTAGTCCGGATGTTCCATTTGATAGCTTAACATCAGGCAGTCAATGGCTTCGTTTTAAGAGAGAATGGCAAAACCATAATCAAaggaggag                       >  2:204867/1‑146 (MQ=255)
                                                                                                                              gCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGa                     <  1:1691839/148‑1 (MQ=255)
                                                                                                                               cACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGAt                    >  1:1119084/1‑148 (MQ=255)
                                                                                                                                   ttttGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTaa                <  2:3281625/148‑1 (MQ=255)
                                                                                                                                    tttGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAAt               >  2:1275838/1‑148 (MQ=255)
                                                                                                                                    tttGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAAt               >  1:394150/1‑148 (MQ=255)
                                                                                                                                       gttgGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCAc            >  2:1533102/1‑148 (MQ=255)
                                                                                                                                         tgGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACgg          <  2:2420076/148‑1 (MQ=255)
                                                                                                                                           gcgcTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACGGTa        >  1:815233/1‑148 (MQ=255)
                                                                                                                                             gcTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACGGTAAc      <  2:737818/148‑1 (MQ=255)
                                                                                                                                               tCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACGGTAACAt    >  1:2118719/1‑148 (MQ=255)
                                                                                                                                                cATTAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACGGTAACAtt   >  2:261414/1‑148 (MQ=255)
                                                                                                                                                cATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACGGTAACAtt   <  2:1161909/148‑1 (MQ=255)
                                                                                                                                                 aTGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACGGTAACAttt  >  2:2157638/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAACGGCGCGTCCGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGCGATCCATTGCGGCGAGGGCTTTTTCCCGATCGCCCAACAGGAAGACCGGCACGTTTTGTTGGCGCTCATGAAGCTTACCGATCAATTGTCTGACGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACGGTAACATTT  >  CP009273/4329954‑4330246

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: ‑