Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 425,348 A→C 8.7% intergenic (‑39/+18) PP_0349 ← / ← PP_0350 membrane protein/ferrichrome‑iron receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947425,3480AC 8.7% 180.1 / 5.2 64intergenic (‑39/+18)PP_0349/PP_0350membrane protein/ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base A (27/31);  new base C (3/3);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTC  >  NC_002947/425202‑425494
                                                                                                                                                  |                                                                                                                                                  
gAAGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAg                                                                                                                                                  >  8:67684/1‑149 (MQ=255)
 aaGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAgg                                                                                                                                                 >  6:169113/1‑149 (MQ=255)
     tAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAggcggc                                                                                                                                             <  5:171722/149‑1 (MQ=255)
       aGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGccc                                                                                                                                           <  5:317217/149‑1 (MQ=255)
       aGGGGCAGGCAGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGTTGGGGCTCTTCATGACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGGAAggcgggc                                                                                                                                            <  2:106129/148‑3 (MQ=255)
         gggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCCTTGCGGccccg                                                                                                                                         <  1:135771/149‑2 (MQ=255)
         gggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCa                                                                                                                                         <  7:67684/149‑1 (MQ=255)
               gCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGctcac                                                                                                                                   >  5:70125/1‑149 (MQ=255)
               gCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGctcac                                                                                                                                   >  6:317187/1‑149 (MQ=255)
                      aCTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAAc                                                                                                                            >  3:60420/1‑149 (MQ=255)
                       gtGCCAGCAACAGCAGGAACAGGGTACAGACGAGGCTGCTACAGGTGTGGATCAAGGTCCAGCGGCGGATGGTGGGGCTCTTCATCACTGTACCGATTTTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACt                                                                                                                           <  6:325394/148‑1 (MQ=255)
                       cTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACt                                                                                                                           >  5:317757/1‑149 (MQ=255)
                         gCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGAGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACtt                                                                                                                          <  4:157074/148‑1 (MQ=255)
                          ccAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGt                                                                                                                        <  6:30544/149‑1 (MQ=255)
                            aGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGa                                                                                                                              >  3:261319/1‑141 (MQ=255)
                                    cagGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCt                                                                                                              >  3:133225/1‑149 (MQ=255)
                                         aCAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTAGGTTGGAACCCGGGGCCGCAAGGCGCCCCCCTCTCACGCAGAACATGTAAGTTGCGCTGGCGa                                                                                                         >  8:152967/1‑149 (MQ=255)
                                          cAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCCAGGCGGCCCCCGCTAACACAGAACATGTAACTGGCGCAGGCgg                                                                                                         >  8:114053/1‑147 (MQ=255)
                                            gggTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACaa                                                                                                      <  4:133225/149‑1 (MQ=255)
                                             ggTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAAc                                                                                                     >  8:59934/1‑149 (MQ=255)
                                               tACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGc                                                                                                   >  3:326995/1‑149 (MQ=255)
                                                      cAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCGGCTCACTGAGAAGGTGTAAGTGGCGCTGGCGACAGCGCTgcgcgc                                                                                            >  8:219153/1‑147 (MQ=255)
                                                             ctCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGcc                                                                                      <  6:317757/148‑1 (MQ=255)
                                                              tCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCg                                                                                     >  1:257462/1‑148 (MQ=255)
                                                              tCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCACCACACTCAGAACATGTAAGTGGCGCTGGCGACAACGCTGCGCTCAACGCCg                                                                                     >  6:7414/1‑148 (MQ=255)
                                                                        gATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCAGCAAGGCGGCCCCAGCTCACTCAGAACTTTTAAGTGGCGCTTGCGCCAACTCTGCGCTCATCGCCGTAGTAGCac                                                                            >  5:163777/1‑146 (MQ=255)
                                                                           cAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGcc                                                                                      <  8:136940/134‑1 (MQ=255)
                                                                           cAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGcc                                                                                      >  7:136940/1‑134 (MQ=255)
                                                                              ggACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGa                                                                    <  4:60420/149‑1 (MQ=255)
                                                                              ggACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGa                                                                    <  5:169113/149‑1 (MQ=255)
                                                                                     cggcggATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACaa                                                             <  6:196277/149‑1 (MQ=255)
                                                                                      ggcggATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAg                                                            <  7:219153/149‑1 (MQ=255)
                                                                                        cggATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTa                                                          <  6:70125/149‑1 (MQ=255)
                                                                                        cggATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTa                                                          <  7:152967/149‑1 (MQ=255)
                                                                                           aTGGTGGGACTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTa                                                          >  3:298675/1‑146 (MQ=255)
                                                                                            tggtggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagaga                                                      <  4:272410/149‑1 (MQ=255)
                                                                                             ggtggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagaga                                                      <  7:114053/148‑1 (MQ=255)
                                                                                                 gggCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTaa                                                 <  4:326995/149‑1 (MQ=255)
                                                                                                     tctTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCt                                             >  5:149771/1‑149 (MQ=255)
                                                                                                      ctTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCtt                                            >  5:161895/1‑149 (MQ=255)
                                                                                                       ttCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCCGCTCACTCAGAACGTGTTAGTGGCGCTGGCGACAACGCGGCGCTCATCGCCGGAGTAGCAGTATTAGGAg                                                                   >  4:225822/1‑125 (MQ=255)
                                                                                                       ttCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAg                                                                   <  3:225822/125‑1 (MQ=255)
                                                                                                          atcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTc                                        <  7:59934/149‑1 (MQ=255)
                                                                                                             aCGGTACGTATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCTGTAGTAGGAGTCACAAGAAGAGAGGTAATCCTTGTTcagc                                     <  6:163777/149‑1 (MQ=255)
                                                                                                             aCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTcag                                      <  5:317187/148‑1 (MQ=255)
                                                                                                              cGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTcagca                                    <  4:298675/149‑1 (MQ=255)
                                                                                                               ggTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTcagcag                                   >  4:63021/1‑149 (MQ=255)
                                                                                                               ggTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTcagcag                                   <  6:149771/149‑1 (MQ=255)
                                                                                                                   ccGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtg                               <  5:7414/149‑1 (MQ=255)
                                                                                                                     gATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATc                                               <  8:236921/131‑1 (MQ=255)
                                                                                                                     gATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATc                                               >  7:236921/1‑131 (MQ=255)
                                                                                                                      aTTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAg                                                            >  8:260828/1‑117 (MQ=255)
                                                                                                                      aTTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAg                                                            <  7:260828/117‑1 (MQ=255)
                                                                                                                            gtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTg                      >  5:238839/1‑149 (MQ=255)
                                                                                                                             tCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGa                     <  3:94285/149‑1 (MQ=255)
                                                                                                                             tCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGa                     <  2:257462/149‑1 (MQ=255)
                                                                                                                             tCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGAAGTAGCAGTAGTAGGAGTTACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGa                     >  2:52889/1‑149 (MQ=255)
                                                                                                                                    tgagctggggCCGCCTTGCGGCCCCGGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgac              >  7:212059/7‑149 (MQ=255)
                                                                                                                                      agctggggCCGCCTTGCGGCCCCGGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATc                                               <  8:172618/110‑1 (MQ=255)
                                                                                                                                      agctggggCCGCCTTGCGGCCCCGGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATc                                               >  7:172618/5‑114 (MQ=255)
                                                                                                                                       cccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccga           >  2:226548/1‑149 (MQ=255)
                                                                                                                                        ccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccgac          <  3:63021/149‑1 (MQ=255)
                                                                                                                                          ggggCCGCCTTGCGGCCCCGGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGAcgc        <  8:212059/149‑1 (MQ=255)
                                                                                                                                                gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTc  >  3:239452/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAAGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTC  >  NC_002947/425202‑425494

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: ‑