Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,550,663 A→G 8.5% D302G (GAC→GGC)  PP_4881 → iron ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,550,6630AG 8.5% 156.4 / 4.6 60D302G (GAC→GGC) PP_4881iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (31/23);  new base G (3/2);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

CACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTT  >  NC_002947/5550517‑5550809
                                                                                                                                                   |                                                                                                                                                      
cACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGAcc                                                                                                                                                      >  4:353362/1‑149 (MQ=255)
 aCTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACcg                                                                                                                                                     >  2:277547/1‑149 (MQ=255)
   tATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACcgcg                                                                                                                                                   >  6:307929/1‑149 (MQ=255)
         cGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑g                                                                                                                                             >  2:15665/1‑148 (MQ=255)
         cGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑Ga                                                                                                                                            <  6:129739/149‑1 (MQ=255)
               aGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGt                                                                                                                                      <  1:277547/149‑1 (MQ=255)
                ggCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGACAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGtt                                                                                                                                     <  1:15665/149‑1 (MQ=255)
                   gTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTTGATGGCCAGCTAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAATCCGTTcg‑‑‑a                                                                                                                               <  6:82025/149‑1 (MQ=255)
                   gTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTcg‑‑‑a                                                                                                                               <  3:128615/149‑1 (MQ=255)
                   gTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATAACCCAGACAACCGCCGAGTACACACTGCACAAGGCCATGGTTTCAGGCCGCGGACA‑GAAAccatgca                                                                                                                                   >  7:69634/1‑143 (MQ=255)
                      gTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCA‑GAAAccatgc                                                                                                                                    >  1:277085/1‑140 (MQ=255)
                      gTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCA‑GAAAccatgc                                                                                                                                    <  2:277085/144‑5 (MQ=255)
                             ccAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCACACCACCGCCGACTACCCAATGCACAACGGCATGGTTTCTGACCGCCGCCT‑GAAACCGTTCG‑‑‑ACGAGCTgcgc                                                                                                                     >  3:47495/1‑149 (MQ=255)
                                         cGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGGAAACCGTTCG‑‑‑ACGAGCTGCGCccgccg                                                                                                               >  5:237600/1‑144 (MQ=255)
                                         cGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGAtc                                                                                                         >  1:83148/1‑149 (MQ=255)
                                          gAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGCCCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGAtct                                                                                                        >  7:139012/1‑149 (MQ=255)
                                           aaaGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCAGCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGAtctc                                                                                                       >  3:276028/1‑149 (MQ=255)
                                            aagaagCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCg                                                                                                      <  2:338766/149‑1 (MQ=255)
                                               aagCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGTCCGCGGCCT‑GAAACCGTTcg‑‑‑ac                                                                                                                              <  2:320942/122‑1 (MQ=255)
                                               aagCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTcg‑‑‑ac                                                                                                                              >  1:320942/1‑122 (MQ=255)
                                                  cccAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCt                                                                                                <  8:139012/149‑1 (MQ=255)
                                                    cAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGa                                                                                              >  3:104649/1‑149 (MQ=255)
                                                     aGGTCCTGCTCATCTGGATGGCCTGCGATGAAGGC‑CAGCGCGTGATGACCGAGACCACCGCCGAGTTCCCACTGCACAAGGGCATGGTTTCTGAGCGCTGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGGCTGCTGAc                                                                                             <  4:47495/149‑1 (MQ=255)
                                                     aGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCACAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTcg‑‑‑a                                                                                                                               >  1:53830/1‑116 (MQ=255)
                                                     aGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTcg‑‑‑a                                                                                                                               <  2:53830/115‑1 (MQ=255)
                                                       gCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTCATCACCCAGACCACCGCCGAGTACCCACTGCACAACGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCt                                                                                           >  8:262787/1‑149 (MQ=255)
                                                          cTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTggg                                                                                        <  2:83148/149‑1 (MQ=255)
                                                             cTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTggg                                                                                        <  6:237600/147‑1 (MQ=255)
                                                             cTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCaa                                                                                     >  2:111939/1‑149 (MQ=255)
                                                                aaCTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCA‑GAAAccatgc                                                                                                                                    >  1:210644/1‑98 (MQ=38)
                                                                aaCTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCA‑GAAAccatgc                                                                                                                                    <  2:210644/102‑5 (MQ=38)
                                                                     gATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                              >  3:309601/1‑148 (MQ=255)
                                                                           cAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCa                                                                       <  3:353362/149‑1 (MQ=255)
                                                                                   aagGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTCGGCATTGCCCAggt                                                                            >  6:39659/1‑135 (MQ=255)
                                                                                   aagGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTCGGCATTGCCCAggt                                                                            <  5:39659/136‑2 (MQ=255)
                                                                                     ggC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGa                                                             <  5:307929/149‑1 (MQ=255)
                                                                                         cAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCGCCCTCGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGa                                                             <  1:111939/149‑1 (MQ=255)
                                                                                                 gATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCg                                                  >  7:112057/1‑149 (MQ=255)
                                                                                                  aTCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCg                                                  <  4:309601/148‑1 (MQ=255)
                                                                                                       ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGc                                            >  5:363218/1‑149 (MQ=255)
                                                                                                             caccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtgat                                      >  7:85169/1‑149 (MQ=255)
                                                                                                             caccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtgat                                      <  7:262787/149‑1 (MQ=255)
                                                                                                                          ccACTGCACAAGGGCATGGTTTCAGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGCCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTCCTCTGATGACTCCCGCCCTg                         >  6:28781/1‑149 (MQ=255)
                                                                                                                           cACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGt                        <  8:112057/149‑1 (MQ=255)
                                                                                                                            aCTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTc                       <  4:276028/149‑1 (MQ=255)
                                                                                                                            aCTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTc                       >  1:210523/1‑149 (MQ=255)
                                                                                                                                 acaAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAAc                                                                 <  8:105675/102‑1 (MQ=255)
                                                                                                                                 acaAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAAc                                                                 >  7:105675/1‑102 (MQ=255)
                                                                                                                                 acaAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAAc                  >  1:232070/1‑149 (MQ=255)
                                                                                                                                 acaAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAAc                  >  7:121457/1‑149 (MQ=255)
                                                                                                                                    aGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCgg               >  2:79268/1‑149 (MQ=255)
                                                                                                                                       gCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGc            <  6:363218/149‑1 (MQ=255)
                                                                                                                                       gCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGc            >  6:371555/1‑149 (MQ=255)
                                                                                                                                         aTGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCg          >  6:204196/1‑149 (MQ=255)
                                                                                                                                          tGGTTTGTGACCGCGGCGT‑GATACCGTTCG‑‑‑ACGAGCTGCGCCGGCCGTAGATCTGGCCAGCTGACCTGTGCAATGCCGAGGAAGCCATGGAACTTGTGGGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGTAACGGTGcag          <  4:104649/148‑3 (MQ=255)
                                                                                                                                          tGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGGCGAAGATCTCGCCAGCTGACCTGGGCAATGGCGAGGTAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtgat                                      <  1:142433/120‑1 (MQ=255)
                                                                                                                                           ggTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCgg         >  1:213035/1‑148 (MQ=255)
                                                                                                                                                 tGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCtt  >  7:173905/1‑149 (MQ=255)
                                                                                                                                                 tGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCtt  >  7:86610/1‑149 (MQ=255)
                                                                                                                                                 tGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCtt  >  7:91560/1‑149 (MQ=255)
                                                                                                                                                   aCCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGc                                                           >  6:182357/1‑90 (MQ=255)
                                                                                                                                                   aCCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGc                                                           <  5:182357/90‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGC‑CAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCT‑GAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTT  >  NC_002947/5550517‑5550809

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