Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,550,661 T→A 10.1% S301S (TCT→TCA 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,6610TA10.1% 157.8 / 5.1 57S301S (TCT→TCAPP_4881iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (27/24);  new base A (4/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTT  >  NC_002947/5550515‑5550809
                                                                                                                                                   |                                                                                                                                                        
gTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCt                                                                                                                                                          <  3:65297/147‑1 (MQ=255)
  cACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGAcc                                                                                                                                                      >  4:353362/1‑149 (MQ=255)
   aCTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACcg                                                                                                                                                     >  2:277547/1‑149 (MQ=255)
     tATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACcgcg                                                                                                                                                   >  6:307929/1‑149 (MQ=255)
           cGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTg                                                                                                                                             >  2:15665/1‑148 (MQ=255)
           cGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGa                                                                                                                                            <  6:129739/149‑1 (MQ=255)
                 aGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGt                                                                                                                                      <  1:277547/149‑1 (MQ=255)
                  ggCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGACAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGtt                                                                                                                                     <  1:15665/149‑1 (MQ=255)
                     gTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTTGATGGCCAGCTAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAATCCGTTcg‑‑‑a                                                                                                                               <  6:82025/149‑1 (MQ=255)
                     gTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTcg‑‑‑a                                                                                                                               <  3:128615/149‑1 (MQ=255)
                     gTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATAACCCAGACAACCGCCGAGTACACACTGCACAAGGCCATGGTTTCAGGCCGCGGACAGAAAccatgca                                                                                                                                   >  7:69634/1‑143 (MQ=255)
                        gTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatgc                                                                                                                                    >  1:277085/1‑140 (MQ=255)
                        gTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatgc                                                                                                                                    <  2:277085/144‑5 (MQ=255)
                               ccAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCACACCACCGCCGACTACCCAATGCACAACGGCATGGTTTCTGACCGCCGCCTGAAACCGTTCG‑‑‑ACGAGCTgcgc                                                                                                                     >  3:47495/1‑149 (MQ=255)
                                           cGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGGAAACCGTTCG‑‑‑ACGAGCTGCGCccgccg                                                                                                               >  5:237600/1‑144 (MQ=255)
                                           cGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGAtc                                                                                                         >  1:83148/1‑149 (MQ=255)
                                            gAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGCCCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGAtct                                                                                                        >  7:139012/1‑149 (MQ=255)
                                             aaaGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCAGCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGAtctc                                                                                                       >  3:276028/1‑149 (MQ=255)
                                              aagaagCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCg                                                                                                      <  2:338766/149‑1 (MQ=255)
                                                 aagCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGTCCGCGGCCTGAAACCGTTcg‑‑‑ac                                                                                                                              <  2:320942/122‑1 (MQ=255)
                                                 aagCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTcg‑‑‑ac                                                                                                                              >  1:320942/1‑122 (MQ=255)
                                                    cccAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCt                                                                                                <  8:139012/149‑1 (MQ=255)
                                                      cAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGa                                                                                              >  3:104649/1‑149 (MQ=255)
                                                       aGGTCCTGCTCATCTGGATGGCCTGCGATGAAGGCCAGCGCGTGATGACCGAGACCACCGCCGAGTTCCCACTGCACAAGGGCATGGTTTCTGAGCGCTGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGGCTGCTGAc                                                                                             <  4:47495/149‑1 (MQ=255)
                                                       aGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCACAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTcg‑‑‑a                                                                                                                               >  1:53830/1‑116 (MQ=255)
                                                       aGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTcg‑‑‑a                                                                                                                               <  2:53830/115‑1 (MQ=255)
                                                         gCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTCATCACCCAGACCACCGCCGAGTACCCACTGCACAACGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCt                                                                                           >  8:262787/1‑149 (MQ=255)
                                                            cTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTggg                                                                                        <  2:83148/149‑1 (MQ=255)
                                                               cTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTggg                                                                                        <  6:237600/147‑1 (MQ=255)
                                                               cTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCaa                                                                                     >  2:111939/1‑149 (MQ=255)
                                                                  aaCTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatgc                                                                                                                                    >  1:210644/1‑98 (MQ=38)
                                                                  aaCTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatgc                                                                                                                                    <  2:210644/102‑5 (MQ=38)
                                                                       gATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                              >  3:309601/1‑148 (MQ=255)
                                                                             cAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCa                                                                       <  3:353362/149‑1 (MQ=255)
                                                                                     aagGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTCGGCATTGCCCAggt                                                                            >  6:39659/1‑135 (MQ=255)
                                                                                     aagGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTCGGCATTGCCCAggt                                                                            <  5:39659/136‑2 (MQ=255)
                                                                                       ggCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGa                                                             <  5:307929/149‑1 (MQ=255)
                                                                                           cAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGCCCTCGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGa                                                             <  1:111939/149‑1 (MQ=255)
                                                                                                   gATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCg                                                  >  7:112057/1‑149 (MQ=255)
                                                                                                    aTCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCg                                                  <  4:309601/148‑1 (MQ=255)
                                                                                                         ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGc                                            >  5:363218/1‑149 (MQ=255)
                                                                                                               caccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtgat                                      >  7:85169/1‑149 (MQ=255)
                                                                                                               caccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtgat                                      <  7:262787/149‑1 (MQ=255)
                                                                                                                            ccACTGCACAAGGGCATGGTTTCAGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGCCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTCCTCTGATGACTCCCGCCCTg                         >  6:28781/1‑149 (MQ=255)
                                                                                                                             cACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGt                        <  8:112057/149‑1 (MQ=255)
                                                                                                                              aCTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTc                       <  4:276028/149‑1 (MQ=255)
                                                                                                                              aCTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTc                       >  1:210523/1‑149 (MQ=255)
                                                                                                                                   acaAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAAc                                                                 >  7:105675/1‑102 (MQ=255)
                                                                                                                                   acaAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAAc                                                                 <  8:105675/102‑1 (MQ=255)
                                                                                                                                   acaAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAAc                  >  7:121457/1‑149 (MQ=255)
                                                                                                                                   acaAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAAc                  >  1:232070/1‑149 (MQ=255)
                                                                                                                                      aGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCgg               >  2:79268/1‑149 (MQ=255)
                                                                                                                                         gCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGc            >  6:371555/1‑149 (MQ=255)
                                                                                                                                         gCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGc            <  6:363218/149‑1 (MQ=255)
                                                                                                                                           aTGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCg          >  6:204196/1‑149 (MQ=255)
                                                                                                                                            tGGTTTGTGACCGCGGCGTGATACCGTTCG‑‑‑ACGAGCTGCGCCGGCCGTAGATCTGGCCAGCTGACCTGTGCAATGCCGAGGAAGCCATGGAACTTGTGGGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGTAACGGTGcag          <  4:104649/148‑3 (MQ=255)
                                                                                                                                            tGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGGCGAAGATCTCGCCAGCTGACCTGGGCAATGGCGAGGTAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtgat                                      <  1:142433/120‑1 (MQ=255)
                                                                                                                                             ggTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCgg         >  1:213035/1‑148 (MQ=255)
                                                                                                                                                   tGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCtt  >  7:173905/1‑149 (MQ=255)
                                                                                                                                                   tGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCtt  >  7:86610/1‑149 (MQ=255)
                                                                                                                                                   tGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCtt  >  7:91560/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
GTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTT  >  NC_002947/5550515‑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: