Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,713,742 T→G 6.1% intergenic (+11/‑91) tsaB → / → PP_1508 endopeptidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,713,7420TG 6.1% 201.4 / 4.3 66intergenic (+11/‑91)tsaB/PP_1508endopeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base G (2/2);  total (35/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 = 3.88e-01

AACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATCGCA  >  NC_002947/1713594‑1713888
                                                                                                                                                     |                                                                                                                                                  
aaCCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTtcgt                                                                                                                                                    <  2:58699/149‑1 (MQ=255)
  cccGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTtcgtcg                                                                                                                                                  <  7:275155/149‑1 (MQ=255)
   ccGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGa                                                                                                                                                 <  6:151902/149‑1 (MQ=255)
          ccTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCg                                                                                                                                           <  6:222684/148‑1 (MQ=255)
           cTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCCCGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGcc                                                                                                                                         >  5:327433/1‑149 (MQ=255)
           cTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGcc                                                                                                                                         >  7:188943/1‑149 (MQ=255)
                ccAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGGCGATAATCGACGAaca                                                                                                                                    <  3:293406/149‑1 (MQ=255)
                 cAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAacac                                                                                                                                   >  6:98424/1‑149 (MQ=255)
                  aaTGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACg                                                                                                                                  <  6:277919/149‑1 (MQ=255)
                  aaTGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACg                                                                                                                                  >  4:415697/1‑149 (MQ=255)
                     ccgctGGACCATCTCAGCCGCGCCCCCGTTGCCTGCGCCCGTGGAGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTa                                                                                                                               <  4:336408/146‑1 (MQ=255)
                      ctctGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGGCGATAATCGACGAacagcact                                                                                                                               <  7:286364/148‑6 (MQ=255)
                       tctGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTtcgt                                                                                                                                                    <  6:424365/126‑1 (MQ=255)
                       tctGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTtcgt                                                                                                                                                    >  5:424365/1‑126 (MQ=255)
                       tctGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAaca                                                                                                                                    <  7:356809/142‑1 (MQ=255)
                       tctGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAaca                                                                                                                                    >  8:356809/1‑142 (MQ=255)
                       tctGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTaa                                                                                                                              >  7:362926/1‑148 (MQ=255)
                       tatGGACATTATCAGCCTGGCCACCTTTTCCTGGGCCCGTGGCGAGGCAATCCTCGCCGAGCAGGCGCAACACGTTTATCTGCGCCATAATGTAACCACCCCCAAGAAGCACAGAGTGCTGTtagtag                                                                                                                                                  >  5:404997/3‑126 (MQ=255)
                           gACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAAcc                                                                                                                          <  3:397180/148‑1 (MQ=255)
                            aCATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCtt                                                                                                                        <  5:380038/149‑1 (MQ=255)
                             cATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCttt                                                                                                                       <  1:215339/149‑1 (MQ=255)
                                 ctcAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTgcg                                                                                                                   >  3:341110/1‑149 (MQ=255)
                                 ctcAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTgcg                                                                                                                   >  1:293553/1‑149 (MQ=255)
                                  tcAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTgcgc                                                                                                                  >  1:298359/1‑149 (MQ=255)
                                  tcAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTgcgc                                                                                                                  >  1:315808/1‑149 (MQ=255)
                                   cAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTgcgc                                                                                                                  >  4:116125/1‑148 (MQ=255)
                                          gCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCtgtg                                                                                                          <  4:218181/149‑1 (MQ=255)
                                             aCCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTc                                                                                                       <  2:415084/149‑1 (MQ=255)
                                                  ggCCTGGGGGCGTGGCGAGGGAAGCGTGGCCGAGCAGGCGCATGAGTTTTATCTGCGTGATAATGTAGCCTCTCCGATGATGCTCTGAGTGGTGTTCGTTGATTATCGGCATTCTGTCTATAACCTTTTGCGCGATCTGTGGTCCa                                                                                                     <  4:54417/145‑1 (MQ=255)
                                                   gCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGtt                                                                                                  <  3:116125/148‑1 (MQ=255)
                                                    ccTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTAt                                                                                                <  2:298359/149‑1 (MQ=255)
                                                           ccGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCaa                                                                                         <  3:91455/149‑1 (MQ=255)
                                                            cGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAg                                                                                        <  4:341110/149‑1 (MQ=255)
                                                              tGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCa                                                                                      >  1:144040/1‑149 (MQ=255)
                                                                  gAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAg                                                                                  <  2:290456/149‑1 (MQ=255)
                                                                   aGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGc                                                                                 >  5:134587/1‑149 (MQ=255)
                                                                       aaTCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACg                                                                             >  6:89590/1‑149 (MQ=255)
                                                                        aTCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACgg                                                                            <  3:415697/149‑1 (MQ=255)
                                                                              gCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAActct                                                                      >  6:163625/1‑149 (MQ=255)
                                                                              gCCGAGCAGGCGCAACCGGTGTATCTGCGCGATAATGTAGCCACGCCCAATAAGCACTGAGTGCTGGTCGTCTATTATCGCCGAACACGATAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACAGAActc                                                                       <  1:225932/148‑1 (MQ=255)
                                                                               ccGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAActct                                                                      >  1:308084/1‑149 (MQ=255)
                                                                               ccGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAActct                                                                      >  7:412531/1‑148 (MQ=255)
                                                                                  aGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGAtg                                                                  >  3:358438/1‑149 (MQ=255)
                                                                                       gcgcAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGAtgctgc                                                              >  6:282124/1‑148 (MQ=255)
                                                                                       gcgcAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGAGATCTTTGGTCCAGTTATCACTCAAGCATTAGCGAAGGAACTCTGAtgctgc                                                              >  6:323106/1‑148 (MQ=255)
                                                                                         ggacccGGTTTATCTGCGCGATAATGTAGCGCGGCCCAAGACGCACTGAGTGCTGTACCTCGATTATCGCCGCACCCGCTAAAACCTGTTGCGCGCTCTGTGGGCCAGTTATCACTCAGGCATTACCTACGCACCTCGGATGCTGCTa                                                            >  6:111979/5‑148 (MQ=255)
                                                                                             ccGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTg                                                       <  6:134587/149‑1 (MQ=255)
                                                                                               ggTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGcc                                                     <  2:144040/149‑1 (MQ=255)
                                                                                                gTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCca                                                    >  5:310207/1‑149 (MQ=255)
                                                                                                 tttATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCca                                                    >  5:152816/1‑148 (MQ=255)
                                                                                                    aTCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCcatcat                                                >  5:300/1‑149 (MQ=255)
                                                                                                     tCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCAtt                                               >  2:367740/1‑149 (MQ=255)
                                                                                                         cgcgATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTg                                              <  6:327433/146‑1 (MQ=255)
                                                                                                                   tAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTc                                 <  4:202328/149‑1 (MQ=255)
                                                                                                                     gCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTcat                               <  2:293553/149‑1 (MQ=255)
                                                                                                                             ccAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTcatgccatgc                        <  5:163625/149‑1 (MQ=255)
                                                                                                                             ccAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCAATCATTGGTCCTGAGTGATcatgccat                          >  4:384674/1‑147 (MQ=255)
                                                                                                                               aagaagCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACAAATTGCGCGATATGTGGTCCAGTTATCACTCAAGCATTAGAGACGGAACTCTGATGCTGCTAAATTGCCATCATTGCTCCTGAGTGCTCATGCCATGCg                       >  2:273110/1‑148 (MQ=255)
                                                                                                                               aagaagCACTGAGTGCTGTTCGTCGATTATCGACGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTCGTCCTGAGTGCTCATGCCATGCg                       >  2:267007/1‑148 (MQ=255)
                                                                                                                                    gCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGa                 >  4:336999/1‑149 (MQ=255)
                                                                                                                                      aCTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACg               >  6:307651/1‑149 (MQ=255)
                                                                                                                                         gAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTcatgccat                          <  6:337140/135‑1 (MQ=255)
                                                                                                                                         gAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTcatgccat                          >  5:337140/1‑135 (MQ=255)
                                                                                                                                          aGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTc                                 >  6:283219/1‑127 (MQ=255)
                                                                                                                                          aGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTc                                 <  5:283219/127‑1 (MQ=255)
                                                                                                                                           gTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTtc          <  5:89590/149‑1 (MQ=255)
                                                                                                                                              gtGTTCGGCGATAATCGACGAACACGCTAAACCCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATCCTGCTAAATTGCCATCATTGGTCCTGAGAGCACCTGCCATGCGTATCGAAGGATTATCa       >  5:434402/2‑149 (MQ=255)
                                                                                                                                              gtGTTCGGCGATAATCGACGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCa       >  1:27635/2‑149 (MQ=255)
                                                                                                                                              cTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCa       <  5:282124/149‑1 (MQ=255)
                                                                                                                                               tGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCAt      <  3:336999/149‑1 (MQ=255)
                                                                                                                                                gTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATc     >  2:82739/1‑149 (MQ=255)
                                                                                                                                                    gtcgATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATCGCa  <  5:307651/148‑1 (MQ=255)
                                                                                                                                                     tcgATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTc                                           <  1:325161/106‑1 (MQ=255)
                                                                                                                                                     tcgATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTc                                           >  2:325161/1‑106 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATCGCA  >  NC_002947/1713594‑1713888

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: