Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,333,094 T→A 6.4% intergenic (+34/+50) PP_3804 → / ← yaiW hypothetical protein/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,333,0940TA 6.4% 185.3 / 3.0 62intergenic (+34/+50)PP_3804/yaiWhypothetical protein/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base A (2/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

ATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCGGGTAATCTTCGGGCTCTTCAGC  >  NC_002947/4332950‑4333234
                                                                                                                                                |                                                                                                                                             
aTCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGcc                                                                                                                                           >  2:190318/1‑149 (MQ=255)
 tCAAGGCGGCCGGGGCCATGCTGGTGGTGGGTGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGccc                                                                                                                                          <  8:76120/149‑1 (MQ=255)
   aaGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGccc                                                                                                                                          >  7:76120/1‑147 (MQ=255)
     ggcggcCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCgtgt                                                                                                                                      <  3:73097/149‑1 (MQ=255)
     ggcggcCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCgtgt                                                                                                                                      >  7:35123/1‑149 (MQ=255)
           cGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                                <  1:190318/149‑1 (MQ=255)
                 cAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAAc                                                                                                                          >  8:64179/1‑149 (MQ=255)
                   aGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAAccc                                                                                                                        <  2:323734/149‑1 (MQ=255)
                         tggtggtggAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCaa                                                                                                                  <  3:356050/149‑1 (MQ=255)
                              gtggAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATc                                                                                                             <  6:106062/149‑1 (MQ=255)
                               tggAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATc                                                                                                             <  6:186312/148‑1 (MQ=255)
                               tggAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATc                                                                                                             >  5:186312/1‑148 (MQ=255)
                                ggAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATcgc                                                                                                           >  5:283948/1‑149 (MQ=255)
                                  agagTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCgg                                                                                                         <  7:123827/149‑1 (MQ=255)
                                    agTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGtt                                                                                                       >  6:97036/1‑149 (MQ=255)
                                     gTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTtg                                                                                                      <  7:214602/149‑1 (MQ=255)
                                     gTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTtg                                                                                                      <  8:300924/149‑1 (MQ=255)
                                        aaccTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGGGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTtgct                                                                                                    <  3:13148/146‑1 (MQ=255)
                                         cccTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCggg                                                                                                                               <  8:284354/120‑1 (MQ=255)
                                         cccTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCggg                                                                                                                               >  7:284354/1‑120 (MQ=255)
                                         cccTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTtgctgc                                                                                                  <  7:220666/149‑1 (MQ=255)
                                           cTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCa                                                                                                <  6:283948/149‑1 (MQ=255)
                                                gATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCGAAGCCCgtgc                                                                                                                                      >  3:198407/1‑105 (MQ=255)
                                                gATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCGAAGCCCgtgc                                                                                                                                      >  3:203018/1‑105 (MQ=255)
                                                gATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCGAAGCCCgtgc                                                                                                                                      <  4:198407/106‑2 (MQ=255)
                                                gATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCGAAGCCCgtgc                                                                                                                                      <  4:203018/106‑2 (MQ=255)
                                                 aTGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTc                                                                                                                                              >  6:44224/1‑97 (MQ=255)
                                                 aTGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTc                                                                                                                                              <  5:44224/97‑1 (MQ=255)
                                                   gCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACg                                                                                        >  1:65339/1‑149 (MQ=255)
                                                     gATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCt                                                                                      >  6:299080/1‑149 (MQ=255)
                                                          cGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATg                                                                                 >  6:317266/1‑149 (MQ=255)
                                                          cGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATg                                                                                 <  5:317266/149‑1 (MQ=255)
                                                           ggggCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGc                                                                                >  6:101858/1‑149 (MQ=255)
                                                             ggCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑c                                                                             <  7:64179/149‑1 (MQ=255)
                                                                     gCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCCCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATa                                                                     <  7:253999/149‑1 (MQ=255)
                                                                      ccAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCa                                                                                                <  3:170903/122‑1 (MQ=255)
                                                                      ccAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCa                                                                                                >  4:170903/1‑122 (MQ=255)
                                                                       cAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAgc                                                                   >  5:201155/1‑149 (MQ=255)
                                                                         gCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAgcgc                                                                 <  6:271795/149‑1 (MQ=255)
                                                                           gTTGACAAAGTGTTTTTGGAATTGTTCTGGCGCTGATATGACAGGCCCGCGCGGTGGATGGCACGAGATTAGCACGTGTTCGCGGGTAAACCCGATAAAGGATCGAGGTTGCTGCAATGCTAACGCTTCATGCA‑CCGCTGATAGctct                                                                >  4:133821/1‑148 (MQ=255)
                                                                                 cAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCAc                                                         >  6:227312/1‑149 (MQ=255)
                                                                                 cAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCAc                                                         <  4:139008/149‑1 (MQ=255)
                                                                                        ttGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCAcccgcccg                                                  >  1:334424/1‑149 (MQ=255)
                                                                                           gggCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGc                                                <  8:35123/148‑1 (MQ=255)
                                                                                              cATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAAc                                            <  5:101858/149‑1 (MQ=255)
                                                                                               aTTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTTGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAAc                                            >  7:297806/1‑148 (MQ=255)
                                                                                                       ggcgCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTg                                   <  5:97036/149‑1 (MQ=255)
                                                                                                           cTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCa                               <  8:194728/149‑1 (MQ=255)
                                                                                                              aTCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCt                            >  6:123553/1‑149 (MQ=255)
                                                                                                                   aCAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCACCCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTg                                   <  1:69395/138‑1 (MQ=255)
                                                                                                                   aCAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTg                                   >  2:69395/1‑137 (MQ=255)
                                                                                                                     aGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCggg                      >  4:89901/1‑148 (MQ=255)
                                                                                                                       gCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCGGGTaa                   >  1:352722/1‑149 (MQ=255)
                                                                                                                         ccGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCg                        >  2:237150/1‑142 (MQ=255)
                                                                                                                         ccGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCg                        <  1:237150/142‑1 (MQ=255)
                                                                                                                            cgcgGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCGGGTAATCTTc              >  2:250142/1‑149 (MQ=255)
                                                                                                                              cgGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCGGGTAATCTTCgg            <  4:91490/149‑1 (MQ=255)
                                                                                                                              cgGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCGGGTAATCTTCgg            <  5:227312/149‑1 (MQ=255)
                                                                                                                                     tGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCGGGTAATCTTCGGGCTCTTc     <  6:201155/149‑1 (MQ=255)
                                                                                                                                        cACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCGGGTAATCTTCGGGCTCTTCAGc  >  5:219604/1‑149 (MQ=255)
                                                                                                                                         aCGGGCTTCGCCCGTGTTCACGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCGGGTAATCTTCGGGCTCTTCa    <  5:111664/146‑1 (MQ=255)
                                                                                                                                         aCGGGCTTCGCCCGTGTTCACGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCGGGTAATCTTCGGGCTCTTCa    >  6:111664/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                             
ATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGCCAGGCTAACGCTTCATGCA‑CCGCTGATAGCGCTCATCCACCCGCCCGGCAAACCAGGCCGTGGTCAGCTTGCGGGTAATCTTCGGGCTCTTCAGC  >  NC_002947/4332950‑4333234

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