Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,552,913 T→G 5.2% intergenic (+576/+85) PP_4882 → / ← PP_5717 iron ABC transporter permease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,552,9130TG 5.2% 252.7 / 6.1 96intergenic (+576/+85)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/49);  new base G (3/2);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

TACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTA  >  NC_002947/5552769‑5553059
                                                                                                                                                |                                                                                                                                                  
tACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGc                                                                                                                                                <  1:114924/149‑1 (MQ=255)
 aCTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGc                                                                                                                                                >  6:57998/1‑148 (MQ=255)
   tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCaaa                                                                                                                                             >  4:62779/1‑149 (MQ=255)
     aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCaaa                                                                                                                                             <  8:145323/147‑1 (MQ=255)
     aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGc                                                                                                                                           >  7:238093/1‑149 (MQ=255)
     aTTATCGAGCTCTCATAAGTCACTGCGCATCAGATGGTACGTCTGGTGTACAGCCGGTTCACCAACAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCTAACCAGAAGCCAAGAACTCCAGCTCCGGCCTTAGCATCGAACAGc                                                                                                                                           >  6:88448/1‑149 (MQ=18)
       tttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGccc                                                                                                                                         >  5:81670/1‑149 (MQ=255)
        ttCGGGCTCTCATAAGTAACGGCGCATCAGTGGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGccc                                                                                                                                         <  6:222247/148‑1 (MQ=18)
         tCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                        <  3:296493/148‑1 (MQ=255)
               tctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTg                                                                                                                                 >  3:79034/1‑149 (MQ=255)
                 tcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGGTTTGTTAGCGAACCCCGGCAACAggg                                                                                                                               >  2:131986/1‑149 (MQ=17)
                       gTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCATTGCTTAGCAAAGCCCGGAAACTGGGGattaa                                                                                                                         >  3:197817/1‑147 (MQ=255)
                            tGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCATCGCAAAGCCCGGAACCTGGGGTTTACGCtgg                                                                                                                    >  6:91292/1‑149 (MQ=255)
                                 aTCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgct                                                                                                               >  5:316436/1‑149 (MQ=255)
                                 aTCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgct                                                                                                               <  1:131986/149‑1 (MQ=255)
                                  tCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgctg                                                                                                              >  1:46221/1‑149 (MQ=255)
                                    aGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgctgc                                                                                                             <  6:150051/148‑1 (MQ=255)
                                     gTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                           >  7:134778/1‑149 (MQ=255)
                                      ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGGTTCGCAAAGCCCGGAAACTGGGGTTTACGGGGGTGGCTGGttt                                                                                                          >  8:227988/1‑148 (MQ=17)
                                        ggTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTCGCAAACCCCGGAAACTGGGGTTTGCGCTGGTGGCTGCTTgtg                                                                                                        >  6:287529/1‑147 (MQ=255)
                                        ggTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                         <  3:285997/148‑1 (MQ=255)
                                        ggTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                         <  7:232868/148‑1 (MQ=255)
                                            cGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                    <  3:5496/149‑1 (MQ=255)
                                            cGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                    <  5:287529/149‑1 (MQ=255)
                                                  gtgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGa                                                                                                <  6:81670/147‑1 (MQ=255)
                                                  gtgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCGGGTGGCTGCATGAGTGTTTTGaa                                                                                               >  6:127025/1‑148 (MQ=255)
                                                   tgtACAGCCGGTTCACGAAGTAGCAGGATTATTCTGAACAAGGTCTTGAGGAATTCCGAACCTGAATCCAAGAACGGCAGGTCCGGGCTTTGCTTCGCAAAGACCGTAAACTGGTGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGaaaa                                                                                             <  5:88448/149‑1 (MQ=17)
                                                    gtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCCGAAACTGGGGTTTCAGCGGGTGGCGGCTTGAGTGTTTAGaaaa                                                                                             >  5:259433/1‑148 (MQ=255)
                                                     tACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAAc                                                                                            <  5:263979/148‑1 (MQ=255)
                                                         gCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                       <  7:211943/149‑1 (MQ=255)
                                                          ccGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATggg                                                                                      <  3:144637/149‑1 (MQ=255)
                                                           gggTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgctgc                                                                                                             <  4:30587/124‑1 (MQ=255)
                                                            ggTTCACCAAGAAGCAGGCTTATGCTGATCAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTGTTCGCAAAGCCCGGTAACTTGGGTTTCAGCTGGTGGCTGCTTTAGTGTTTTGTAAACAAAggggt                                                                                    >  3:295095/1‑147 (MQ=17)
                                                            ggTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGtt                                                                                    >  7:270644/1‑149 (MQ=255)
                                                                cACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCa                                                                                 <  5:127025/148‑1 (MQ=255)
                                                                   cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                               <  5:259014/99‑1 (MQ=255)
                                                                   cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGTATAAAATGGGTTGGATTcg                                                                             >  6:241973/1‑148 (MQ=255)
                                                                   cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCt                                                                             >  6:126036/1‑149 (MQ=255)
                                                                   cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCCCAAACACCCGAAACATGGGTTTCAGCAGGTGCATGCATGAGTGTTTATAAAACAAATGGTTTCaatat                                                                             >  3:184247/1‑145 (MQ=17)
                                                                   cAAGAAGCAGCCTTATGATGAACAAGGTCTTGAGGAATTAAGAACAAGAAGCCAAGACCGGCATGTACGAGCTATGTTTCGCAAAGCCCGGAAACTGGGGTTTCTGCTGGTGGCTGCTTGAGTGTTTATGAAACAAAGGGTTGCAttgt                                                                             >  4:160409/1‑147 (MQ=14)
                                                                     agaagCATGCTTATTCTGAACAAGGGCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGGTTAGAAAACAATGGGTTGCATTCTGt                                                                           <  8:134778/149‑1 (MQ=255)
                                                                     agaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                           >  5:147986/1‑149 (MQ=255)
                                                                      gaagCAGGCTTATGCTGAACAAGGGCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGtt                                                                          <  3:62779/149‑1 (MQ=255)
                                                                        agCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                           <  7:228745/114‑1 (MQ=255)
                                                                        agCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                           <  7:233025/114‑1 (MQ=255)
                                                                        agCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                           >  8:228745/1‑114 (MQ=255)
                                                                        agCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                           >  8:233025/1‑114 (MQ=255)
                                                                         gCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGGTTCGCAAAGCCCGGAAACTGGGGTTTAAGCTGGTGGCGGGTTGAGTGTTTTGAAAACAATGGGTTGCATTATGTTTGa                                                                       >  8:63456/1‑149 (MQ=14)
                                                                           aGGCTTATGCTGAACAAGGTCTTGAGTAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCATAGCAAAACCCCGAAAAAGGGGTTTAAGCAGGTGGCTGCTTGAGTGTTTATAAAACAATGGGTTTCAATATGttttca                                                                      >  5:287383/1‑145 (MQ=17)
                                                                           aGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGAc                                                                      <  5:241973/148‑1 (MQ=255)
                                                                               ttATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAggg                                                                 <  5:57998/149‑1 (MQ=255)
                                                                                 aTGCTTAACAAGGTCTTTTGGAATTCCTAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGAAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAgg                                                                  <  4:16542/146‑1 (MQ=17)
                                                                                 aTGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGc                                                                <  4:79034/148‑1 (MQ=255)
                                                                                  tGCTGACCAAGGTCTTGAGGAATTCCTAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCtt                                                              <  7:227988/149‑1 (MQ=255)
                                                                                       aaCAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCt                                                          <  6:316436/148‑1 (MQ=255)
                                                                                       aaCAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTc                                                         <  7:63456/149‑1 (MQ=255)
                                                                                       aaCAAGGGCTTGTGGAATACCGAACCAGAAGCCAAGACCTGCAGGTCCGGGCTTTGCTTCGCAACGCCCGGAAACTGGGGTTTCAGCTGTTGGCTGCTTGATTGTTTAGAAAACAATGTGTTGCATTCTGTTTGACCAGGGCTTACCTc                                                         <  6:137516/149‑1 (MQ=18)
                                                                                        aCAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTc                                                         <  8:270644/148‑1 (MQ=255)
                                                                                         cAAGGTCTTGATGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGtt                                                                                    <  5:65139/120‑1 (MQ=255)
                                                                                         cAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                       >  7:243249/1‑149 (MQ=255)
                                                                                           aGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTCCTTCGCAAAGCCCGGAAACTgggg                                                                                                                              <  7:43505/76‑1 (MQ=255)
                                                                                              tCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTAct                                                  <  2:46221/149‑1 (MQ=255)
                                                                                              tCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGAACAGGGCTTACCTCAAGTAct                                                  >  2:53450/1‑149 (MQ=255)
                                                                                                 tGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGa                                                                       <  3:231588/125‑1 (MQ=255)
                                                                                                 tGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGa                                                                       >  4:231588/1‑125 (MQ=255)
                                                                                                  gacgAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAgg                                                                  <  5:284833/126‑1 (MQ=255)
                                                                                                  gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTCGCAAAGCCCCGAAACAGGGGTTTAAGCTGGTGGCTGCTTGAGTGTTTTTAAAAAAATGGGTTGCATTATGTTTTACCAgg                                                                  >  6:284833/1‑129 (MQ=17)
                                                                                                  gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCATTTCTTAGCAAACCCCGGCAACTGGGGTTTAAGCTGGTGGCTGCGTGAGTGTTTAGAAAAAAAAGGGTTGCATTCTTTTTGACCAGGGCTTACCTCAAGTAcaaga                                               >  8:245103/1‑144 (MQ=18)
                                                                                                         tCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAggg                                                                 <  5:15286/123‑1 (MQ=255)
                                                                                                          ccGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaaa                                         >  5:238156/1‑146 (MQ=255)
                                                                                                          ccGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaaa                                         <  6:238156/146‑1 (MQ=255)
                                                                                                           cGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTGTTCGCAAAGCCCGGAAACTGGGGTTTCAGCGGGTGGCTGCGTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATACAAAGaa                                      >  7:116971/1‑148 (MQ=21)
                                                                                                           cGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTGAGCTGGGGGCTGCGTGAGTGTTTAGAAAACAATGGGGTGGATTCTGTTTGGCCAGGGCTTACCTCAAGTACTCTATCCAAAGaa                                      >  2:87160/1‑148 (MQ=21)
                                                                                                             aaCCAGAAGCCAAGACCTGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaa                                   <  8:238093/149‑1 (MQ=255)
                                                                                                               ccAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCGATCCAAAGaaaaaa                                  >  3:42800/1‑148 (MQ=255)
                                                                                                                cAGAAGCCAAGCCCGGCATGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAACGAATGGGGTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                <  4:197817/149‑1 (MQ=21)
                                                                                                                cAGAAGCCAAGACCGGCAGGTCCGGGCTTTTGTTAGCAAAGCCCCGAAACTGGGGTTTGAGCTGGTGGCGGCGTGAGTGTTTAGAAAAAAAAGGGTTTCATTGTGTTTTACCAGGGCGTAGCGGAATTACTGTATCCAAAGaaaaaaa                                 >  4:101725/1‑148 (MQ=17)
                                                                                                                cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                <  7:25383/149‑1 (MQ=255)
                                                                                                                cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                <  1:124992/149‑1 (MQ=255)
                                                                                                                 aGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGa                              <  3:32205/149‑1 (MQ=255)
                                                                                                                   aaGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGaa                             <  5:126036/149‑1 (MQ=255)
                                                                                                                    aGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTCCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGaa                             <  7:245103/148‑1 (MQ=255)
                                                                                                                       cAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAAt                         <  8:243249/149‑1 (MQ=255)
                                                                                                                          gACCGGCAGGTCCGGGCTTTGCTTAGCAAAGCCCGGAAACTGGGGTTTCCTCTGGTGGCTGCGTGATTGTTTAGAAAACCATGGGTTGGAGTATGTTTTACCAGGGCTTAACTCAAGTACTCTATCCAAAGAAAAAAAAGAACAAtacg                      >  6:240673/1‑146 (MQ=21)
                                                                                                                          gACCGGCAGGTCCGGGATTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAg                      <  4:295095/149‑1 (MQ=255)
                                                                                                                             cGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTg                    >  7:239539/1‑148 (MQ=255)
                                                                                                                                agtttCCGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGaa                             >  5:5110/5‑136 (MQ=255)
                                                                                                                                   gTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAg             <  3:145723/149‑1 (MQ=255)
                                                                                                                                      cGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                       >  7:56083/1‑104 (MQ=255)
                                                                                                                                      cGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                       <  8:56083/104‑1 (MQ=255)
                                                                                                                                      cGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCa           >  6:79181/1‑148 (MQ=255)
                                                                                                                                      cGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaa                                          >  8:23022/1‑117 (MQ=255)
                                                                                                                                       gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                       <  8:50831/103‑1 (MQ=255)
                                                                                                                                       gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                       >  7:50831/1‑103 (MQ=255)
                                                                                                                                        ggCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaaa        >  1:210432/1‑149 (MQ=255)
                                                                                                                                             tGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTTCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAAttt     <  1:53450/147‑1 (MQ=255)
                                                                                                                                                ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                <  4:88696/117‑1 (MQ=255)
                                                                                                                                                ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                >  3:88696/1‑117 (MQ=255)
                                                                                                                                                ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTa  >  5:27216/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTA  >  NC_002947/5552769‑5553059

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: