Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,711,401 T→A 7.7% N244K (AAT→AAA ghrB → glyoxylate/hydroxypyruvate reductase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,711,4010TA 7.7% 1002.7 / 57.7 298N244K (AAT→AAAghrBglyoxylate/hydroxypyruvate reductase B
Reads supporting (aligned to +/- strand):  ref base T (156/119);  new base A (15/8);  total (171/127)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 302 total aligned reads displayed.

GATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCG  >  CP009273/3711252‑3711550
                                                                                                                                                     |                                                                                                                                                     
gATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAAt                                                                                                                                                       >  2:2223606/1‑150 (MQ=255)
gATACTCTGTTACAACAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAAt                                                                                                                                                       >  1:985731/1‑150 (MQ=255)
 aTACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATg                                                                                                                                                      >  1:2401938/1‑150 (MQ=255)
  tACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGc                                                                                                                                                     <  2:722378/150‑1 (MQ=255)
  tACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGc                                                                                                                                                     <  1:831134/150‑1 (MQ=255)
    ctctGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCAc                                                                                                                                                   >  1:24383/1‑150 (MQ=255)
          tACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCg                                                                                                                                             <  2:1080852/150‑1 (MQ=255)
           aCAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCgc                                                                                                                                            <  1:271659/150‑1 (MQ=255)
                 gTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGATACGCTTCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCAGGTGGTTGACGAAAATGCACTGATCGCAGCAtt                                                                                                                                      <  1:1435524/150‑1 (MQ=255)
                 gTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCAtt                                                                                                                                      >  1:3085534/1‑150 (MQ=255)
                 gTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCAtt                                                                                                                                      <  1:3770543/150‑1 (MQ=255)
                 gTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCAtt                                                                                                                                      >  1:905739/1‑150 (MQ=255)
                    aGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGTTGACGAAAATGCCCTGATCGCAGCCTTGCAg                                                                                                                                  >  2:104960/1‑150 (MQ=255)
                      aTTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGa                                                                                                                                 >  2:1786975/1‑150 (MQ=255)
                        ttCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGaaa                                                                                                                               <  1:2022858/150‑1 (MQ=255)
                        ttCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGaaa                                                                                                                               <  2:42913/150‑1 (MQ=255)
                          cGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAgg                                                                                                                             <  2:3583560/150‑1 (MQ=255)
                           gTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGc                                                                                                                            <  1:2177398/150‑1 (MQ=255)
                           gTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAAAGCACTGATCGCAGCATTGCAGAAAGGc                                                                                                                            >  2:2919998/1‑150 (MQ=255)
                            tttGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCAGGTGGTTGACGAAAATGCACTGCTCGCAGCATTGCAGAAAGGCg                                                                                                                           >  1:1664206/1‑150 (MQ=255)
                               gCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGaaa                                                                                                                        <  1:1921286/150‑1 (MQ=255)
                                ccTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAAt                                                                                                                       <  1:3394748/150‑1 (MQ=255)
                                  tGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTc                                                                                                                     <  1:3587545/150‑1 (MQ=255)
                                     tCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACg                                                                                                                  >  1:979568/1‑150 (MQ=255)
                                      ccTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAAAGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGc                                                                                                                 >  1:122504/1‑150 (MQ=255)
                                            gTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCgg                                                                                                           >  2:3171799/1‑150 (MQ=255)
                                              tAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGc                                                                                                         >  2:2253733/1‑150 (MQ=255)
                                              tAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGc                                                                                                         >  2:81525/1‑150 (MQ=255)
                                                aCTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAAAGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTg                                                                                                       >  1:482196/1‑150 (MQ=255)
                                                   gatgaGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGAt                                                                                                    <  2:125787/150‑1 (MQ=255)
                                                    atgaGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATg                                                                                                   <  1:1029653/150‑1 (MQ=255)
                                                    atgaGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATg                                                                                                   >  1:128225/1‑150 (MQ=255)
                                                    atgaGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATg                                                                                                   <  1:3249271/150‑1 (MQ=255)
                                                          cGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCg                                                                                             >  1:2949315/1‑150 (MQ=255)
                                                          cGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCg                                                                                             >  1:1525935/1‑150 (MQ=255)
                                                           gCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGa                                                                                            <  2:3430102/150‑1 (MQ=255)
                                                            catcatCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGaa                                                                                           >  1:1900187/1‑150 (MQ=255)
                                                             atcatcTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAAAGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAAc                                                                                          <  1:4078631/150‑1 (MQ=255)
                                                                atcTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAag                                                                                       >  1:1080354/1‑150 (MQ=255)
                                                                             aGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGt                                                                          >  2:3098770/1‑150 (MQ=255)
                                                                                 cttttCGCCATAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGTAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAt                                                                      <  1:3185143/147‑1 (MQ=255)
                                                                                 cAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAt                                                                      >  1:1865194/1‑150 (MQ=255)
                                                                                 cAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAt                                                                      >  1:1559092/1‑150 (MQ=255)
                                                                                 cAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAt                                                                      <  1:2897283/150‑1 (MQ=255)
                                                                                  aaTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAtt                                                                     >  1:3309491/1‑150 (MQ=255)
                                                                                  aaTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAAAGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAtt                                                                     >  2:3231156/1‑150 (MQ=255)
                                                                                  aaTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAAAGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGAtt                                                                     >  2:1217497/1‑150 (MQ=255)
                                                                                   aTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTc                                                                    <  1:248144/150‑1 (MQ=255)
                                                                                   aTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTc                                                                    >  1:22880/1‑150 (MQ=255)
                                                                                   aTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTc                                                                    <  1:3858956/150‑1 (MQ=255)
                                                                                   aTTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTc                                                                    <  1:274216/150‑1 (MQ=255)
                                                                                        ccAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGt                                                               <  1:2106429/150‑1 (MQ=255)
                                                                                         cAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGtt                                                              >  2:270760/1‑150 (MQ=255)
                                                                                         cAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGtt                                                              >  1:1486183/1‑150 (MQ=255)
                                                                                           aaaTGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGc                                                            >  2:1392137/1‑150 (MQ=255)
                                                                                            aaTGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGct                                                           <  2:1493948/150‑1 (MQ=255)
                                                                                               gAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCtctc                                                        >  1:668104/1‑150 (MQ=255)
                                                                                               gAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCtctc                                                        >  2:1771809/1‑150 (MQ=255)
                                                                                                aaaTCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCa                                                       >  2:3249476/1‑150 (MQ=255)
                                                                                                aaaTCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCa                                                       >  2:4028292/1‑150 (MQ=255)
                                                                                                aaaTCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCa                                                       >  2:2546143/1‑150 (MQ=255)
                                                                                                 aaTCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCaa                                                      >  2:3275943/1‑150 (MQ=255)
                                                                                                  aTCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAAt                                                     >  1:2139429/1‑150 (MQ=255)
                                                                                                    cctccGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATgg                                                   <  1:1638399/150‑1 (MQ=255)
                                                                                                    cctccGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATgg                                                   >  2:2350323/1‑150 (MQ=255)
                                                                                                    cctccGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATgg                                                   <  1:374204/150‑1 (MQ=255)
                                                                                                         gCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAAc                                              >  2:2955250/1‑150 (MQ=255)
                                                                                                         gCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAAc                                              <  2:187512/150‑1 (MQ=255)
                                                                                                          ccATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCaatg                                             <  2:2833954/150‑3 (MQ=255)
                                                                                                          ccATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAAcg                                             <  1:2527024/150‑1 (MQ=255)
                                                                                                           cATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAAcgt                                            >  1:3906055/1‑150 (MQ=255)
                                                                                                            aTTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAAcgtc                                           <  1:2360220/150‑1 (MQ=255)
                                                                                                               ttCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACgtcgtc                                        >  2:4004984/1‑150 (MQ=255)
                                                                                                               ttCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACgtcgtc                                        <  1:733876/150‑1 (MQ=255)
                                                                                                                   ttAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAg                                    >  1:1749099/1‑150 (MQ=255)
                                                                                                                     aaTGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCACTGGCCAACGTCGTCGCAGTa                                  >  2:1994611/1‑150 (MQ=255)
                                                                                                                        gCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCg                               <  2:2120513/150‑1 (MQ=255)
                                                                                                                        gCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCg                               <  2:1624061/150‑1 (MQ=255)
                                                                                                                        gCCGGACGTGGCCCGGTGGTTGACGAAAAAGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCg                               <  1:741020/150‑1 (MQ=255)
                                                                                                                           ggACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCat                            <  2:4337672/150‑1 (MQ=255)
                                                                                                                              cGTGGCCCGGTGGTTGACGAAAAAGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATAtt                         >  2:60567/1‑150 (MQ=255)
                                                                                                                               gTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTg                        >  2:2157267/1‑150 (MQ=255)
                                                                                                                                tGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTgg                       <  1:2815897/150‑1 (MQ=255)
                                                                                                                                 ggCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAACTTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGa                      >  2:2457767/1‑150 (MQ=255)
                                                                                                                                  gCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGAt                     <  1:2505182/150‑1 (MQ=255)
                                                                                                                                  gCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGAt                     <  1:2264026/150‑1 (MQ=255)
                                                                                                                                    ccGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCt                   <  2:697040/150‑1 (MQ=255)
                                                                                                                                     cGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTg                  <  2:3941466/150‑1 (MQ=255)
                                                                                                                                        tggtTGACGAAAAAGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCa               >  1:2640906/1‑150 (MQ=255)
                                                                                                                                          gtTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCAcc             >  1:1153623/1‑150 (MQ=255)
                                                                                                                                            tGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCa           >  2:4170092/1‑150 (MQ=255)
                                                                                                                                            tGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATAGTGGATCTGCCACCCa           >  1:1432475/1‑150 (MQ=255)
                                                                                                                                             gACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCAt          >  2:3644063/1‑150 (MQ=255)
                                                                                                                                                gAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATgag       >  1:3512684/1‑150 (MQ=255)
                                                                                                                                                  aaaTGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAc     >  1:1373744/1‑150 (MQ=255)
                                                                                                                                                  aaaTGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAc     >  1:3642153/1‑150 (MQ=255)
                                                                                                                                                  aaaTGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAc     <  1:516545/150‑1 (MQ=255)
                                                                                                                                                   aaTGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAcg    <  1:3450512/150‑1 (MQ=255)
                                                                                                                                                   aaTGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAcg    >  1:3847389/1‑150 (MQ=255)
                                                                                                                                                     tGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGAcgcg  <  2:1769474/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGGACGTGGCCCGGTGGTTGACGAAAATGCACTGATCGCAGCATTGCAGAAAGGCGAAATTCACGCTGCCGGGCTGGATGTCTTCGAACAAGAGCCACTGTCCGTAGATTCGCCGTTGCTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCG  >  CP009273/3711252‑3711550

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: