Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,222,655 C→G 7.4% intergenic (+111/+262) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,222,6550CG 7.4% 235.3 / 3.4 84intergenic (+111/+262)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base C (3/74);  new base G (5/2);  total (8/76)
Fisher's exact test for biased strand distribution p-value = 3.57e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAT  >  CP009273/2222508‑2222800
                                                                                                                                                    |                                                                                                                                                 
aCTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCa                                                                                                                                                  <  1:1681603/149‑1 (MQ=255)
   ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggg                                                                                                                                               <  2:3941406/149‑1 (MQ=255)
    gtggtgACTTGATTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAgggg                                                                                                                                              <  2:2620343/149‑1 (MQ=255)
      ggtgACTGGCGTATTTAAGATTATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                            <  2:3598770/149‑1 (MQ=255)
      ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                            <  1:3316108/149‑1 (MQ=255)
      ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                            <  1:70156/149‑1 (MQ=255)
      ggtgACTGGCGTATTTAAGATGATATTTTAAAAGTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCgggggg                                                                                                                                            >  1:390898/1‑148 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGgggggggggggc                                                                                                                                        >  2:4358589/1‑147 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcc                                                                                                                                        <  2:3146449/149‑1 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcc                                                                                                                                        <  2:1636317/149‑1 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcc                                                                                                                                        <  2:3964668/149‑1 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcc                                                                                                                                        <  2:3964693/149‑1 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcc                                                                                                                                        <  2:3964718/149‑1 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcc                                                                                                                                        <  1:1345066/149‑1 (MQ=255)
           cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGgggggggggggcc                                                                                                                                       >  1:4754380/1‑146 (MQ=255)
           cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCa                                                                                                                                       <  1:211967/149‑1 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATCTAATCACTGGCGCAAGGAGCGCGGCgggggggggcgc                                                                                                                                      >  1:803867/1‑145 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGgggggggggggggc                                                                                                                                      >  2:896062/1‑145 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGgggggggggggcca                                                                                                                                      >  2:194306/1‑145 (MQ=255)
             ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTGTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATCTTCACTGGCGCAAGGAGCGCGgggggggggggggga                                                                                                                                     >  2:4596176/1‑144 (MQ=255)
             ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATCATCACTGGCGCAAGGAGCGCGgggggggggggccca                                                                                                                                     >  1:665277/1‑144 (MQ=255)
                 tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCAACcgc                                                                                                                                  >  2:2702765/1‑149 (MQ=255)
                   tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgc                                                                                                                               <  2:4430913/149‑1 (MQ=255)
                    ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgcc                                                                                                                              <  1:1872346/149‑1 (MQ=255)
                    ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgcc                                                                                                                              <  1:1372171/149‑1 (MQ=255)
                        gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGccc                                                                                                                          <  1:3403183/149‑1 (MQ=255)
                         atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                         <  2:3616317/149‑1 (MQ=255)
                          tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGccccc                                                                                                                        <  1:4173705/149‑1 (MQ=255)
                          tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCACAAGGAGCGCGCAGGGGGCGGCCAAGCGCCGCCGccccc                                                                                                                        <  1:1406267/149‑1 (MQ=255)
                              atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCCGGGGGGGGCCAAccccccccccccccgcc                                                                                                                     >  2:4222136/1‑146 (MQ=255)
                                  tAAAATTAATTCAGGTCATCAGGTCCGAACATAACGCGAATATTTCAGTCTCTCATCCTGTGGCGCTCCTGTCCTGTGCTGTGCTTCATATACCCACGGGCGCAAGGAGCGCGCCGGGGGCGGCCAACCGCCGCCGCCCCCTGCAcccc                                                                                                                <  1:4839632/149‑1 (MQ=21)
                                       ttaattaaTTTCATAAGGTCCTAAAATAACTATAATATTTAAGTCTCTCATCATGTTGCGCTCCTGTCTTGTGCATTGCTTAATATAATCACTGGCGCAAGTAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGc                                                                                                           <  2:1200603/149‑1 (MQ=25)
                                           ttaatGTCATCAGGTCCTAAAATAACGAGAATATTTAAGTCTATCATCCTGTTGCGGTCCTTTCTTGTGCTTTGCTTCATATAATCACTGGCGCACGGTGCGCGCTGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                       <  1:301295/149‑1 (MQ=21)
                                           ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCAAAcccccccccccccccccccccccccccag                                                                                                        >  1:409203/1‑142 (MQ=255)
                                                     caGGTCCTAAAATAACGATAATTTTTCTGTCTCTCATCCTGTTGCGCTCCGGTCATGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACaaaa                                                                                             <  2:208122/149‑1 (MQ=21)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATACTCACTGGCGCACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAAt                                                                                            <  2:1505845/149‑1 (MQ=255)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGCGGCCAAGCGCCGCCCCCCCCAGCCCCCCCCGGCGCGGGCGAACaaaaa                                                                                            >  2:402194/1‑148 (MQ=255)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAAt                                                                                            <  2:4686228/149‑1 (MQ=255)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAAt                                                                                            <  2:3591351/149‑1 (MQ=255)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAAt                                                                                            <  2:4686305/149‑1 (MQ=255)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAAt                                                                                            <  2:3430464/149‑1 (MQ=255)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCTTGTGCATTGATTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAAt                                                                                            <  2:169156/149‑1 (MQ=255)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGGTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAAt                                                                                            <  2:2759372/149‑1 (MQ=255)
                                                           cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGCCCAACACCCCCCCCCCCCGCCCCCCCCCGCCCCGGGGGACaaaaaacacc                                                                                       >  2:4181881/1‑143 (MQ=255)
                                                               aaTAACTATAATATTTCAGTCTCTCATCGTGTTGCGCTCCTTTCTTGTGCTTTGATTAATATAATCAGGGGCGAACGGAGCGCGCTGGGGGCGCCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgc                                                                                   <  2:3527686/149‑1 (MQ=18)
                                                               aaTAACGATAATATTTCAGTATCTCATCCTGTTGCGCTCCTGTCTTGTGAATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgc                                                                                   <  2:431081/149‑1 (MQ=255)
                                                               aaTAACGAGAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgc                                                                                   <  2:431374/149‑1 (MQ=255)
                                                                          tatTTCTGTCTCTCATCCTGGTGCGCTCCTGTCATGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                        <  2:2459319/149‑1 (MQ=255)
                                                                             ttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGAACGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                     <  2:1458623/149‑1 (MQ=255)
                                                                             ttCAGTCTCTCATCCTGTGGCGCGCCTGTCATGGGCATTGCTTCATATAATCACGGGCGCAAGGTGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                     <  2:4943893/149‑1 (MQ=18)
                                                                               cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGAATTGATCAATATAATCACTGGCGCACGTCGCGCGCCGGGGGCGGCCAACCGCCGCCCCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCtt                                                                   <  2:1722511/149‑1 (MQ=18)
                                                                               cAGTCTCTCACCCGGTTGCGCTCCTTCAATGGGCATTGCTCAATAAAATCACGGGCGCACGTCGCGCGCGGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCtt                                                                   <  2:1722525/149‑1 (MQ=17)
                                                                                    tctcATCGTGTTGCGCTCCTTTATGGTGCTGTGCTTCATATAACCAGTGGCGCAAGGCGCGCGCTGGGGGCGGCCAAGCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATccc                                                              <  1:4550868/149‑1 (MQ=14)
                                                                                    tctcACCCTGTTGCGCTCCTTTAATGTGATTTGCTTAAATTAAGCACGGGCGCATGGTGCGCGCTGGGGGGGGCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATccc                                                              <  1:3222662/149‑1 (MQ=14)
                                                                                             gtggcgcTCCTGTCTGGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCGCAGGGGGCGGCCAATCGCGGCCGCCCCGTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                     <  1:2986471/146‑1 (MQ=17)
                                                                                             gTTGCGCTCTTTTCTTGTGCTTTGTTTAATATAATCCGTGGCGCATGGCGGGCGCGGGGGGGGGCAAAGCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                     <  1:4134538/149‑1 (MQ=14)
                                                                                             gTTGCGCTCTTGTCTTGTGCTTTGATTCATATAATCCGTGGCGCATGTTGCGCGCTGGGGGGGGCAAATCGCGGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                     <  1:4134539/149‑1 (MQ=17)
                                                                                             gTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                     <  2:4206026/149‑1 (MQ=255)
                                                                                              ttGGGCTCGTGTCTTGTGATTTGCTTCATATAATCAGTGGCGCACGGAGCGCGCTGGGGGGGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGTGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                    <  1:5179180/149‑1 (MQ=14)
                                                                                                gcgcTCGTGTCTTTTGCTTGGCTTCATATAATCGCGGGCGCAGGGGGGGGGGGGGGGGCGCCCATGCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCg                                                  <  1:4864280/149‑1 (MQ=14)
                                                                                                gcgcTCCTGTCTTGTGATTTGCTTCATATAAGCACTGGCGAACGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCGCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCg                                                  <  1:1061807/149‑1 (MQ=17)
                                                                                                   cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                               <  2:5016423/149‑1 (MQ=255)
                                                                                                   cTCCTGTCATGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                               <  2:5016976/149‑1 (MQ=255)
                                                                                                    tCCTGTCTTGTGCTTTGCTTCATATAATCAGTGGCGCATGGTGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                              <  1:356328/149‑1 (MQ=17)
                                                                                                    tCCTGTCTTGTGATTTGCTTCATATAATCACTGGCGCACGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                              <  1:4989128/149‑1 (MQ=14)
                                                                                                      cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                            <  2:3370544/149‑1 (MQ=255)
                                                                                                             gtgCTTTGCTTAATATAATCACGGGCGCATGGCGCGCGCGGGGGGGGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                    <  1:4924610/149‑1 (MQ=14)
                                                                                                             gtgCATTGATTCATATAATCACGGGCGCACGGCGCGCGCTGGGGGCGGCCACCCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                    <  1:4924586/149‑1 (MQ=14)
                                                                                                              tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                    <  2:3450058/149‑1 (MQ=255)
                                                                                                               gttttGCTTAATATATTCCCTGGCTAAGGGTGGGGGCTGGGGGGGCCAAACCGCCGCCGCCCCCGCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                   <  1:4101019/146‑1 (MQ=11)
                                                                                                               gctttGATTAATATAATCTGGGGCGCATGGGGCGCGCGGGGGGCGCCACCGCGCCCCCCCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                   <  1:1382130/146‑1 (MQ=11)
                                                                                                                ctttGCTTCATATAACCCCGGGCGCACGGGGCGCGCGGGGGGCGGCCACGCGCCGCCGCCCCCTCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                  <  1:3205802/147‑1 (MQ=14)
                                                                                                                ctttGCTTCATAAAACCCGGGGCGAACGGCGGGCGCGGGGGGGGCCCACCCGCCGCCCCCCCCTCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                  <  1:544812/147‑1 (MQ=11)
                                                                                                                ctttGCTTAATATAATCCGGGGCGCAGGGCGCGCGCTGGGGGCGGCCACGCGCCGCCCCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                  <  1:3205858/147‑1 (MQ=11)
                                                                                                                cagtGCTTCATATAATCACTGGCGCACGGAGCGCGCAGGGGGGGGCCAAGCGCGGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                  <  2:2226180/146‑1 (MQ=14)
                                                                                                                cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                  <  2:1257551/149‑1 (MQ=255)
                                                                                                                cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                  <  2:2429283/149‑1 (MQ=255)
                                                                                                                cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                  <  2:3874073/149‑1 (MQ=255)
                                                                                                                cATTGATTAATATAATCACGGGCGCACGGAGCGCGCTGGGGGCGGCAAACCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                  <  1:3205837/149‑1 (MQ=14)
                                                                                                                   tGCTTCATATAATCACTGGCGCAAGGAGCGCACAGGGGGCGGCCAAACGCCCCCGCCCCCAGCACCCCCCGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGCCCACCGCTTAAACCATACGGCGAACCCGTCGGGCGCGCGGTAACa                               >  1:2797810/1‑149 (MQ=255)
                                                                                                                          tataATCACTGGCGCACGGAGGGCGCGGGGGGCGCCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                        <  1:3255615/149‑1 (MQ=255)
                                                                                                                          tataATCACGGGCGAATGGCGCGCGCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                        <  1:4201943/149‑1 (MQ=11)
                                                                                                                          tataAGCACGGGCGCACGGCGCGCGCGGGGGGCGGCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                        <  1:4201932/149‑1 (MQ=11)
                                                                                                                          aataACCACTGGCCCATGGCGCGCGCGGGGGGCGCCCAAGCGCCGCCGCCCCCGCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                        <  1:1430267/148‑1 (MQ=11)
                                                                                                                                   tGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc               <  2:5198229/149‑1 (MQ=32)
                                                                                                                                    ggCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGcc              <  2:4424567/149‑1 (MQ=17)
                                                                                                                                     gcgcATGGCGCGCGGGGGGGGCGCCCAAGCGCCGCCGCCCCCTCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc             <  1:4369717/149‑1 (MQ=11)
                                                                                                                                     gcgcACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc             <  2:3159425/149‑1 (MQ=17)
                                                                                                                                     gcgcAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc             <  2:2837071/149‑1 (MQ=32)
                                                                                                                                            gAGCGCCCCGGGGGCGCCCACCCGCCGCCGCCCCCGCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc      <  1:4962648/149‑1 (MQ=11)
                                                                                                                                              gggcgcTGGGGGCGGCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac    <  1:1411899/147‑1 (MQ=12)
                                                                                                                                              gcgctcTGGGGGCGGCCTCTGCCGGCCGCCCCATGCACCCCCGGGCTCTGGCGAAAAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac    <  1:1412180/149‑1 (MQ=11)
                                                                                                                                              gcgcgcTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac    <  1:2560463/149‑1 (MQ=17)
                                                                                                                                               cgcgcTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca   <  1:275718/149‑1 (MQ=17)
                                                                                                                                                gcgcGGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt  <  1:4254489/149‑1 (MQ=17)
                                                                                                                                                gcgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt  <  1:3874711/149‑1 (MQ=32)
                                                                                                                                                    |                                                                                                                                                 
ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAT  >  CP009273/2222508‑2222800

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: