Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,825,763 C→A 53.1% A159A (GCG→GCT astC ← succinylornithine transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,825,7630CA53.1% 5.2 / inf 275A159A (GCG→GCTastCsuccinylornithine transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (70/59);  new base A (89/57);  total (159/116)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 280 total aligned reads displayed.

AGACCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTT  >  CP009273/1825614‑1825911
                                                                                                                                                     |                                                                                                                                                    
  aCCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGgagc                                                                                                                                                    <  1:534028/150‑1 (MQ=255)
   ccTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCa                                                                                                                                                   >  1:1517337/1‑150 (MQ=255)
    cTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCaa                                                                                                                                                  >  1:283475/1‑150 (MQ=255)
     ttGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCaaa                                                                                                                                                 <  2:4734376/150‑1 (MQ=255)
     ttGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCaaa                                                                                                                                                 <  2:2466139/150‑1 (MQ=255)
       gTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAAt                                                                                                                                               >  2:3523151/1‑150 (MQ=255)
       gTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTACCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAAt                                                                                                                                               >  1:3102222/1‑150 (MQ=255)
        tAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATc                                                                                                                                              <  2:3896860/150‑1 (MQ=255)
         aaaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATcc                                                                                                                                             >  2:2302435/1‑150 (MQ=255)
         aaaaaCGCGGTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATcc                                                                                                                                             <  1:3499814/150‑1 (MQ=255)
                cgTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGACTATCCGCCGGCAGTGGCGCAAAATCCTGTGAAt                                                                                                                                      >  2:4452550/1‑150 (MQ=255)
                  ttGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAg                                                                                                                                    >  1:3533112/1‑150 (MQ=255)
                       ggctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAggctgg                                                                                                                               <  2:1596503/150‑1 (MQ=255)
                           ggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGc                                                                                                                           >  1:2288527/1‑150 (MQ=255)
                           ggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGc                                                                                                                           >  1:98303/1‑150 (MQ=255)
                             caccacAGCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGccc                                                                                                                         <  1:4083502/150‑1 (MQ=255)
                              accacACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTCGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCa                                                                                                                        >  2:2185776/1‑150 (MQ=255)
                              accacACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCa                                                                                                                        >  2:2011235/1‑150 (MQ=255)
                              accacACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCa                                                                                                                        >  1:918295/1‑150 (MQ=255)
                               ccacACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCAc                                                                                                                       >  1:1929605/1‑150 (MQ=255)
                                 acacCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGACCAccc                                                                                                                     >  2:606658/1‑150 (MQ=255)
                                    ccGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCa                                                                                                                  <  1:4178434/150‑1 (MQ=255)
                                      gCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACt                                                                                                                >  2:80819/1‑150 (MQ=255)
                                       ccTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATAAACTTATATTAAGCATGACGAATATACGACGGCAGTGGAGCAAAATACAGTGAATAGGCAGGCTGCCCACCCGCAATg                                                                                                               >  1:3287931/1‑150 (MQ=255)
                                           cccccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAg                                                                                                           <  2:1808056/150‑1 (MQ=255)
                                            ccccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGt                                                                                                          <  2:1970942/150‑1 (MQ=255)
                                             cccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTa                                                                                                         <  2:1249821/150‑1 (MQ=255)
                                              ccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTaa                                                                                                        <  2:2782559/150‑1 (MQ=255)
                                                tGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAAc                                                                                                      <  2:4059056/150‑1 (MQ=255)
                                                 ggATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACa                                                                                                     <  2:18985/150‑1 (MQ=255)
                                                  gATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAg                                                                                                    >  1:1016670/1‑150 (MQ=255)
                                                   aTGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGc                                                                                                   <  2:2627957/150‑1 (MQ=255)
                                                   aTGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGc                                                                                                   <  1:4576796/150‑1 (MQ=255)
                                                    tGGGTTCGACAATCACCTCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCg                                                                                                  <  2:499261/150‑1 (MQ=255)
                                                    tGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCg                                                                                                  <  1:4097290/150‑1 (MQ=255)
                                                     gggTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGt                                                                                                 <  2:3926909/150‑1 (MQ=255)
                                                       gTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGc                                                                                               <  2:2105395/150‑1 (MQ=255)
                                                         tCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGa                                                                                             >  2:2949833/1‑150 (MQ=255)
                                                           gACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGAcc                                                                                           >  1:1682976/1‑150 (MQ=255)
                                                           gACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGAcc                                                                                           <  1:55284/150‑1 (MQ=255)
                                                            aCAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCa                                                                                          >  1:3373189/1‑150 (MQ=255)
                                                             cAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCAt                                                                                         <  2:4213323/150‑1 (MQ=255)
                                                              aaTCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATg                                                                                        >  1:4026855/1‑150 (MQ=255)
                                                              aaTCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATg                                                                                        >  1:4372125/1‑150 (MQ=255)
                                                                 cACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGaaa                                                                                     >  1:1338321/1‑150 (MQ=255)
                                                                 cACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGaaa                                                                                     <  1:1606657/150‑1 (MQ=255)
                                                                 cACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGaaa                                                                                     >  1:2419394/1‑150 (MQ=255)
                                                                   ccGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACg                                                                                   >  2:3600863/1‑150 (MQ=255)
                                                                      cacaGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAt                                                                                >  1:2192656/1‑150 (MQ=255)
                                                                       acaGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtt                                                                               >  2:3762615/1‑150 (MQ=255)
                                                                       acaGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtt                                                                               <  2:3830071/150‑1 (MQ=255)
                                                                        caGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAttt                                                                              >  2:2557644/1‑150 (MQ=255)
                                                                         aGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtttt                                                                             <  2:4692862/150‑1 (MQ=255)
                                                                         aGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtttt                                                                             <  1:4416799/150‑1 (MQ=255)
                                                                         aGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCTAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtttt                                                                             <  1:205167/150‑1 (MQ=255)
                                                                         aGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtttt                                                                             >  2:2679144/1‑150 (MQ=255)
                                                                         aGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtttt                                                                             <  1:3937916/150‑1 (MQ=255)
                                                                          ggTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAttttt                                                                            <  2:864776/150‑1 (MQ=255)
                                                                            tAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGa                                                                          <  1:3828822/150‑1 (MQ=255)
                                                                            tAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGa                                                                          >  2:1069543/1‑150 (MQ=255)
                                                                             agagTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGaa                                                                         >  2:570147/1‑150 (MQ=255)
                                                                                  cgtcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCa                                                                    >  1:1816678/1‑150 (MQ=255)
                                                                                    tcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACg                                                                  <  1:2157351/150‑1 (MQ=255)
                                                                                      aaTCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGAt                                                                >  1:4133060/1‑150 (MQ=255)
                                                                                      aaTCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGAt                                                                <  2:2293590/150‑1 (MQ=255)
                                                                                       aTCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATg                                                               >  1:1078774/1‑150 (MQ=255)
                                                                                       aTCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATg                                                               >  2:3453233/1‑150 (MQ=255)
                                                                                          aGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCACACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCg                                                            <  1:2204770/150‑1 (MQ=255)
                                                                                                 tGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTAt                                                     >  2:4326741/1‑150 (MQ=255)
                                                                                                 tGGCAGAGTTAATATCGTTATATGCAGCATGACGAAGATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTAt                                                     <  1:729622/150‑1 (MQ=255)
                                                                                                  ggCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATg                                                    >  1:3417218/1‑150 (MQ=255)
                                                                                                        gTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGcc                                              >  2:1143509/1‑150 (MQ=255)
                                                                                                         ttAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCg                                             >  1:1247000/1‑150 (MQ=255)
                                                                                                            atatCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAg                                          >  1:1891/1‑150 (MQ=255)
                                                                                                            atatCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAg                                          >  2:4459004/1‑150 (MQ=255)
                                                                                                             tatCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGc                                         <  1:3741195/150‑1 (MQ=255)
                                                                                                                   tataTGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGt                                   <  2:941434/150‑1 (MQ=255)
                                                                                                                    atatGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTg                                  <  1:2926562/150‑1 (MQ=255)
                                                                                                                    atatGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTg                                  >  2:1864196/1‑150 (MQ=255)
                                                                                                                     tatGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGa                                 <  1:2521858/150‑1 (MQ=255)
                                                                                                                           gcaTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAt                           <  1:309578/150‑1 (MQ=255)
                                                                                                                           gcaTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAt                           <  2:2377066/150‑1 (MQ=255)
                                                                                                                            caTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAtt                          >  2:1572489/1‑150 (MQ=255)
                                                                                                                             aTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAttt                         >  2:3967489/1‑150 (MQ=255)
                                                                                                                             aTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAttt                         >  2:4589046/1‑150 (MQ=255)
                                                                                                                             aTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAttt                         >  2:1507477/1‑150 (MQ=255)
                                                                                                                                aCGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcg                      <  2:4363143/150‑1 (MQ=255)
                                                                                                                                aCGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcg                      <  2:3108419/150‑1 (MQ=255)
                                                                                                                                 cGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcgc                     >  1:4264718/1‑150 (MQ=255)
                                                                                                                                  gAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcgcg                    >  2:4536729/1‑150 (MQ=255)
                                                                                                                                   aaTATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACgcgc                   <  2:1417272/150‑1 (MQ=255)
                                                                                                                                          gccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGtttt            <  2:1335648/150‑1 (MQ=255)
                                                                                                                                           ccgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTa           <  1:1648451/150‑1 (MQ=255)
                                                                                                                                           ccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTa           <  2:1560434/150‑1 (MQ=255)
                                                                                                                                            cgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAg          >  1:1736092/1‑150 (MQ=255)
                                                                                                                                             ggCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgc         >  1:2567342/1‑150 (MQ=255)
                                                                                                                                             ggCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgc         >  1:2995877/1‑150 (MQ=255)
                                                                                                                                               cAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgcgc       <  2:2419394/150‑1 (MQ=255)
                                                                                                                                               cAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgcgc       >  1:3752970/1‑150 (MQ=255)
                                                                                                                                                aGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGcgcc      >  2:669157/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGACCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTT  >  CP009273/1825614‑1825911

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: