Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,326,469 T→C 8.2% intergenic (‑38/+65) PP_4686 ← / ← hmuV putative lipoprotein/hemin import ATP‑binding protein HmuV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,326,4690TC 8.2% 247.5 / 7.2 85intergenic (‑38/+65)PP_4686/hmuVputative lipoprotein/hemin import ATP‑binding protein HmuV
Reads supporting (aligned to +/- strand):  ref base T (41/37);  new base C (2/5);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGCGACCTTCACTGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCGATACCGTACACCGCCTTCAGGGC  >  NC_002947/5326321‑5326614
                                                                                                                                                    |                                                                                                                                                 
cgtgcgACCTTCACTGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCt                                                                                                                                                   <  3:44708/149‑1 (MQ=255)
 gtgcgACCTTCACTGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGGCGCATCGGCAGTTACTCCATGTCATCCACACGTTAAATGGCCTg                                                                                                                                                  >  6:366251/1‑149 (MQ=255)
       ccTTCACTGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCgg                                                                                                                                            <  6:135700/149‑1 (MQ=255)
        cTTCACTGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCggg                                                                                                                                           <  8:295826/149‑1 (MQ=255)
        cTTCACTGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGAAAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCggg                                                                                                                                           <  1:348363/149‑1 (MQ=255)
            aCTGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGgtgcg                                                                                                                                       <  8:213106/149‑1 (MQ=255)
            aCTGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGgggcg                                                                                                                                       <  1:15059/149‑1 (MQ=255)
            aCTGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGgggcg                                                                                                                                       >  3:249867/1‑149 (MQ=255)
              tGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGgggcggg                                                                                                                                     >  4:109180/1‑149 (MQ=255)
                ctctGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGcc                                                                                                                                   >  7:117661/1‑149 (MQ=255)
                 tctGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCATATTGCCAGTACCAGGCACAGCGAGAAACCGGACATCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCg                                                                                                                                  >  5:269884/1‑149 (MQ=255)
                     ccAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCgg                                                                                                                                 <  6:168008/146‑1 (MQ=255)
                     ccAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCCGCCCCCGGTGCAGGccatt                                                                                                                                <  6:295698/147‑4 (MQ=255)
                      cAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCgggg                                                                                                                                          >  6:338909/1‑136 (MQ=255)
                      cAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCgggg                                                                                                                                          <  5:338909/136‑1 (MQ=255)
                       aGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCCGCCCCCGGTGCAGGccatttaaa                                                                                                                            <  7:330324/149‑8 (MQ=255)
                       aGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAACGGCCCGCCCCCGGTGCAGGccatttaaa                                                                                                                            <  1:145926/149‑8 (MQ=255)
                             gCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCCGCCCCCGGTGCAGGccat                                                                                                                                 >  1:80689/1‑136 (MQ=255)
                             gCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCCGCCCCCGGTGCAGGccat                                                                                                                                 <  2:80689/138‑3 (MQ=255)
                              cAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGAc                                                                                                                        >  2:207699/1‑146 (MQ=255)
                              cAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGAc                                                                                                                        <  1:207699/146‑1 (MQ=255)
                              cAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAg                                                                                                                     <  4:255226/149‑1 (MQ=255)
                               aaCGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAgg                                                                                                                    >  7:68609/1‑149 (MQ=255)
                                 cggcggCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTc                                                                                                                  <  2:126483/149‑1 (MQ=255)
                                          gCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCcaca                                                                                                         <  2:242510/149‑1 (MQ=255)
                                          gCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCcaca                                                                                                         <  8:38674/149‑1 (MQ=255)
                                              gCCTGACAGCCGCCTATTGCCAGTACCAGGCTCAGCGAGAAACCGGACAGCAGCGGGTGTCTCATCTGCAGTTACTCCATGTCATCCACTCTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAg                                                                                                      <  6:269884/148‑1 (MQ=255)
                                              gCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGgggcgggc                                                                                                                                    >  7:305662/1‑118 (MQ=255)
                                              gCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGgggcgggc                                                                                                                                    <  8:305662/118‑1 (MQ=255)
                                              gCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAg                                                                                                      <  8:68609/148‑1 (MQ=255)
                                              gCCTGACAGCCGCCTATTGCCAGGACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACATTAAATGGCCTGCACCGGGGGCGGCCCGGTCCCAAGACAAGGTCGCTACCACAAAgc                                                                                                     >  7:171930/1‑148 (MQ=255)
                                                cTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGcc                                                                                                   >  4:143651/1‑149 (MQ=255)
                                                   aCAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTg                                                                                                <  1:210873/149‑1 (MQ=255)
                                                   aCAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTg                                                                                                >  5:227502/1‑149 (MQ=255)
                                                    cAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGc                                                                                               <  6:227502/149‑1 (MQ=255)
                                                       ccgccTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTCCCACAAAGGCCGTGCACg                                                                                            <  8:212022/149‑1 (MQ=255)
                                                          ccTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGtt                                                                                         <  3:143651/149‑1 (MQ=255)
                                                          ccTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGtt                                                                                         >  4:361247/1‑149 (MQ=255)
                                                             aTTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACg                                                                                            >  4:107295/1‑143 (MQ=255)
                                                             aTTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACg                                                                                            <  3:107295/143‑1 (MQ=255)
                                                                 ccAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAAt                                                                                  >  1:357254/1‑149 (MQ=255)
                                                                 ccAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAAt                                                                                  >  1:181095/1‑149 (MQ=255)
                                                                   aGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCCGCCCCCGGTGCAGGccatt                                                                                                                                <  6:119693/101‑4 (MQ=255)
                                                                   aGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCCGCCCCCGGTGCAGGccatt                                                                                                                                >  5:119693/1‑98 (MQ=255)
                                                                    gTACCAGGCACAGCGAGAAACCGGACAGTAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCAACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCaa                                                                               >  3:216269/1‑149 (MQ=255)
                                                                    gTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCaa                                                                               >  6:63651/1‑149 (MQ=255)
                                                                     tACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAAc                                                                              >  3:227029/1‑149 (MQ=255)
                                                                        cAGGCACAGCGAGAAACCGGACACCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCTAGACAAGGTCGCATCCACAAAGGCCGTCCACGGTTCGAGAATCAACgct                                                                           >  4:64043/1‑147 (MQ=255)
                                                                          ggCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGt                                                                          >  4:234194/1‑148 (MQ=255)
                                                                           gCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTg                                                                                                <  7:48659/125‑1 (MQ=255)
                                                                           gCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTg                                                                                                >  8:48659/1‑125 (MQ=255)
                                                                             acaGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGtgat                                                                       >  3:225253/1‑148 (MQ=255)
                                                                                 cGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTgatgat                                                                    >  1:86565/1‑147 (MQ=255)
                                                                                 cGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTgatgat                                                                    >  2:106362/1‑147 (MQ=255)
                                                                                      aaCCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCgg                                                              <  4:249867/148‑1 (MQ=255)
                                                                                       aCCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGAt                                                            <  2:125290/149‑1 (MQ=255)
                                                                                              agcagcGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGCGGGCGGGCCGGTCGAAAGACAAGGACGCTTCCCCAAACGCCGTGCCCGGTTCGCGAATCAACGGGTGATGATCAGCGGag                                                            >  4:206440/1‑141 (MQ=255)
                                                                                                cagcGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGt                                                                          >  1:127615/1‑126 (MQ=255)
                                                                                                cagcGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGt                                                                          <  2:127615/126‑1 (MQ=255)
                                                                                                cagcGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCt                                                   <  2:357254/149‑1 (MQ=255)
                                                                                                    gggTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCggg                                               <  3:109180/149‑1 (MQ=255)
                                                                                                    gggTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGGTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCggg                                               <  2:181095/149‑1 (MQ=255)
                                                                                                           gCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGcc                                        >  4:344813/1‑149 (MQ=255)
                                                                                                            cATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCt                                       <  7:220412/149‑1 (MQ=255)
                                                                                                              tCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGc                                     <  8:117661/149‑1 (MQ=255)
                                                                                                                 gCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCAc                                   >  7:93111/1‑148 (MQ=255)
                                                                                                                       aCTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTAc                             <  8:93111/148‑1 (MQ=255)
                                                                                                                        cTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGa                                                             >  4:277227/1‑115 (MQ=255)
                                                                                                                        cTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGt                           >  1:287013/1‑149 (MQ=255)
                                                                                                                        cTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGt                           >  6:9368/1‑149 (MQ=255)
                                                                                                                         tCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGTCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTc                          >  2:213594/1‑149 (MQ=255)
                                                                                                                           cgtGTCATGCAGACTTTTAAAGGACTGGAGCGGGGGCGGGCCGGTCGCATGTCAATGTTGCTTCCTCAAAGGCGGTTCAGGGTTCTAGAATCATCGTGTGATGATCAGCGGATGCCCGCGCTCTGGGTGCGCCTGGACGAGTACGTc                          <  4:227029/145‑1 (MQ=255)
                                                                                                                           cATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCg                         <  5:63651/148‑1 (MQ=255)
                                                                                                                            aTGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCGa                        >  2:187042/1‑148 (MQ=255)
                                                                                                                                   ccACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCTGGGTGCGCCTGCACCAGTACGTCGCTACCGTa                 >  6:279153/1‑148 (MQ=255)
                                                                                                                                   ccACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCGATACCGTa                 >  7:43415/1‑148 (MQ=255)
                                                                                                                                     acacTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCGATACCGTacac              >  8:279366/1‑149 (MQ=255)
                                                                                                                                        cTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCGATACCGTACACCGc           <  3:125450/149‑1 (MQ=255)
                                                                                                                                         tttAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGt                           <  4:220566/132‑1 (MQ=255)
                                                                                                                                         tttAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGt                           >  3:220566/1‑132 (MQ=255)
                                                                                                                                         tttAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCGATACCGTACACCGcc          >  5:192340/1‑149 (MQ=255)
                                                                                                                                             aaTGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCGATACCGTACACCGCCTTCa      <  5:9368/149‑1 (MQ=255)
                                                                                                                                               tGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCGATACCGTACACCGCCTTCAgg    <  8:43415/149‑1 (MQ=255)
                                                                                                                                                 gCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCGATACCGTACACCGCCTTCAGGGc  >  7:202544/1‑149 (MQ=255)
                                                                                                                                                   cTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCAc                                   <  2:260044/114‑1 (MQ=255)
                                                                                                                                                   cTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCAc                                   >  1:260044/1‑114 (MQ=255)
                                                                                                                                                    tGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGtga                                                                        >  5:180487/1‑76 (MQ=255)
                                                                                                                                                    tGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGtga                                                                        <  6:180487/76‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGTGCGACCTTCACTGCTCTGCCAGGCGGGCAACGGCGGCAGGCTGGCCTGACAGCCGCCTATTGCCAGTACCAGGCACAGCGAGAAACCGGACAGCAGCGGGTGTCGCATCGGCAGTTACTCCATGTCATCCACACTTTAAATGGCCTGCACCGGGGGCGGGCCGGTCGCAAGACAAGGTCGCTTCCACAAAGGCCGTGCACGGTTCGAGAATCAACGGGTGATGATCAGCGGATGCCCGCGCTCGGGGTGCGCCTGCACCAGTACGTCGATACCGTACACCGCCTTCAGGGC  >  NC_002947/5326321‑5326614

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