Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,688,062 C→A 9.3% D281Y (GAT→TAT)  purL ← phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,688,0620CA 9.3% 744.3 / 37.0 203D281Y (GAT→TAT) purLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (58/126);  new base A (2/17);  total (60/143)
Fisher's exact test for biased strand distribution p-value = 6.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 205 total aligned reads displayed.

GCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTT  >  CP009273/2687914‑2688211
                                                                                                                                                    |                                                                                                                                                     
gCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCg                                                                                                                                                      <  2:453273/150‑1 (MQ=255)
 cTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGt                                                                                                                                                     <  1:3835193/150‑1 (MQ=255)
     cGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCg                                                                                                                                                 >  1:2642539/1‑150 (MQ=255)
     cGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCg                                                                                                                                                 >  1:2642520/1‑150 (MQ=255)
          ttGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCg                                                                                                                                            >  1:673873/1‑150 (MQ=255)
          ttGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCg                                                                                                                                            <  2:3695930/150‑1 (MQ=255)
          ttGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATAGTAGCGGCCCg                                                                                                                                            >  2:2198520/1‑150 (MQ=255)
           tGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGt                                                                                                                                           >  2:2115472/1‑150 (MQ=255)
            gcgccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGtt                                                                                                                                          <  2:2485046/150‑1 (MQ=255)
             cgccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATAGTAGCGGCCCGttt                                                                                                                                         <  2:609850/150‑1 (MQ=255)
              gccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTc                                                                                                                                        <  1:2548011/150‑1 (MQ=255)
               ccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCg                                                                                                                                       <  1:2987869/150‑1 (MQ=255)
                cgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGt                                                                                                                                      <  2:1644128/150‑1 (MQ=255)
                  cgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTATCATCAGAATATGCGCCGGTTCCTGATGGACATCGTAGCGGCCCGTTTCGTgg                                                                                                                                    >  1:3079611/1‑150 (MQ=255)
                     ccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCa                                                                                                                                 <  1:849800/150‑1 (MQ=255)
                                    tCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGccc                                                                                                                  <  2:119602/150‑1 (MQ=255)
                                    tCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGccc                                                                                                                  <  2:1223001/150‑1 (MQ=255)
                                       cgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACt                                                                                                               <  1:2784463/150‑1 (MQ=255)
                                        ggATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACtt                                                                                                              >  1:1855943/1‑150 (MQ=255)
                                        ggATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATAGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACtt                                                                                                              <  2:2561816/150‑1 (MQ=255)
                                        ggATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATAGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACtt                                                                                                              <  2:2894649/150‑1 (MQ=255)
                                          aTTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCa                                                                                                            <  1:1093726/150‑1 (MQ=255)
                                          aTTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCa                                                                                                            >  2:1838301/1‑150 (MQ=255)
                                               aCCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAAcc                                                                                                       <  2:2826880/150‑1 (MQ=255)
                                                    cgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCa                                                                                                  >  2:1524998/1‑150 (MQ=255)
                                                          cGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcgg                                                                                            >  2:560267/1‑150 (MQ=255)
                                                            gTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcggcg                                                                                          <  2:3346563/150‑1 (MQ=255)
                                                             tCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcggcgg                                                                                         <  1:1395131/150‑1 (MQ=255)
                                                              cgccgcGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACggcggc                                                                                        <  2:3444930/150‑1 (MQ=255)
                                                               gccgcGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGgcggcg                                                                                       <  2:2034379/150‑1 (MQ=255)
                                                               gccgcGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGgcggcg                                                                                       <  1:3699830/150‑1 (MQ=255)
                                                               gccgcGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGgcggcg                                                                                       <  2:2933358/150‑1 (MQ=255)
                                                                ccgcGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGt                                                                                      <  2:4019898/150‑1 (MQ=255)
                                                                ccgcGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGt                                                                                      <  1:2408224/150‑1 (MQ=255)
                                                                 cgcgCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGtt                                                                                     >  2:786797/1‑150 (MQ=255)
                                                                   cgcCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTAt                                                                                   >  1:2064886/1‑150 (MQ=255)
                                                                   cgcCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTAt                                                                                   <  2:344079/150‑1 (MQ=255)
                                                                     cccGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCt                                                                                 <  1:1197868/150‑1 (MQ=255)
                                                                         gCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATAGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTtat                                                                             <  2:2870878/150‑1 (MQ=255)
                                                                           cACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATaa                                                                           >  2:3120152/1‑150 (MQ=255)
                                                                           cACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATaa                                                                           <  2:2367985/150‑1 (MQ=255)
                                                                               gagaAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCag                                                                       >  1:2379202/1‑150 (MQ=255)
                                                                               gagaAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCag                                                                       >  1:559297/1‑150 (MQ=255)
                                                                                   aaTCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagagag                                                                   <  1:4055390/150‑1 (MQ=255)
                                                                                   aaTCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagagag                                                                   <  2:2003776/150‑1 (MQ=255)
                                                                                      cGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAAc                                                                >  2:3066399/1‑150 (MQ=255)
                                                                                          gTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGACCGTGa                                                            >  2:1690517/1‑150 (MQ=255)
                                                                                          gTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGa                                                            >  2:2553901/1‑150 (MQ=255)
                                                                                            cGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATc                                                          <  2:1712019/150‑1 (MQ=255)
                                                                                            cGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATc                                                          >  2:2981036/1‑150 (MQ=255)
                                                                                              ggtggtTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATAGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTg                                                        >  2:3919019/1‑150 (MQ=255)
                                                                                                tggtTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGc                                                      <  2:2908742/150‑1 (MQ=255)
                                                                                                 ggtTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCg                                                     <  2:3722624/150‑1 (MQ=255)
                                                                                                  gtTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGt                                                    <  1:819769/150‑1 (MQ=255)
                                                                                                  gtTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGt                                                    <  1:2282724/150‑1 (MQ=255)
                                                                                                     gtgAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGt                                                 >  2:754069/1‑150 (MQ=255)
                                                                                                      tgAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATAGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGtt                                                <  1:471342/150‑1 (MQ=255)
                                                                                                       gAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGttt                                               <  2:149033/150‑1 (MQ=255)
                                                                                                        aGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTc                                              >  2:397178/1‑150 (MQ=255)
                                                                                                         gTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCg                                             >  2:321715/1‑150 (MQ=255)
                                                                                                         gTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCg                                             <  1:877943/150‑1 (MQ=255)
                                                                                                          tttCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGa                                            <  2:2146269/150‑1 (MQ=255)
                                                                                                             cGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAg                                         <  1:3601457/150‑1 (MQ=255)
                                                                                                              gACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATAGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGt                                        <  2:3606202/150‑1 (MQ=255)
                                                                                                              gACTTTCATCAGAAGATGCGCCGGTTCCTGATGGAAATAGTAGCGGCCCGTTTCGTGGTCAGCAACGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGt                                        <  2:3605554/150‑1 (MQ=255)
                                                                                                                cTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGACCCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTAt                                      <  1:943993/150‑1 (MQ=255)
                                                                                                                cTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTAt                                      >  1:499533/1‑150 (MQ=255)
                                                                                                                cTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCATGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTAt                                      <  1:2182345/150‑1 (MQ=255)
                                                                                                                 tttCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTAtt                                     <  1:3933861/150‑1 (MQ=255)
                                                                                                                  ttCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTAttt                                    <  1:3845059/150‑1 (MQ=255)
                                                                                                                   tcatcaGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTAtttt                                   <  2:281951/150‑1 (MQ=255)
                                                                                                                     atcaGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTg                                 >  1:164265/1‑150 (MQ=255)
                                                                                                                      tcaGAATATGCGCCGGTTCCTGATGGAAATAGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGa                                <  1:1239111/150‑1 (MQ=255)
                                                                                                                        aGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatc                              >  1:3589449/1‑150 (MQ=255)
                                                                                                                         gAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatca                             <  1:2115472/150‑1 (MQ=255)
                                                                                                                          tatatGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCCGAACCTTCCATTACGGCGGCGTTATCTTTATTAGCAGAGTGAACGTGATCTTGCGTGGTTTCGAAAGTATTTTTGatcat                            <  2:2321696/149‑1 (MQ=255)
                                                                                                                           atatGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcatc                           <  2:860810/150‑1 (MQ=255)
                                                                                                                           atatGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcatc                           >  1:3472207/1‑150 (MQ=255)
                                                                                                                           atatGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGAAAAGTAGCGGCCCACTTCAGAACCTTCAATTACGGCGGCGTTATCTTTATAAGCAGAGAGAAAGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcatc                           >  1:666903/1‑150 (MQ=255)
                                                                                                                             atGCGCCGGGTCCTGTTGGAACTCGTTGCGGCCCGTTTCGTGGTCTGCTAAGTTGCGGCCCACTTCAGACCCTTCCATTACGGCGGCGTTATCTTTATAAGCTGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGagcatcgt                         <  2:7303/150‑3 (MQ=255)
                                                                                                                              tGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTg                        <  1:2904819/150‑1 (MQ=255)
                                                                                                                                cgcCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGaa                      <  1:1439010/150‑1 (MQ=255)
                                                                                                                                cgcCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGaa                      >  2:3650715/1‑150 (MQ=255)
                                                                                                                                    ggTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGc                  >  2:2444614/1‑150 (MQ=255)
                                                                                                                                        ccTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGAtt              <  1:3412818/150‑1 (MQ=255)
                                                                                                                                        ccTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGAtt              <  1:2243338/150‑1 (MQ=255)
                                                                                                                                         cTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGAttt             <  1:1039906/150‑1 (MQ=255)
                                                                                                                                         cTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGAttt             <  2:1050865/150‑1 (MQ=255)
                                                                                                                                         cTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGAttt             <  2:2727251/150‑1 (MQ=255)
                                                                                                                                          tGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTc            <  1:2167496/150‑1 (MQ=255)
                                                                                                                                          tGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTc            <  1:1700042/150‑1 (MQ=255)
                                                                                                                                          tGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTGATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTc            >  1:356503/1‑150 (MQ=255)
                                                                                                                                           gATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCg           <  1:2287254/150‑1 (MQ=255)
                                                                                                                                             tGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGgc         <  2:2288625/150‑1 (MQ=255)
                                                                                                                                              ggAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGgct        <  2:3869652/150‑1 (MQ=255)
                                                                                                                                              ggAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGgct        >  1:2214916/1‑150 (MQ=255)
                                                                                                                                                 aaTCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGgctgct     <  2:281404/150‑1 (MQ=255)
                                                                                                                                                 aaTCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGgctgct     >  2:3010652/1‑150 (MQ=255)
                                                                                                                                                  aTCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGctgctg    <  1:2730574/150‑1 (MQ=255)
                                                                                                                                                    cGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGtt  <  1:3909385/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTT  >  CP009273/2687914‑2688211

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: