Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,711,448 T→A 9.9% L260Q (CTG→CAG)  ghrB → glyoxylate/hydroxypyruvate reductase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,711,4480TA 9.9% 558.8 / 42.8 192L260Q (CTG→CAG) ghrBglyoxylate/hydroxypyruvate reductase B
Reads supporting (aligned to +/- strand):  ref base T (84/89);  new base A (12/7);  total (96/96)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 196 total aligned reads displayed.

TGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTT  >  CP009273/3711302‑3711595
                                                                                                                                                  |                                                                                                                                                   
tgatgaGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTg                                                                                                                                                    >  1:739078/1‑148 (MQ=255)
 gatgaGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTgg                                                                                                                                                   <  1:3276031/148‑1 (MQ=255)
 gatgaGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTgg                                                                                                                                                   >  2:2665985/1‑148 (MQ=255)
     agaCGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGt                                                                                                                                               <  1:22567/148‑1 (MQ=255)
      gaCGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTc                                                                                                                                              >  2:947093/1‑148 (MQ=255)
       aCGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCt                                                                                                                                             >  1:1489054/1‑148 (MQ=255)
       aCGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCt                                                                                                                                             <  2:3090717/148‑1 (MQ=255)
       aCGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCt                                                                                                                                             >  1:1166882/1‑148 (MQ=255)
           atcatcTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGa                                                                                                                                         >  1:3329949/1‑148 (MQ=255)
           atcatcTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGa                                                                                                                                         <  2:2529414/148‑1 (MQ=255)
                cTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAag                                                                                                                                    <  2:739078/148‑1 (MQ=255)
                cTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAag                                                                                                                                    >  1:3494150/1‑148 (MQ=255)
                 tGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAaga                                                                                                                                   <  2:2159391/148‑1 (MQ=255)
                        cgcAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTg                                                                                                                            <  2:1722741/148‑1 (MQ=255)
                        cgcAGAACAATTCGCCAAAATGAAATCCTCCGCCATGTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTg                                                                                                                            <  2:400869/148‑1 (MQ=255)
                         gcAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCAGGATGTCTTCGAACAAGAGCCACTGt                                                                                                                           <  2:2628853/148‑1 (MQ=255)
                          cAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTc                                                                                                                          <  1:1671257/148‑1 (MQ=255)
                          cAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTc                                                                                                                          >  2:1677481/1‑148 (MQ=255)
                               cAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAg                                                                                                                     >  1:327415/1‑148 (MQ=255)
                                aaTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGa                                                                                                                    <  1:659681/148‑1 (MQ=255)
                                 aTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACTCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAt                                                                                                                   >  1:2775097/1‑148 (MQ=255)
                                 aTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAt                                                                                                                   <  1:422866/148‑1 (MQ=255)
                                 aTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAt                                                                                                                   <  2:482260/148‑1 (MQ=255)
                                 aTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCAGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAt                                                                                                                   <  2:1547154/148‑1 (MQ=255)
                                  ttCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCAGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAtt                                                                                                                  >  1:539683/1‑148 (MQ=255)
                                    cGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCg                                                                                                                <  1:3289917/148‑1 (MQ=255)
                                     gCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGc                                                                                                               >  1:2753616/1‑148 (MQ=255)
                                      ccAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGcc                                                                                                              >  2:1914293/1‑148 (MQ=255)
                                      ccAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGcc                                                                                                              >  1:217466/1‑148 (MQ=255)
                                       cAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCg                                                                                                             <  2:775014/148‑1 (MQ=255)
                                       cAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCg                                                                                                             >  2:1784235/1‑148 (MQ=255)
                                        aaaaTGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGt                                                                                                            <  1:3520824/148‑1 (MQ=255)
                                         aaaTGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCCGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACACGAGCCACTGTCCGTAGATTCGCCGtt                                                                                                           >  2:3354688/1‑148 (MQ=255)
                                          aaTGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTg                                                                                                          <  2:2809461/148‑1 (MQ=255)
                                              aaaTCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGctct                                                                                                      <  1:654909/148‑1 (MQ=255)
                                              aaaTCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCAGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGctct                                                                                                      >  1:953141/1‑148 (MQ=255)
                                               aaTCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCtctc                                                                                                     >  2:1278015/1‑148 (MQ=255)
                                                aTCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCa                                                                                                    <  2:1940380/148‑1 (MQ=255)
                                                  cctccGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAAt                                                                                                  >  2:1029215/1‑148 (MQ=255)
                                                  cctccGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCAGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAAt                                                                                                  >  2:3409043/1‑148 (MQ=255)
                                                     ccGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTAGCTCTCAATGGc                                                                                               <  1:1411565/148‑1 (MQ=255)
                                                      cGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGcc                                                                                              >  2:1376885/1‑148 (MQ=255)
                                                      cGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGcc                                                                                              >  2:1376416/1‑148 (MQ=255)
                                                       gCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCa                                                                                             >  1:3521642/1‑148 (MQ=255)
                                                        ccATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCaa                                                                                            >  1:658373/1‑148 (MQ=255)
                                                         cATTTTCATTAATTCCGGACGTGGACCGGTGGTAGACGAAAATTCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTAGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCaaa                                                                                           >  1:2083107/1‑147 (MQ=255)
                                                          aTTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAAcg                                                                                          <  2:531699/148‑1 (MQ=255)
                                                           ttttCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAAcgt                                                                                         <  1:1029749/148‑1 (MQ=255)
                                                           ttttCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAAcgt                                                                                         >  2:1713580/1‑148 (MQ=255)
                                                            tttCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAAcgtc                                                                                        >  1:516406/1‑148 (MQ=255)
                                                            tttCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAAcgtc                                                                                        >  1:1238365/1‑148 (MQ=255)
                                                               cATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACgtcgtc                                                                                     >  1:3536660/1‑148 (MQ=255)
                                                                 ttAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGc                                                                                   <  2:3200864/148‑1 (MQ=255)
                                                                     tGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTa                                                                               >  1:3063730/1‑148 (MQ=255)
                                                                     tGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTa                                                                               >  1:3320390/1‑148 (MQ=255)
                                                                     tGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTa                                                                               >  1:3063454/1‑148 (MQ=255)
                                                                      gCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTAc                                                                              >  2:222188/1‑148 (MQ=255)
                                                                      gCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTAc                                                                              >  1:1024516/1‑148 (MQ=255)
                                                                         ggACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGc                                                                           <  2:1907803/148‑1 (MQ=255)
                                                                           aCGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCat                                                                         <  2:1999322/148‑1 (MQ=255)
                                                                            cGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCata                                                                        >  1:890422/1‑148 (MQ=255)
                                                                                 cccGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGa                                                                   <  1:1784235/148‑1 (MQ=255)
                                                                                     gtggtTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTg                                                               <  1:1914293/148‑1 (MQ=255)
                                                                                     gtggtTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTg                                                               <  2:2841066/148‑1 (MQ=255)
                                                                                     gtggtTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTg                                                               <  2:2098787/148‑1 (MQ=255)
                                                                                         ttGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCAc                                                           >  2:3132928/1‑148 (MQ=255)
                                                                                                  aTGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAc                                                  >  2:3538245/1‑148 (MQ=255)
                                                                                                  aTGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAc                                                  >  2:3072899/1‑148 (MQ=255)
                                                                                                   tGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAcg                                                 >  1:880766/1‑148 (MQ=255)
                                                                                                    gCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAcgc                                                >  1:1091265/1‑148 (MQ=255)
                                                                                                     cACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAcgcg                                               <  1:2886423/148‑1 (MQ=255)
                                                                                                      aCTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGt                                              <  1:2510705/148‑1 (MQ=255)
                                                                                                        tGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTa                                            >  1:2842311/1‑148 (MQ=255)
                                                                                                             gcagcaTTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTatggca                                       <  2:1168171/148‑1 (MQ=255)
                                                                                                                   ttGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGgccg                                 >  2:3012676/1‑148 (MQ=255)
                                                                                                                         aaaGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCtgtg                           >  1:2343135/1‑148 (MQ=255)
                                                                                                                         aaaGGCGAAAGTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCtgtg                           <  2:2190146/148‑1 (MQ=255)
                                                                                                                         aaaGGCGAAAGTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCtgtg                           <  1:1576357/148‑1 (MQ=255)
                                                                                                                           aGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGcc                         >  2:1756438/1‑148 (MQ=255)
                                                                                                                           aGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGcc                         >  2:1060084/1‑148 (MQ=255)
                                                                                                                            ggCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCg                        >  1:1298239/1‑148 (MQ=255)
                                                                                                                            ggCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCg                        >  2:1356051/1‑148 (MQ=255)
                                                                                                                               gAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTgg                     <  1:589008/148‑1 (MQ=255)
                                                                                                                               gAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTgg                     <  2:611002/148‑1 (MQ=255)
                                                                                                                               gAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTgg                     <  1:2889802/148‑1 (MQ=255)
                                                                                                                                  aTTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATa                  <  1:2595068/148‑1 (MQ=255)
                                                                                                                                   ttCACGCTGCCGGGCAGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATaa                 >  2:2939105/1‑148 (MQ=255)
                                                                                                                                   ttCACGCTGCCGGGCAGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATaa                 >  2:1162041/1‑148 (MQ=255)
                                                                                                                                    tCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAAt                <  2:2550481/148‑1 (MQ=255)
                                                                                                                                    tCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAAt                >  2:870264/1‑148 (MQ=255)
                                                                                                                                    tCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGATGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATaag                >  1:284731/1‑147 (MQ=255)
                                                                                                                                     cACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAAtt               <  1:827993/148‑1 (MQ=255)
                                                                                                                                      aCGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAAttt              <  2:1238365/148‑1 (MQ=255)
                                                                                                                                       cGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGACCGCCTGTGCCGTGGATAATttg             >  1:3181453/1‑148 (MQ=255)
                                                                                                                                              gggCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATg      <  2:1489054/148‑1 (MQ=255)
                                                                                                                                              gggCAGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATg      <  1:3408388/148‑1 (MQ=255)
                                                                                                                                                gCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCg    <  2:890422/148‑1 (MQ=255)
                                                                                                                                                 cTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGt   <  2:647722/148‑1 (MQ=255)
                                                                                                                                                  tGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGtt  >  2:2486153/1‑148 (MQ=255)
                                                                                                                                                  tGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGtt  <  1:1043447/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTT  >  CP009273/3711302‑3711595

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: