Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 58,750 T→A 10.7% intergenic (‑661/+27) PP_5424 ← / ← PP_0050 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294758,7500TA10.7% 219.3 / 15.8 84intergenic (‑661/+27)PP_5424/PP_0050hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base A (5/4);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-01

AGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTCCA  >  NC_002947/58604‑58896
                                                                                                                                                    |                                                                                                                                                   
aGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGctgctg                                                                                                                                                    <  1:154747/148‑1 (MQ=255)
   aGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCgg                                                                                                                                                 <  3:227658/148‑1 (MQ=255)
    gggATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGc                                                                                                                                               >  2:318179/1‑149 (MQ=255)
     gtttttCACGGATTTTATGTAAAATGAACCCATTAATTTAAAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGGCACTAAAAAGCCCGCAGCAGCGGGCtt                                                                                                                                             <  5:192262/146‑1 (MQ=255)
     ggATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCt                                                                                                                                              >  2:224585/1‑149 (MQ=255)
      gATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGGCACTAAAAAGCCCGCAGCAGCGGGCtt                                                                                                                                             <  2:187396/149‑1 (MQ=255)
        tttCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCtttt                                                                                                                                           <  2:271233/149‑1 (MQ=255)
           ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAg                                                                                                                                        >  5:345505/1‑149 (MQ=255)
           ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                         <  3:161120/148‑1 (MQ=255)
            tcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCtttt                                                                                                                                           <  7:278309/145‑1 (MQ=255)
            tcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTCTAATTCGCCAGACACTAAAAAGACCGCTGCTGCGGCCtttt                                                                                                                                           >  8:278309/1‑145 (MQ=255)
             gtgCTTTTCACAGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTg                                                                                                                                      >  1:133426/3‑149 (MQ=255)
              ggCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                         >  4:205412/1‑145 (MQ=255)
              ggCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATCCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                         <  3:205412/145‑1 (MQ=255)
                      acagAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCaaa                                                                                                                              >  5:333877/4‑148 (MQ=255)
                        tgAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACaa                                                                                                                           >  3:76043/1‑149 (MQ=255)
                         gAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAg                                                                                                                                        <  5:346371/135‑1 (MQ=255)
                         gAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAg                                                                                                                                        >  6:346371/1‑135 (MQ=255)
                          aaaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCggg                                                                                                                                                <  6:162535/126‑1 (MQ=255)
                          aaaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCggg                                                                                                                                                >  5:162535/1‑126 (MQ=255)
                          aaaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAAt                                                                                                                          <  1:318179/148‑1 (MQ=255)
                            aaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAg                                                                                                                       <  7:116935/149‑1 (MQ=255)
                              gCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCTCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTa                                                                                                                     <  1:60221/149‑1 (MQ=255)
                                   ccTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCgg                                                                                                                >  8:280738/1‑149 (MQ=255)
                                    cTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTTGTGCCTGCAAACAATCAGTAGCCGGt                                                                                                               <  1:9598/149‑1 (MQ=255)
                                    cTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGt                                                                                                               <  6:345505/149‑1 (MQ=255)
                                          ttCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAAt                                                                                                         >  5:340733/1‑149 (MQ=255)
                                            cacaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACTA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGt                                                                                                       >  6:340842/1‑149 (MQ=255)
                                            cacaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGt                                                                                                       >  5:73752/1‑149 (MQ=255)
                                             acaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGt                                                                                                       >  5:321386/1‑148 (MQ=255)
                                               aaaaaaaaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGGGAGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAAAAATCAGTAGCCGGTGAGAATGt                                                                                                       >  1:280068/1‑148 (MQ=255)
                                               aaaaaaaaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCa                                                                                                    >  5:220285/1‑149 (MQ=255)
                                               aaaaaaaaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCa                                                                                                    <  4:56570/149‑1 (MQ=255)
                                                aaaaaaaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCAc                                                                                                   <  7:280738/149‑1 (MQ=255)
                                                 aaaaaaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATAAGTAGCCGGTGAGAATGATCAc                                                                                                   >  6:241142/1‑148 (MQ=255)
                                                         gCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACAcc                                                                                          >  8:91956/1‑149 (MQ=255)
                                                            aaGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                         <  7:150926/99‑1 (MQ=255)
                                                            aaGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                         >  8:150926/1‑99 (MQ=255)
                                                             aGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTg                                                                                      <  6:321386/149‑1 (MQ=255)
                                                              gATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTgg                                                                                     <  6:340733/149‑1 (MQ=255)
                                                                 tCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                         <  6:296716/94‑1 (MQ=255)
                                                                   cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGTCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGTCAATg                                                                                <  5:241142/149‑1 (MQ=255)
                                                                   cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAAt                                                                                 <  2:133426/148‑1 (MQ=255)
                                                                   cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATg                                                                                >  1:69884/1‑149 (MQ=255)
                                                                     aTTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACTA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGc                                                                              <  5:340842/149‑1 (MQ=255)
                                                                     aTTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGc                                                                              <  2:142941/149‑1 (MQ=255)
                                                                         cACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGccac                                                                           <  2:239673/148‑1 (MQ=255)
                                                                             gggTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGaa                                                                      <  6:220285/149‑1 (MQ=255)
                                                                                   ggCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGccac                                                                           >  1:331915/1‑138 (MQ=255)
                                                                                   ggCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGccac                                                                           <  2:331915/138‑1 (MQ=255)
                                                                                       tATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGc                                                            >  2:345786/1‑148 (MQ=255)
                                                                                          cACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTa                                                                                                                     >  8:233981/1‑89 (MQ=255)
                                                                                          cACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTa                                                                                                                     <  7:233981/89‑1 (MQ=255)
                                                                                           aCCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTc                                                       <  6:333877/149‑1 (MQ=255)
                                                                                             cTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATAAGTAGCCGTTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                                      >  6:116620/1‑148 (MQ=255)
                                                                                                    gtgAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAg                                               <  4:76043/148‑1 (MQ=255)
                                                                                                      gAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGccac                                                                           <  6:239502/119‑1 (MQ=255)
                                                                                                      gAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGccac                                                                           >  5:239502/1‑119 (MQ=255)
                                                                                                       aaCCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTga                                                                                                             <  8:71881/84‑1 (MQ=255)
                                                                                                       aaCCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTga                                                                                                             >  7:71881/1‑84 (MQ=255)
                                                                                                         ccTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAAc                                          >  1:307123/1‑148 (MQ=255)
                                                                                                         ccTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGAGCCTGCAAACAATCAGTAGCCGGTGACAATGTTCACGATCACACACCTGGCAATGGCCACCAGAACAAAGTCGCCATTCAAATACAGCCAAc                                          >  5:144234/1‑148 (MQ=255)
                                                                                                                  ggTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCaa                                           >  8:199578/1‑140 (MQ=255)
                                                                                                                  ggTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCaa                                           <  7:199578/140‑1 (MQ=255)
                                                                                                                    tgtAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCAGCCATTCACATACAGCCAACGATGATCa                                  <  1:21317/148‑1 (MQ=255)
                                                                                                                    tgtAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGAATTTTAGTGCCTCCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCCCTCTCAATGGCCACCAGAACATACTCGCCATTCACATACAGCCAACGATGATAACg                                >  8:42429/1‑149 (MQ=255)
                                                                                                                      tAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGa                               <  7:91956/148‑1 (MQ=255)
                                                                                                                       aaTTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAgg                             >  1:119891/1‑149 (MQ=255)
                                                                                                                            gCCCGACCCTAACACCCCCGCGGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGgc                         <  2:280068/148‑1 (MQ=255)
                                                                                                                            gCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGgcg                        <  2:69884/149‑1 (MQ=255)
                                                                                                                              cAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCa                      >  6:53450/1‑149 (MQ=255)
                                                                                                                                 gcacTAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGaa                                                                      >  6:67669/2‑99 (MQ=255)
                                                                                                                                 gcacTAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGaa                                                                      <  5:67669/98‑1 (MQ=255)
                                                                                                                                   acTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAAGGTTCACGATCACACCGATGGCAATGGCCACCAGAACATAGTCGGCATTCACATACAGCCAACGATGATCACGAGGTGGAGCATACAg                 >  3:18251/1‑149 (MQ=255)
                                                                                                                                    cTAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAgg                >  6:323705/1‑149 (MQ=255)
                                                                                                                                     tAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGt               >  4:160262/1‑149 (MQ=255)
                                                                                                                                     tAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGt               >  3:334995/1‑149 (MQ=255)
                                                                                                                                     tAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGTCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAAAATAGTCGCCAGTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGt               >  1:247227/1‑149 (MQ=39)
                                                                                                                                      aaaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGGATACAGGTg              >  7:277938/1‑149 (MQ=255)
                                                                                                                                       aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAg                 >  7:197712/1‑145 (MQ=255)
                                                                                                                                       aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAg                 <  8:197712/145‑1 (MQ=255)
                                                                                                                                       aaaaGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTgg             >  6:49115/1‑149 (MQ=255)
                                                                                                                                       aaaaGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTgg             >  3:20443/1‑149 (MQ=255)
                                                                                                                                          aGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGgcgc          <  4:334995/149‑1 (MQ=255)
                                                                                                                                            cccGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                                      >  6:70191/1‑103 (MQ=255)
                                                                                                                                                   ctgcGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGTTCACGAGGTGGCGCATACAGGTGGCGCGCTTTCca  <  1:224585/148‑1 (MQ=255)
                                                                                                                                                    tgcGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTcc   >  1:338134/1‑146 (MQ=255)
                                                                                                                                                    tgcGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTcc   <  2:338134/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTCCA  >  NC_002947/58604‑58896

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