Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,552,915 C→A 9.6% intergenic (+578/+83) 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,9150CA 9.6% 248.5 / 7.5 99intergenic (+578/+83)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/50);  new base A (9/2);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 2.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACA  >  NC_002947/5552769‑5553061
                                                                                                                                                  |                                                                                                                                                  
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                                                                                                                 <  1:131986/149‑1 (MQ=255)
                                 aTCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgct                                                                                                                 >  5:316436/1‑149 (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                                                                                                             >  8:233025/1‑114 (MQ=255)
                                                                        agCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                             >  8:228745/1‑114 (MQ=255)
                                                                        agCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                             <  7:228745/114‑1 (MQ=255)
                                                                        agCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                             <  7:233025/114‑1 (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                                           <  6:238156/146‑1 (MQ=255)
                                                                                                          ccGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaaa                                           >  5:238156/1‑146 (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                                  <  1:124992/149‑1 (MQ=255)
                                                                                                                cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                  <  7:25383/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                                                         >  7:50831/1‑103 (MQ=255)
                                                                                                                                       gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                         <  8:50831/103‑1 (MQ=255)
                                                                                                                                        ggCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaaa          >  1:210432/1‑149 (MQ=255)
                                                                                                                                             tGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTTCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAAttt       <  1:53450/147‑1 (MQ=255)
                                                                                                                                                ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                  >  3:88696/1‑117 (MQ=255)
                                                                                                                                                ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                  <  4:88696/117‑1 (MQ=255)
                                                                                                                                                ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTa    >  5:27216/1‑147 (MQ=255)
                                                                                                                                                  cGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACa  <  6:259433/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACA  >  NC_002947/5552769‑5553061

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: