Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 395,846 C→A 10.1% intergenic (+28/+176) fdhA → / ← PP_0329 formaldehyde dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947395,8460CA10.1% 258.0 / 15.7 88intergenic (+28/+176)fdhA/PP_0329formaldehyde dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/40);  new base A (2/7);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

GGTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGC  >  NC_002947/395704‑395989
                                                                                                                                               |                                                                                                                                               
ggTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAg                                                                                                                                           <  6:149081/149‑1 (MQ=255)
  tAGGTGTGCATGTGATCAACCTGTATCAGGCTCAGGAAGGGTATGGCGAGTTTGATGCATGTGTGCCTAAGAAATTCGTTATCGACCCTCACAAAATGTGGTGTGCGGCGTAATCAGCTTCGAGTAAAAAGATCCCCTCGCATCGAggg                                                                                                                                         <  1:327724/149‑1 (MQ=255)
    ggTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGCGAGGGGc                                                                                                                                       <  8:343267/149‑1 (MQ=255)
     gtgtGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGAGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCt                                                                                                                                      >  3:179265/1‑149 (MQ=255)
      tgtgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCtt                                                                                                                                     <  1:158821/149‑1 (MQ=255)
      tgtgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGCGAGGGGCtc                                                                                                                                     <  4:127094/149‑2 (MQ=255)
      tgtgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGCGAGGGGCtc                                                                                                                                     <  4:132461/149‑2 (MQ=255)
       gtgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAAAGCCCCTCGATGCGAgggg                                                                                                                                        <  4:106143/145‑1 (MQ=255)
       gtgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAAAGCCCCTCGATGCGAgggg                                                                                                                                        >  3:106143/1‑145 (MQ=255)
        tgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCtttt                                                                                                                                   <  6:178961/149‑1 (MQ=255)
        tgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAAAGCCCCTCGATGCGAGGGGCtct                                                                                                                                    <  2:199078/148‑3 (MQ=255)
          cAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGCGAGGGGCtcttt                                                                                                                                  <  3:202443/148‑1 (MQ=255)
          cAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCACACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCttttt                                                                                                                                  <  5:126535/148‑1 (MQ=255)
            ggTGATCAACCTGGATCATGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTTCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCTAGTAAAAAGAGCCCCTCTCATCGAGGGGCTTTTTGTc                                                                                                                               <  5:348601/149‑1 (MQ=255)
            ggTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGCGAGGGGCtcttttta                                                                                                                               <  5:28886/149‑4 (MQ=255)
             gTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTc                                                                                                                               >  3:293422/1‑148 (MQ=255)
               gATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAACGAGCCCCTCGCATCGAGGCGCTGTTTGTCCAt                                                                                                                            >  6:260751/1‑149 (MQ=255)
               gATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCAt                                                                                                                            >  4:329889/1‑149 (MQ=255)
                   aaCCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCt                                                                                                                         >  3:206643/1‑148 (MQ=255)
                    aCCTGGATGATGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGTAATTCGTTATCGACCCTCACAAAATGTGGGGTGCGTCGTAATCAGCCTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCGTTTTGTCCATTCtta                                                                                                                       <  6:237019/149‑3 (MQ=255)
                     ccTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGa                                                                                                                       >  7:233822/1‑148 (MQ=255)
                          aTCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGa                                                                                                                  >  6:354004/1‑148 (MQ=255)
                          aTCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGa                                                                                                                  >  2:345769/1‑148 (MQ=255)
                          aTCAGCCGCCGGAAGGGTATGCCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGAAATCGACCCGCACAAAAAGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGAc                                                                                                                 <  5:102582/149‑1 (MQ=255)
                           tCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACa                                                                                                                >  1:151265/1‑149 (MQ=255)
                                  ccGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGc                                                                                                         >  4:318884/1‑149 (MQ=255)
                                       aGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAg                                                                                                    <  7:20121/149‑1 (MQ=255)
                                        gggTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAgg                                                                                                   <  7:280060/149‑1 (MQ=255)
                                        gggTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAgg                                                                                                   <  5:359318/149‑1 (MQ=255)
                                          gTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACAtctc                                                                                                            >  6:310924/1‑138 (MQ=255)
                                          gTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACAtctc                                                                                                            <  5:310924/138‑1 (MQ=255)
                                           tATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGc                                                                                                <  4:206643/149‑1 (MQ=255)
                                               gCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGt                                                                                            >  6:95174/1‑149 (MQ=255)
                                               gCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCACTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGTCCGGAGGTGCTTg                                                                                             >  5:4508/1‑148 (MQ=255)
                                                   gTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATa                                                                                        <  1:345769/149‑1 (MQ=255)
                                                    tttGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCAt                                                                                         >  2:186754/1‑147 (MQ=255)
                                                    tttGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATaa                                                                                       <  5:260751/149‑1 (MQ=255)
                                                     ttGATGCAGGTTTGCCGAAGAAATTCGTTATCGTCCCGCACAAAATTTGGTGTTCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGTGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCTGAGGTGCTTGtaataa                                                                                       <  6:379412/148‑1 (MQ=255)
                                                       gATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTATATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGc                                                                                     >  1:293472/1‑149 (MQ=255)
                                                       gATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGc                                                                                     <  3:318884/148‑1 (MQ=255)
                                                        aTGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGCGGTGCGGCGTAATCACCGTCGAGTAAAAAGAGCCCCTCGCAACGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGCTGCTTGTCATAAGCCt                                                                                   >  2:11469/1‑149 (MQ=255)
                                                         tGCAGGTGTGCCGAAGAAATTTGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCtt                                                                                  >  2:155184/1‑149 (MQ=255)
                                                          gCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCtt                                                                                  >  4:359070/1‑148 (MQ=255)
                                                            aGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCa                                                                               <  5:95174/149‑1 (MQ=255)
                                                                   ccGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTg                                                                                             <  3:248914/128‑1 (MQ=255)
                                                                   ccGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTg                                                                                             >  4:248914/1‑128 (MQ=255)
                                                                    cGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtttt                                                                       >  8:152005/1‑149 (MQ=255)
                                                                    cGAAGAAATTCGTTATCGACCCGCACAAAATGAGGGGTGCGGCGTAATAAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGAAAAGACATCTCGGCCGCAGGTGCTTCTCATAAGCCTTTAATCGCtttt                                                                       >  4:270548/1‑149 (MQ=255)
                                                                    agaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCTGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATAGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtttt                                                                       <  6:53651/148‑1 (MQ=255)
                                                                     gaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTcc                                                                                                                              >  6:250056/1‑93 (MQ=255)
                                                                     gaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTcc                                                                                                                              <  5:250056/93‑1 (MQ=255)
                                                                         aaaTTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCttt                                                                                                                                    <  3:205454/83‑1 (MQ=255)
                                                                         aaaTTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCttt                                                                                                                                    >  4:205454/1‑83 (MQ=255)
                                                                         aaaTTCGTTATCGACCCGCACAAAATGTGGCGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGt                                                                  >  1:32096/1‑149 (MQ=255)
                                                                            ttCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGa                                                               <  5:354004/149‑1 (MQ=255)
                                                                               gTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAgcg                                                            >  2:17583/1‑149 (MQ=255)
                                                                               gTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAgcg                                                            >  3:73183/1‑149 (MQ=255)
                                                                               gTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAgcg                                                            <  4:267895/149‑1 (MQ=255)
                                                                                       cccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGa                                                               <  6:305334/138‑1 (MQ=255)
                                                                                       cccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGATTGTCATAAGCCTTTCATCGCTTTTGAAGTCGa                                                               >  5:305334/1‑138 (MQ=255)
                                                                                        ccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGca                                                   >  5:310154/1‑149 (MQ=255)
                                                                                        ccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGca                                                   <  6:352619/149‑1 (MQ=255)
                                                                                            acaAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATa                                               >  8:81453/1‑149 (MQ=255)
                                                                                              aaaaTGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATAcc                                             <  2:293472/149‑1 (MQ=255)
                                                                                              aaaaTGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATAcc                                             >  6:74900/1‑149 (MQ=255)
                                                                                                aaTGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCtc                                           <  3:270548/149‑1 (MQ=255)
                                                                                                aaTGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCtc                                           <  7:152005/149‑1 (MQ=255)
                                                                                                  tgtgGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCtcat                                         <  8:4966/149‑1 (MQ=255)
                                                                                                    tgGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCtcatc                                        <  8:233822/148‑1 (MQ=255)
                                                                                                    tgGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCtcatc                                        <  2:151265/148‑1 (MQ=255)
                                                                                                        gTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCtt                                                                                              >  1:107866/1‑90 (MQ=255)
                                                                                                        gTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCtt                                                                                              <  2:107866/90‑1 (MQ=255)
                                                                                                            ggcgTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGaa                               <  7:81453/149‑1 (MQ=255)
                                                                                                            ggcgTAATCAGCGGCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGaa                               <  2:312145/149‑1 (MQ=255)
                                                                                                             gcgTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAg                              >  5:19136/1‑149 (MQ=255)
                                                                                                              cgTAATCAGCGTCGAGTAAAACGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGc                             <  6:363788/149‑1 (MQ=255)
                                                                                                                tAATCAGCGGCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCt                           <  3:280627/149‑1 (MQ=255)
                                                                                                                 aaTCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTg                          >  2:172254/1‑149 (MQ=255)
                                                                                                                   aTCAGCGTCGAGTAAAAAGAGCCCATCGCATCTAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGTAGGTGATTGTCATAAGCATTTCATCTCTTTTGAAGTCGAGAGAAGCCAGCACATACCTCATAAGTTCAGAATCCTgt                         >  8:322558/1‑148 (MQ=255)
                                                                                                                    tCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGt                        <  7:322558/149‑1 (MQ=255)
                                                                                                                    tCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGt                        >  5:14631/1‑149 (MQ=255)
                                                                                                                      aGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGa                      <  6:19136/149‑1 (MQ=255)
                                                                                                                            gAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGcc                >  6:249662/1‑149 (MQ=255)
                                                                                                                              ggacaaaaAGCCCCTCGATGCGAGGGGCTCTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCt              >  4:159557/5‑149 (MQ=255)
                                                                                                                              gTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCt              >  1:155978/1‑149 (MQ=255)
                                                                                                                                        ccccTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAgg    >  2:135500/1‑149 (MQ=255)
                                                                                                                                         cccTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAggg   <  4:179265/149‑1 (MQ=255)
                                                                                                                                           cTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGc  >  3:349383/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GGTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGC  >  NC_002947/395704‑395989

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