Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,611,581 C→A 6.9% intergenic (+162/‑5) deoD → / → yjjJ purine‑nucleoside phosphorylase/putative protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,611,5810CA 6.9% 1248.3 / 43.2 327intergenic (+162/‑5)deoD/yjjJpurine‑nucleoside phosphorylase/putative protein kinase
Reads supporting (aligned to +/- strand):  ref base C (185/119);  new base A (16/7);  total (201/126)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 329 total aligned reads displayed.

CAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGC  >  CP009273/4611433‑4611730
                                                                                                                                                    |                                                                                                                                                     
cAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCa                                                                                                                                                      >  2:314143/1‑150 (MQ=255)
cAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCa                                                                                                                                                      >  2:1383761/1‑150 (MQ=255)
 aaGGCGATTGCCTTGTGAAGCCGGAGCGTGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCAt                                                                                                                                                     <  1:213550/150‑1 (MQ=255)
 aaGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCAt                                                                                                                                                     >  1:3501711/1‑150 (MQ=255)
 aaGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCAt                                                                                                                                                     <  2:3883795/150‑1 (MQ=255)
   ggCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTa                                                                                                                                                   <  2:2284422/150‑1 (MQ=255)
     cGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAAt                                                                                                                                                 >  1:4817005/1‑150 (MQ=255)
         tGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCGATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATgagc                                                                                                                                             <  1:2849939/150‑1 (MQ=255)
          gCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATgagcg                                                                                                                                            <  1:2722878/150‑1 (MQ=255)
                  aaGCCGGAGCGGGAGCCTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGctgact                                                                                                                                    <  2:4831150/150‑1 (MQ=255)
                  aaGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGctgact                                                                                                                                    <  1:1651568/150‑1 (MQ=255)
                    gCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGctgactga                                                                                                                                  <  2:3688388/150‑1 (MQ=255)
                       ggAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCt                                                                                                                               >  2:4001521/1‑150 (MQ=255)
                        gagaGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATATTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAATTATATATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCtt                                                                                                                              >  2:3508116/1‑150 (MQ=255)
                        gAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCtt                                                                                                                              >  1:4654592/1‑150 (MQ=255)
                        gAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCtt                                                                                                                              >  1:3533476/1‑150 (MQ=255)
                        gAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCtt                                                                                                                              >  1:4654353/1‑150 (MQ=255)
                        gAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCtt                                                                                                                              <  1:1462654/150‑1 (MQ=255)
                           cGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTa                                                                                                                           >  1:1232008/1‑150 (MQ=255)
                           cGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTa                                                                                                                           <  2:2029151/150‑1 (MQ=255)
                            gggAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTAc                                                                                                                          >  1:71611/1‑150 (MQ=255)
                            gggAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTAc                                                                                                                          >  1:1006645/1‑150 (MQ=255)
                             ggAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACt                                                                                                                         >  1:4203070/1‑150 (MQ=255)
                               agaCTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGc                                                                                                                       >  2:1014632/1‑150 (MQ=255)
                                 aCTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCaa                                                                                                                     >  1:3883335/1‑150 (MQ=255)
                                  cTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTAATTAATGAGCGAGCTGACTGATCTTTTACTGCAAg                                                                                                                    >  2:450940/1‑150 (MQ=255)
                                   tGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATAAATGAGCGAGCTGACTGATCTTTTACTGCAAgg                                                                                                                   >  1:239567/1‑150 (MQ=255)
                                   tGCTCCGGATTTTTTGTATCTATTCTTTTTTCTCTCCAGCTTGAATTTTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAgg                                                                                                                   <  1:2180061/150‑1 (MQ=255)
                                     cTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCTCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTAAAGTATCTATTCACGAATCTAATAATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGc                                                                                                                 >  2:5146454/1‑150 (MQ=255)
                                      tCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGcc                                                                                                                <  2:3808761/150‑1 (MQ=255)
                                                   tATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCg                                                                                                   <  2:2621499/150‑1 (MQ=255)
                                                    aTCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCgg                                                                                                  >  1:3778388/1‑150 (MQ=255)
                                                     tCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGa                                                                                                 >  2:2228286/1‑150 (MQ=255)
                                                      cTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGaa                                                                                                >  2:4622016/1‑150 (MQ=255)
                                                       tATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTAATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAAt                                                                                               >  2:2961679/1‑150 (MQ=255)
                                                          tCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTgc                                                                                            <  1:2934868/150‑1 (MQ=255)
                                                            aTTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTgcgc                                                                                          >  1:3067736/1‑150 (MQ=255)
                                                             ttgttCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGcc                                                                                         <  1:3390110/147‑1 (MQ=255)
                                                               tttCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAg                                                                                       >  2:2449916/1‑150 (MQ=255)
                                                               tttCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAg                                                                                       >  2:4045929/1‑150 (MQ=255)
                                                               tttCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAg                                                                                       >  1:2299837/1‑150 (MQ=255)
                                                                   tctcCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTc                                                                                   >  1:4930910/1‑150 (MQ=255)
                                                                      ccAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTgg                                                                                >  1:4722830/1‑150 (MQ=255)
                                                                          cTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAAt                                                                            >  1:534182/1‑150 (MQ=255)
                                                                           ttGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAAtc                                                                           <  1:4556543/150‑1 (MQ=255)
                                                                           ttGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAAtc                                                                           >  1:3278071/1‑150 (MQ=255)
                                                                             gAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAAtcag                                                                         >  1:2630333/1‑150 (MQ=255)
                                                                                tatTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAAtcagtca                                                                      >  2:4128109/1‑150 (MQ=255)
                                                                                tatTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAAtcagtca                                                                      >  1:1928054/1‑150 (MQ=255)
                                                                                 atTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAAtcagtcag                                                                     >  2:4552719/1‑150 (MQ=255)
                                                                                    ttCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCg                                                                  >  1:945509/1‑150 (MQ=255)
                                                                                     tCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGa                                                                 >  2:2772/1‑150 (MQ=255)
                                                                                        cTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTAATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGt                                                              <  1:3238329/150‑1 (MQ=255)
                                                                                            aCTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTtctc                                                          >  1:837745/1‑150 (MQ=255)
                                                                                            aCTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTtctc                                                          >  1:4319782/1‑150 (MQ=255)
                                                                                                  cAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCt                                                    <  2:714259/150‑1 (MQ=255)
                                                                                                    gTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCttg                                                  >  1:1731966/1‑150 (MQ=255)
                                                                                                       aaaTGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCttgttg                                               >  1:2600397/1‑150 (MQ=255)
                                                                                                       aaaTGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCttgttg                                               >  1:2241297/1‑150 (MQ=255)
                                                                                                        aaTGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTAATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGc                                              >  2:2868012/1‑150 (MQ=255)
                                                                                                          tgtgTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCa                                            >  2:1332824/1‑150 (MQ=255)
                                                                                                            tgtTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCaga                                          <  2:1633355/150‑1 (MQ=255)
                                                                                                              ttAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAgaga                                        <  1:1461252/150‑1 (MQ=255)
                                                                                                              ttAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTATCACGCCTAGATGCCAgaga                                        >  1:4396795/1‑150 (MQ=255)
                                                                                                               taataaATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGagaa                                       <  2:4166508/150‑1 (MQ=255)
                                                                                                                aataaATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGagaag                                      <  2:3244893/150‑1 (MQ=255)
                                                                                                                      tCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCggg                                >  2:5062808/1‑150 (MQ=255)
                                                                                                                        tATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTg                              <  2:3067736/150‑1 (MQ=255)
                                                                                                                         aTTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGa                             <  2:4890482/150‑1 (MQ=255)
                                                                                                                          ttCAAGTATCTCTTCACGACTCTCTTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGAt                            <  2:1765259/150‑1 (MQ=255)
                                                                                                                           tCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGAtt                           >  2:103506/1‑150 (MQ=255)
                                                                                                                           tCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGAtt                           <  2:804165/150‑1 (MQ=255)
                                                                                                                            cAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTc                          >  2:876520/1‑150 (MQ=255)
                                                                                                                            cAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTc                          >  2:2103325/1‑150 (MQ=255)
                                                                                                                            cAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTc                          >  1:1959977/1‑150 (MQ=255)
                                                                                                                              aGTATCTATTCACGAATCTATTAATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGc                        >  1:2559659/1‑150 (MQ=255)
                                                                                                                                tATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCtt                      >  2:4209089/1‑150 (MQ=255)
                                                                                                                                tATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCtt                      <  1:3367776/150‑1 (MQ=255)
                                                                                                                                 aTCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCttt                     <  2:1323427/150‑1 (MQ=255)
                                                                                                                                   cTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTgg                   >  2:3469803/1‑150 (MQ=255)
                                                                                                                                   cTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTgg                   >  2:4687841/1‑150 (MQ=255)
                                                                                                                                   cTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTgg                   >  2:1530977/1‑150 (MQ=255)
                                                                                                                                   cTATTCACGAATCTATTAATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTgg                   >  2:542632/1‑150 (MQ=255)
                                                                                                                                      ttCACGAATCTATTAATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTaa                <  2:2559659/150‑1 (MQ=255)
                                                                                                                                       tCACGAATCTATTCATTAATGAGAGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGATTAGGTaaa               >  2:1155278/1‑150 (MQ=255)
                                                                                                                                       tCACGAATCTATTAATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTaaa               <  2:4440906/150‑1 (MQ=255)
                                                                                                                                        cACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAg              >  2:1926810/1‑150 (MQ=255)
                                                                                                                                          cGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCa            >  1:3601757/1‑150 (MQ=255)
                                                                                                                                          cGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCa            >  1:4985466/1‑150 (MQ=255)
                                                                                                                                           gAATCTATTAATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCAc           >  1:1532682/1‑150 (MQ=255)
                                                                                                                                            aaTCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACg          >  1:666168/1‑150 (MQ=255)
                                                                                                                                            aaTCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACg          >  2:3573124/1‑150 (MQ=255)
                                                                                                                                             aTCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACgg         <  2:1900687/150‑1 (MQ=255)
                                                                                                                                              tCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACggg        >  2:1260969/1‑150 (MQ=255)
                                                                                                                                                tATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCa      <  1:4140194/150‑1 (MQ=255)
                                                                                                                                                tATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCa      >  2:4017142/1‑150 (MQ=255)
                                                                                                                                                 aTTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCaa     >  2:1290980/1‑150 (MQ=255)
                                                                                                                                                  ttCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAAc    >  1:3837316/1‑150 (MQ=255)
                                                                                                                                                  ttCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAAc    >  2:1234559/1‑150 (MQ=255)
                                                                                                                                                    cATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAAcgc  >  2:4316220/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGC  >  CP009273/4611433‑4611730

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: