Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,550,672 T→A 7.8% L305Q (CTG→CAG)  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,6720TA 7.8% 172.8 / 4.7 64L305Q (CTG→CAG) PP_4881iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base A (3/2);  total (33/31)
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.26e-01

CGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAGCC  >  NC_002947/5550526‑5550820
                                                                                                                                                   |                                                                                                                                                        
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                                                                                                                                               <  2:277085/144‑5 (MQ=255)
             gTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatgc                                                                                                                                               >  1:277085/1‑140 (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                                                                                                                                               <  2:210644/102‑5 (MQ=38)
                                                       aaCTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatgc                                                                                                                                               >  1:210644/1‑98 (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:262787/149‑1 (MQ=255)
                                                                                                    caccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtgat                                                 >  7:85169/1‑149 (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                             >  1:232070/1‑149 (MQ=255)
                                                                                                                        acaAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAAc                             >  7:121457/1‑149 (MQ=255)
                                                                                                                           aGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCgg                          >  2:79268/1‑149 (MQ=255)
                                                                                                                              gCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGc                       <  6:363218/149‑1 (MQ=255)
                                                                                                                              gCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGc                       >  6:371555/1‑149 (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)
                                                                                                                                          aCCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGc                                                                      <  5:182357/90‑1 (MQ=255)
                                                                                                                                          aCCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGc                                                                      >  6:182357/1‑90 (MQ=255)
                                                                                                                                           ccGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGt                   >  7:129826/1‑140 (MQ=255)
                                                                                                                                           ccGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGt                   <  8:129826/140‑1 (MQ=255)
                                                                                                                                             gcgGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTAc        <  2:232070/149‑1 (MQ=255)
                                                                                                                                                 ccTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAg    <  8:85169/149‑1 (MQ=255)
                                                                                                                                                  cTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAGc   <  8:173905/149‑1 (MQ=255)
                                                                                                                                                   tGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGAAGCGTACCggcac  >  5:336995/1‑145 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
CGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCG‑‑‑ACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAGCC  >  NC_002947/5550526‑5550820

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: