Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,552,914 T→A 5.8% intergenic (+577/+84) 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,9140TA 5.8% 202.1 / 2.6 69intergenic (+577/+84)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/36);  new base A (2/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

TACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACA  >  NC_002947/5552769‑5553061
                                                                                                                                                  |                                                                                                                                                     
tACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGc                                                                                                                                                    <  8:324333/149‑1 (MQ=255)
     aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGGTTCGGaaaac                                                                                                                                               >  4:201687/1‑147 (MQ=255)
     aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGc                                                                                                                                               >  7:183874/1‑149 (MQ=255)
         tCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                            <  2:279043/148‑1 (MQ=255)
          cgcgCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCgg                                                                                                                                           <  6:281404/148‑1 (MQ=255)
                 tcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTAGCAAAGCCCGGaaaaatgg                                                                                                                                   >  8:27449/1‑144 (MQ=255)
                 tcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCATTTCTTAGCAAAGCCCGGAAACTggg                                                                                                                                   >  2:240200/1‑149 (MQ=255)
                  ccgatctCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGATCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                   <  1:230727/143‑1 (MQ=255)
                  cATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                  <  5:216567/149‑1 (MQ=255)
                  cATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                  >  4:350194/1‑149 (MQ=255)
                   aTAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGt                                                                                                                                 <  3:244988/149‑1 (MQ=255)
                   aTAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGt                                                                                                                                 >  1:9156/1‑149 (MQ=255)
                       gTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCa                                                                                                                             <  4:170192/149‑1 (MQ=255)
                       gTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCa                                                                                                                             <  2:15187/149‑1 (MQ=255)
                          aCTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTAGCAAAACCCGGAAAAAGGGGTTTAAGca                                                                                                                          >  4:167151/1‑148 (MQ=17)
                          aCTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                   >  2:230727/1‑140 (MQ=255)
                          aCTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCt                                                                                                                          >  7:22963/1‑149 (MQ=255)
                                 aTCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTc                                                                                                                              <  7:46893/138‑1 (MQ=255)
                                 aTCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTc                                                                                                                              >  8:46893/1‑138 (MQ=255)
                                     gTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                               <  2:9156/149‑1 (MQ=255)
                                     gTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                               <  3:201687/149‑1 (MQ=255)
                                        ggTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                             <  5:316311/148‑1 (MQ=255)
                                            cGTCTTGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAAGAAGGTGTTGAGGAATTGCTACTCAGAAGCCCAGACCGGGAGGTCCGGGCTTTGGTTCGCTATGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                       <  5:19685/149‑1 (MQ=255)
                                                 ggTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTa                                                                                                     >  5:338060/1‑146 (MQ=255)
                                                 ggTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTa                                                                                                     <  6:338060/146‑1 (MQ=255)
                                                  gtgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGGTTCGGAAAGCCCGGCAACTGGGGTTTCAGGTGGTGGCTGCTTGTGTGTTTAGaa                                                                                                  >  3:316142/1‑148 (MQ=17)
                                                  gtgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGaa                                                                                                  <  3:19804/148‑1 (MQ=255)
                                                  gtgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTTGTTAGCAAAACCCCGGAACGGGGGGTTAAACAGGGGGGTGCGTGAGTGTTTTGaa                                                                                                  >  4:19804/1‑148 (MQ=17)
                                                   tgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                  >  2:5591/1‑116 (MQ=255)
                                                   tgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                  <  1:5591/116‑1 (MQ=255)
                                                        aGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGGCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATg                                                                                          <  8:310869/149‑1 (MQ=255)
                                                         gCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                         <  3:350194/149‑1 (MQ=255)
                                                            ggTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCATAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGtt                                                                                      <  6:340820/149‑1 (MQ=255)
                                                                 aCCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAAAAATGGGGTGCAtt                                                                                 >  2:35756/1‑149 (MQ=255)
                                                                    aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTg                                                                              <  8:183874/149‑1 (MQ=255)
                                                                    aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTTATTGGCAAAACCCCCAAAAAGCGGTTTAGGGGGGTGGGTGCGTTTGTGTTTTAAAAAAAAATGGtttctataat                                                                               >  8:124942/1‑142 (MQ=9)
                                                                         gCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgtg                                                                                                          <  6:256821/118‑1 (MQ=255)
                                                                         gCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgtg                                                                                                          >  5:256821/1‑118 (MQ=255)
                                                                           aGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCGGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGAc                                                                        >  2:58061/1‑148 (MQ=255)
                                                                                   gCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCATTTCTTAGCAAAGCCCGGAAAATGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTTTAAAAAAATGGGGTGGATTCTGTTTGAACATGGCGTa                                                               >  6:337069/1‑149 (MQ=17)
                                                                                      gAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCt                                                            >  5:346680/1‑149 (MQ=255)
                                                                                        aCATGGTCTTGTGGAATTCCGAACCTGAAGCCAAGAGCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCGGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTc                                                           <  7:124942/148‑1 (MQ=18)
                                                                                         cAAGGTCTTGAGGATTTCCCAACCAGAAGCCAAGACCGTCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                         <  3:14826/149‑1 (MQ=255)
                                                                                         cAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                         <  6:346680/149‑1 (MQ=255)
                                                                                         cAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                         <  3:167151/149‑1 (MQ=255)
                                                                                         cAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                         <  8:54161/149‑1 (MQ=255)
                                                                                          aaGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACACATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                         <  1:58061/149‑1 (MQ=255)
                                                                                              tCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTAAGCTGGTGGGTGCTTGAGTGTTTAGAAATAAATGGGTTGCATTATGTTTGAACAGGGCTTACCTCAAGttc                                                     >  8:87221/1‑146 (MQ=18)
                                                                                                 tGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGGTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTa                                                 >  5:260394/1‑149 (MQ=255)
                                                                                                 tGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGcc                                                                                                                                              <  3:324896/58‑1 (MQ=17)
                                                                                                 tGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGcc                                                                                                                                              >  4:324896/1‑58 (MQ=17)
                                                                                                  gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTAt                                                >  5:287080/1‑149 (MQ=255)
                                                                                                     gAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCa                                             <  1:240200/149‑1 (MQ=255)
                                                                                                     gAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGTCGCAAAGCCCGGAAACTGGGGTTTGACCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTGTGTTTTACCAGGGCTTAACTCAAGTACTatatac                                              >  8:65856/1‑146 (MQ=18)
                                                                                                      aaTTCCGAACCAGAAGCCAAGACCGGCAGGTCCGTTCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaa                                            <  7:27449/149‑1 (MQ=255)
                                                                                                            gAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTGTTAGCCAAGCCCGGAAACTGGGGTTTGACCTGGTGGGTGCTTGAGTGTTTATAAAACAAATGGTTGGATTGTGTTTGACCAGGGCTTACCGGAAGTTCACGATGCAAAGaaa                                       >  4:80325/1‑148 (MQ=17)
                                                                                                            gAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGGTTCGCAAAGCCCGGAAAAGGGGGTTTAAGCTGGTGGCTGGGTGTGTGTTTAGAAAACAATTGGTTGCATTATGTTTGAACAGGGCTTACCTCAAGGACTCTAACCAAAGaaa                                       >  6:84571/1‑148 (MQ=18)
                                                                                                            gAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTAct                                                    <  7:133483/135‑1 (MQ=255)
                                                                                                                cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAg                                                                                                    <  4:64667/84‑1 (MQ=255)
                                                                                                                cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAg                                                                                                    >  3:64667/1‑84 (MQ=255)
                                                                                                                       ccaGACCGTCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAAt                           <  6:270050/147‑1 (MQ=255)
                                                                                                                               gCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGTGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCaa                   <  7:144978/149‑1 (MQ=255)
                                                                                                                                   gTCCGGGCTTTGCTTCGAAAAGCATGGAAACTGGCGTTTCAGCTGGTGTCTGCTTGAGTGTTTAGAAAACAATTGGTTTCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCa                <  3:83260/148‑1 (MQ=255)
                                                                                                                                     tCCGGGCTTTTCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATTGGTTGCATTCTGTTTTACCAGGGCTTACCTCAAGTAATCTATCCAAATAAAAAAAAGAACAATCATTACAATCat               >  4:83260/1‑147 (MQ=25)
                                                                                                                                       cGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                         <  5:138579/104‑1 (MQ=255)
                                                                                                                                       cGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                         >  6:138579/1‑104 (MQ=255)
                                                                                                                                       cGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCa             >  3:128437/1‑148 (MQ=255)
                                                                                                                                            tttGCTTCGCAAAGCCCGTAAACTGGTGTTTCAGCTGGTGGCTTCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAAttt       <  3:143136/149‑1 (MQ=255)
                                                                                                                                                cTTCTCAAAGCCCGGAAACTGGTGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAATCAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTa    <  7:87221/148‑1 (MQ=255)
                                                                                                                                                 ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACa  <  6:287080/149‑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: