Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,041,297 T→A 6.1% intergenic (+227/‑435) PP_5680 → / → PP_4447 Cro/CI family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,041,2970TA 6.1% 321.4 / 8.7 100intergenic (+227/‑435)PP_5680/PP_4447Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (50/43);  new base A (2/4);  total (53/47)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCATTT  >  NC_002947/5041153‑5041445
                                                                                                                                                |                                                                                                                                                    
tgtTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATgg                                                                                                                                                 <  1:1163280/150‑1 (MQ=255)
    gATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGt                                                                                                                                             >  1:949673/1‑150 (MQ=255)
      tCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCATAAGGCGGtta                                                                                                                                           >  1:1324350/1‑149 (MQ=255)
        cacaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTCAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtct                                                                                                                                         >  2:932986/1‑150 (MQ=255)
         acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtctc                                                                                                                                        >  1:1077718/1‑150 (MQ=255)
         acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtctc                                                                                                                                        >  1:1135846/1‑150 (MQ=255)
         acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCATAAGGCGGttaggg                                                                                                                                        <  1:10214/150‑5 (MQ=255)
         acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCATAAGGCGGttaggg                                                                                                                                        <  1:753362/150‑5 (MQ=255)
          caGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCATAAGGCGGttagggg                                                                                                                                       <  2:1324350/150‑6 (MQ=255)
           agagGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGTCCCCCAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTg                                                                                                                                      >  1:1834034/1‑150 (MQ=255)
            gagGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGc                                                                                                                                     <  2:719535/150‑1 (MQ=255)
            gagGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTATTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGc                                                                                                                                     <  2:1723680/150‑1 (MQ=255)
                gTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTcc                                                                                                                                 <  1:1091497/150‑1 (MQ=255)
                gTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTcc                                                                                                                                 <  2:949673/150‑1 (MQ=255)
                  tAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                                                                               <  1:1212924/150‑1 (MQ=255)
                  tAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                                                                               <  2:1282443/150‑1 (MQ=255)
                        tAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTctgc                                                                                                                         >  1:1324317/1‑150 (MQ=255)
                           ccTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCtgctgc                                                                                                                      >  2:661957/1‑150 (MQ=255)
                             ttAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGTCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCCGCCTCCCAGTCTGCTGCGa                                                                                                                    >  2:146621/1‑150 (MQ=255)
                              tAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGaa                                                                                                                   <  2:1324317/150‑1 (MQ=255)
                                  gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGaa                                                                                                                   >  1:1155287/1‑146 (MQ=255)
                                  gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGaa                                                                                                                   <  2:1155287/146‑1 (MQ=255)
                                    gTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGa                                                                                                             >  2:163126/1‑150 (MQ=255)
                                    gTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGa                                                                                                             >  1:376529/1‑150 (MQ=255)
                                             tGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGt                                                                                                    >  2:1976077/1‑150 (MQ=255)
                                                   cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGcc                                                                                                               <  2:898400/133‑1 (MQ=255)
                                                   cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGcc                                                                                                               >  1:898400/1‑133 (MQ=255)
                                                       acCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAAtt                                                                                          <  2:1249103/150‑1 (MQ=255)
                                                             gCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAttt                                                                                    <  1:661957/150‑1 (MQ=255)
                                                                     gACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGATCAATACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCt                                                                            <  1:636962/150‑1 (MQ=255)
                                                                     gACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCt                                                                            >  1:566853/1‑150 (MQ=255)
                                                                       cATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCACGCCGCAGTACCCTAACCAAGCAACGTCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTgg                                                                          >  1:1498467/1‑150 (MQ=255)
                                                                        aTGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGc                                                                         >  2:1539648/1‑150 (MQ=255)
                                                                         tGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGcc                                                                        <  1:1431047/150‑1 (MQ=255)
                                                                          gACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGaa                                                                                                                   >  1:2113426/1‑106 (MQ=255)
                                                                          gACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGaa                                                                                                                   <  2:2113426/106‑1 (MQ=255)
                                                                              tctTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGAcc                                                                   >  1:2196516/1‑150 (MQ=255)
                                                                               ctTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg                                                                  >  1:617061/1‑150 (MQ=255)
                                                                                 tCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCa                                                                <  2:617061/150‑1 (MQ=255)
                                                                                 tCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCa                                                                <  1:1040790/150‑1 (MQ=255)
                                                                                  aatcccagcgcCTTAAGCCCGCCGACAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCTTCCCAGTCTGCTGCGAAGGCCGTTAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAg                                                               <  1:554028/143‑1 (MQ=255)
                                                                                   aTCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGc                                                                                                                >  1:1380317/1‑100 (MQ=255)
                                                                                   aTCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGc                                                                                                                <  2:1380317/100‑1 (MQ=255)
                                                                                   aTCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGt                                                              >  2:1960202/1‑150 (MQ=255)
                                                                                      ggCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTaga                                                           >  2:2253238/1‑150 (MQ=255)
                                                                                      ggCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCCACTTTTTGACCTGGCCTGACCGCAGTaga                                                           <  1:932986/150‑1 (MQ=255)
                                                                                        cgcgcCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGTGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGc                                                         <  2:338253/150‑1 (MQ=255)
                                                                                        cgcgcCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGc                                                         <  1:1960202/150‑1 (MQ=255)
                                                                                              ttAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGc                                                   <  2:1135846/150‑1 (MQ=255)
                                                                                                 aGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTAt                                                <  1:830689/150‑1 (MQ=255)
                                                                                                   cccGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCa                                              >  1:1885889/1‑150 (MQ=255)
                                                                                                    ccgccgAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCag                                             >  1:2029007/1‑150 (MQ=255)
                                                                                                     cgccgAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGc                                                         <  1:681799/137‑1 (MQ=255)
                                                                                                     cgccgAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGc                                                         >  2:681799/1‑137 (MQ=255)
                                                                                                        cgAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaaccaa                                                                                            <  2:2279139/99‑1 (MQ=255)
                                                                                                        cgAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaaccaa                                                                                            >  1:2279139/1‑99 (MQ=255)
                                                                                                        cgAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGaa                                         >  1:1435112/1‑150 (MQ=255)
                                                                                                               ggCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGt                                  <  1:146621/149‑1 (MQ=255)
                                                                                                               agCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGt                                  <  2:1077718/150‑1 (MQ=255)
                                                                                                                 ccAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTAt                                >  2:219083/1‑150 (MQ=255)
                                                                                                                   agaaaacccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCa                              >  2:532282/7‑150 (MQ=255)
                                                                                                                         cccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTg                                    <  2:2259218/138‑1 (MQ=255)
                                                                                                                         cccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTg                                    >  1:2259218/1‑138 (MQ=255)
                                                                                                                         cccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAatataa                        >  1:847064/1‑148 (MQ=255)
                                                                                                                         cccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCa                        >  2:1747831/1‑150 (MQ=255)
                                                                                                                          ccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGc                                                                                                         <  2:410636/68‑1 (MQ=255)
                                                                                                                          ccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGc                                                                                                         >  1:410636/1‑68 (MQ=255)
                                                                                                                          ccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCTGAGTCCAACCAATTCTa                                                                                       <  2:1231922/86‑1 (MQ=255)
                                                                                                                          ccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTa                                                                                       >  1:1231922/1‑86 (MQ=255)
                                                                                                                          ccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCa                                              <  1:729158/127‑1 (MQ=255)
                                                                                                                          ccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCa                                              >  2:729158/1‑127 (MQ=255)
                                                                                                                          ccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAg                       >  1:2396864/1‑150 (MQ=255)
                                                                                                                          ccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAg                       >  2:122481/1‑150 (MQ=255)
                                                                                                                          ccGAGCAAGACCCCTAACCGCCTTATGGAGGTTTCTCTGCCTCGCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAg                       >  2:25142/1‑150 (MQ=255)
                                                                                                                           cGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCa                                              <  1:35681/126‑1 (MQ=255)
                                                                                                                           cGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCa                                              >  2:35681/1‑126 (MQ=255)
                                                                                                                           cGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGc                      <  2:2196516/150‑1 (MQ=255)
                                                                                                                            gAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCGCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGGGTGCTATCAGAGAAAAGTGGTATCAATATCAGcc                     <  2:123774/150‑1 (MQ=255)
                                                                                                                                  gACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCAc               >  2:595594/1‑150 (MQ=255)
                                                                                                                                     cccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAttttt                                                                                  <  1:373375/80‑1 (MQ=255)
                                                                                                                                     cccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAttttt                                                                                  >  2:373375/1‑80 (MQ=255)
                                                                                                                                     cccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGa            >  1:1993897/1‑150 (MQ=255)
                                                                                                                                      ccgaaCCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTa                                                                                       >  1:2403157/4‑74 (MQ=255)
                                                                                                                                      ccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTa                                                                                       <  2:2403157/74‑1 (MQ=255)
                                                                                                                                      ccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGAc                                                                    >  2:350593/1‑93 (MQ=255)
                                                                                                                                      ccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGAc                                                                    <  1:350593/93‑1 (MQ=255)
                                                                                                                                      agaaaCCGCCATAAGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAat                            <  1:751730/130‑1 (MQ=255)
                                                                                                                                      agaaaCCGCCATAAGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAat                            >  2:751730/4‑133 (MQ=255)
                                                                                                                                       cTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGaaa          >  2:1174408/1‑150 (MQ=255)
                                                                                                                                        tAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGaaaa         >  2:110084/1‑150 (MQ=255)
                                                                                                                                        tAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCAaaaa         <  1:2075597/150‑1 (MQ=255)
                                                                                                                                            cGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGAc                                                                              >  2:1271980/1‑77 (MQ=255)
                                                                                                                                            cGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGAc                                                                              <  1:1271980/77‑1 (MQ=255)
                                                                                                                                            cGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCa     >  2:1596781/1‑150 (MQ=255)
                                                                                                                                            cGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCa     >  2:410258/1‑150 (MQ=255)
                                                                                                                                             gCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCAt    >  2:569952/1‑150 (MQ=255)
                                                                                                                                             gCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCAt    <  2:1498467/150‑1 (MQ=255)
                                                                                                                                              ccTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCAtt   >  2:2326635/1‑150 (MQ=255)
                                                                                                                                              ccTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCAtt   <  2:376529/150‑1 (MQ=255)
                                                                                                                                               cTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGATATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCAttt  <  2:566853/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCATTT  >  NC_002947/5041153‑5041445

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