Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,412,330 G→A 6.0% P614S (CCC→TCC)  PP_4755 ← ferrichrome‑iron receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,412,3300GA 6.0% 194.9 / 2.6 66P614S (CCC→TCC) PP_4755ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base A (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

CCTTGGCTTCCAGGTCAAGGCCACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGATATAGGTCAGTGCCGCGCGGTAGGAGGTTTCGG  >  NC_002947/5412183‑5412476
                                                                                                                                                   |                                                                                                                                                    
ccTTGGCTTCCAGGTCAAGGCCACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATggg                                                                                                                                                      <  6:210878/148‑1 (MQ=255)
ccTTGGCTTCCAGGTCAAGGCCACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGc                                                                                                                                                     >  7:14281/1‑149 (MQ=255)
ccTTGGCTTCCAGGGCAAGGCCACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGc                                                                                                                                                     <  6:189410/149‑1 (MQ=255)
  ttGGCTTCCAGGTCAAGGCCACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCg                                                                                                                                                    >  1:165277/1‑148 (MQ=255)
    ggCTTCCAGGTCAAGGCCACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGaa                                                                                                                                                  <  8:76719/148‑1 (MQ=255)
    ggCTTCCAGGTCAAGGCCACGCACGCGTGACTCGCCCACCGGCTGATGCTCGATGATGCCGCTTGGCAACACCACGGCAATTGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGaa                                                                                                                                                  <  4:61512/148‑1 (MQ=255)
         ccAGGTCTAGGCGACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTGGGCAACACCACGGCAATGGTGACGTTCTCTTGGGGCAGGTCGTAAACAGCGGCGGAGAAGAGGGCGTCCATGCCCATGGGCGAATACt                                                                                                                                              <  3:168378/147‑1 (MQ=255)
                aaGGCCACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCa                                                                                                                                     <  1:98203/149‑1 (MQ=255)
                   gCCACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACtt                                                                                                                                  >  8:21576/1‑149 (MQ=255)
                   gCCACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACt                                                                                                                                   >  1:258774/1‑148 (MQ=255)
                      acgcacgcGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTACATGCCCATGGGCGAATACTTGATGCCCACTTCg                                                                                                                                >  6:328023/1‑148 (MQ=255)
                       cgcacgcGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATTGGCGAATACTTGATGCCCACTTCGTa                                                                                                                              >  4:238114/1‑149 (MQ=255)
                        gcacgcGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATTGGCGAATACTTGATGCCCACTTCGTAt                                                                                                                             <  3:238114/149‑1 (MQ=255)
                                 aCTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGcc                                                                                                                     <  2:165277/148‑1 (MQ=255)
                                  cTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCgg                                                                                                                   <  2:258774/149‑1 (MQ=255)
                                     gCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCt                                                                                                                >  5:296395/1‑149 (MQ=255)
                                     gCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGCTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGt                                                                                                                  <  1:101236/147‑1 (MQ=255)
                                        cACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCg                                                                                                             >  5:271843/1‑149 (MQ=255)
                                                 ctgctCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCTAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAAc                                                                                                      <  1:232787/147‑1 (MQ=255)
                                                  tgctCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACt                                                                                                                                              >  6:32106/1‑106 (MQ=255)
                                                  tgctCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACt                                                                                                                                              <  5:32106/106‑1 (MQ=255)
                                                  tgctCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGcc                                                                                                   >  4:215415/1‑149 (MQ=255)
                                                       gatgatGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTa                                                                                                                              <  6:127547/117‑1 (MQ=255)
                                                       gatgatGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTa                                                                                                                              >  5:127547/1‑117 (MQ=255)
                                                       gatgatGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGc                                                                                                              <  5:36512/133‑1 (MQ=255)
                                                       gatgatGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGc                                                                                                              >  6:36512/1‑133 (MQ=255)
                                                               cGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTgcggcggc                                                                                      >  3:157614/1‑149 (MQ=255)
                                                                  ttGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTgcggcggc                                                                                      >  7:139675/1‑146 (MQ=255)
                                                                    ggCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGa                                                                                                                                           >  4:200645/1‑91 (MQ=255)
                                                                    ggCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGa                                                                                                                                           <  3:200645/91‑1 (MQ=255)
                                                                    ggCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGACGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGAc                                                                                 >  6:186468/1‑149 (MQ=255)
                                                                    ggCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCACATGGACGAATACTTGATGCCCACTTCGTATTGCTTGCAGGTCTGCGCGGAAACGCCCACTTCCGGCGTCGAGAc                                                                                 >  6:331546/1‑149 (MQ=255)
                                                                        acacCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTACATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGAc                                                                             <  5:328023/149‑1 (MQ=255)
                                                                          accacGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCTGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑c                                                                         <  8:14281/149‑1 (MQ=255)
                                                                            cacGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑cca                                                                       >  6:248920/1‑149 (MQ=255)
                                                                            cacGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑cca                                                                       <  6:271843/149‑1 (MQ=255)
                                                                               ggCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑ccacca                                                                    <  6:296395/149‑1 (MQ=255)
                                                                                gCAATGGTGACGTTCTGTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCgg                                                                                                            <  1:114094/110‑1 (MQ=255)
                                                                                gCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGAGAAAACTTGATGCCCACATCGTAATGATTGCCGGTCTGCgg                                                                                                            >  2:114094/1‑110 (MQ=255)
                                                                                 cAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGATTACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACAATg                                                                  >  8:296076/1‑149 (MQ=255)
                                                                                  aaTGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGc                                                                 >  1:16132/1‑149 (MQ=255)
                                                                                    tGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTc                                                               <  3:215415/149‑1 (MQ=255)
                                                                                          cGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAg                                                         >  4:219227/1‑149 (MQ=255)
                                                                                                 ttGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑ccacca                                                                    <  8:95714/131‑1 (MQ=255)
                                                                                                 ttGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑ccacca                                                                    >  7:95714/1‑131 (MQ=255)
                                                                                                   gggTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGACGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGaca                                                  <  5:186468/147‑1 (MQ=255)
                                                                                                   gggTCAGGTCGAAAACAGCGGCGGAGAAGAGTGCGTACATGCCCATGGACGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGTGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACTTAGGTCGaca                                                  <  5:331546/147‑1 (MQ=255)
                                                                                                        aGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATc                                           <  2:16132/149‑1 (MQ=255)
                                                                                                         ggTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCg                                          <  8:139675/149‑1 (MQ=255)
                                                                                                          gTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGt                                         >  4:214145/1‑149 (MQ=255)
                                                                                                           tCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGtt                                        <  3:219227/149‑1 (MQ=255)
                                                                                                             gTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGa                                      <  4:279205/149‑1 (MQ=255)
                                                                                                                 aCAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGa                                   >  4:109933/1‑148 (MQ=255)
                                                                                                                 aCAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGa                                   >  8:120189/1‑148 (MQ=255)
                                                                                                                    gcggcgGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCgaga                                                                                  >  5:200474/1‑100 (MQ=255)
                                                                                                                    gcggcgGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCgaga                                                                                  <  6:200474/100‑1 (MQ=255)
                                                                                                                       gcggAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑GCACCATGCTCACGTAGGTCGacac                                                 <  7:200587/128‑1 (MQ=255)
                                                                                                                       gcggAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGacac                                                 >  8:200587/1‑128 (MQ=255)
                                                                                                                               gagTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCTGTGGGGTGACGGTGTGACCACCATGCTCACGTAGGTCGAAACTTCATCGTTGACGATATAGGTCAGTGc                      >  1:265200/1‑149 (MQ=255)
                                                                                                                               gagTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCAAGTAGGTCGACACTTCATCGTTGACGATATAGGTCAGTGCcg                    >  4:217754/1‑149 (MQ=255)
                                                                                                                                 gTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGATATAGGTCAGTGCcgcg                  >  2:175845/1‑149 (MQ=255)
                                                                                                                                     gTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGATATAGGTCAGTGCCGCGCgg               >  7:118766/1‑148 (MQ=255)
                                                                                                                                         aTGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAgg                                                        <  8:231756/103‑1 (MQ=255)
                                                                                                                                         aTGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAgg                                                        >  7:231756/1‑103 (MQ=255)
                                                                                                                                          tGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGATATAGGTCAGTGCCGCGCGGTaggag         <  7:21576/149‑1 (MQ=255)
                                                                                                                                              cATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGATATAGGTCAGTGCCGCGCGGTAGGAGGttt     >  5:55230/1‑149 (MQ=255)
                                                                                                                                                tGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGATATAGGTCAGTGCCGCGCGGTAGGAGGTTTCg   <  8:118766/149‑1 (MQ=255)
                                                                                                                                                 gggCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGATATAGGTCAGTGCCGCGCGGTAGGAGGTTTCg   >  6:269305/1‑148 (MQ=255)
                                                                                                                                                  ggCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGATATAGGTCAGTGCCGCGCGGTAGGAGGTTTCgg  <  3:109933/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTTGGCTTCCAGGTCAAGGCCACGCACGCGTGACTCGCCCACCGTCTGCTGCTCGATGATGCCGCTTGGCAACACCACGGCAATGGTGACGTTCTCTTGGGTCAGGTCGTAAACAGCGGCGGAGAAGAGTGCGTCCATGCCCATGGGCGAATACTTGATGCCCACTTCGTATTGCTTGCCGGTCTGCGGGGTAACGCCCACTTGCGGCGGCGAGACGGACT‑‑CCACCATGCTCACGTAGGTCGACACTTCATCGTTGACGATATAGGTCAGTGCCGCGCGGTAGGAGGTTTCGG  >  NC_002947/5412183‑5412476

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