Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 247,385 G→T 6.6% S225I (AGC→ATC)  yafL → putative lipoprotein and C40 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913247,3850GT 6.6% 261.9 / 7.2 76S225I (AGC→ATC) yafLputative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (33/38);  new base T (2/3);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTA  >  NC_000913/247242‑247531
                                                                                                                                               |                                                                                                                                                  
gATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAt                                                                                                                                               >  2:737602/1‑149 (MQ=255)
 aTTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATCTGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAtt                                                                                                                                              <  4:694552/149‑1 (MQ=255)
 aTTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAtt                                                                                                                                              >  4:116056/1‑149 (MQ=255)
 aTTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAtt                                                                                                                                              >  4:614947/1‑149 (MQ=255)
     tgGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGcc                                                                                                                                          <  3:476334/149‑1 (MQ=255)
           acaacGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaatgg                                                                                                                                    <  1:174842/149‑4 (MQ=255)
             aacGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATAGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGcc                                                                                                                                                  <  3:664049/133‑1 (MQ=255)
             aacGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATAGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGcc                                                                                                                                                  >  4:664049/1‑133 (MQ=255)
              acGCCCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGACGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTc                                                                                                                                 <  3:614947/149‑1 (MQ=255)
               cGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCt                                                                                                                                <  4:42029/149‑1 (MQ=255)
                           gggAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGc                                                                                                                                           >  3:724580/1‑126 (MQ=255)
                           gggAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGc                                                                                                                                           <  4:724580/126‑1 (MQ=255)
                           gggAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtt                                                                                                                    <  1:753248/149‑1 (MQ=255)
                            ggAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAttt                                                                                                                   >  1:443500/1‑149 (MQ=255)
                               gaTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtttttt                                                                                                                <  2:443500/149‑1 (MQ=255)
                               gaTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtttttt                                                                                                                <  1:814624/149‑1 (MQ=255)
                               gaTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtttttt                                                                                                                <  2:352039/149‑1 (MQ=255)
                                aTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTg                                                                                                               <  2:231939/149‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                           >  4:437988/1‑136 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                           >  2:348989/1‑136 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                           <  1:348989/136‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                           <  3:437988/136‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTgg                                                                                                              <  3:81317/149‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTgg                                                                                                              <  4:670076/149‑1 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                       >  3:296120/1‑110 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                       <  2:114686/110‑1 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                       <  4:296120/110‑1 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                       >  1:114686/1‑110 (MQ=255)
                                                tatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAAc                                                                                                                                      <  4:234369/110‑1 (MQ=255)
                                                tatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAAc                                                                                                                                      >  3:234369/1‑110 (MQ=255)
                                                 atCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGAtttt                                                                                              <  1:145493/149‑1 (MQ=255)
                                                 atCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGAtttt                                                                                              <  3:399548/149‑1 (MQ=255)
                                                   ccACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGa                                                                                            >  4:19409/1‑149 (MQ=255)
                                                     acaGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACg                                                                                          >  4:66727/1‑149 (MQ=255)
                                                        gCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGaa                                                                                       <  3:19409/149‑1 (MQ=255)
                                                            cgAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAgaga                                                                                   <  3:667969/149‑1 (MQ=255)
                                                                    ccGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTg                                                                           <  4:718306/149‑1 (MQ=255)
                                                                      taTCCTATGGGCGTGTATTTGGGCGATGGGCCTTTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCTGGACCATCTTTTGGGCGCTCGCAGGAGTTTAACGGAAGAGCCGATTTTGTa                                                                         >  4:752459/2‑149 (MQ=255)
                                                                      gATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTa                                                                         <  3:116056/149‑1 (MQ=255)
                                                                      gATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTa                                                                         <  2:145194/149‑1 (MQ=255)
                                                                          atatGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGAc                                                                     <  4:787326/149‑1 (MQ=255)
                                                                             tGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAtt                                                                                                                                              >  4:765243/1‑73 (MQ=255)
                                                                             tGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAtt                                                                                                                                              <  3:765243/73‑1 (MQ=255)
                                                                              gggCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGa                                                                                                  >  1:219345/1‑116 (MQ=255)
                                                                              gggCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGa                                                                                                  <  2:219345/116‑1 (MQ=255)
                                                                              gggCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGAt                                                                 >  2:465248/1‑149 (MQ=255)
                                                                                gCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATaa                                                               <  1:745905/149‑1 (MQ=255)
                                                                                gCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCAAGTAACGGCGAAACCCTTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAATTTTTGGGCGGGCGCAGGATTTTGAAGGAAGAGACGATTTTGTAGGACGGATaa                                                               <  3:66727/149‑1 (MQ=255)
                                                                                       tttGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGttt                                                        >  1:512608/1‑149 (MQ=255)
                                                                                          gggCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTAcg                                                     >  1:473196/1‑149 (MQ=255)
                                                                                                  ggCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATcc                                             >  4:115320/1‑149 (MQ=255)
                                                                                                   gCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTAc                                                      >  3:36574/1‑139 (MQ=255)
                                                                                                   gCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTAc                                                      <  4:36574/139‑1 (MQ=255)
                                                                                                   gCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCg                                            >  2:520640/1‑149 (MQ=255)
                                                                                                          aTCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTg                                     >  4:314573/1‑149 (MQ=255)
                                                                                                               gTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGc                                >  3:246301/1‑149 (MQ=255)
                                                                                                                tCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCa                               >  1:560789/1‑149 (MQ=255)
                                                                                                                 cGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAg                              <  1:737602/149‑1 (MQ=255)
                                                                                                                   ccACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGaa                                                                                       <  3:699374/90‑1 (MQ=255)
                                                                                                                   ccACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGaa                                                                                       >  4:699374/1‑90 (MQ=255)
                                                                                                                    cACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGAt                                                                               >  1:550206/1‑97 (MQ=255)
                                                                                                                    cACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGAt                                                                               <  2:550206/97‑1 (MQ=255)
                                                                                                                     aCGTACCGGCGAAACCAGTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAAGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGc                          <  4:676309/149‑1 (MQ=255)
                                                                                                                        tACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTg                       <  4:246301/149‑1 (MQ=255)
                                                                                                                          ccGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGAt                     <  1:465248/149‑1 (MQ=255)
                                                                                                                               gAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCa                               >  4:296466/1‑134 (MQ=255)
                                                                                                                               gAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCa                               <  3:296466/134‑1 (MQ=255)
                                                                                                                                  aCCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCt             >  4:580890/1‑149 (MQ=255)
                                                                                                                                  aCCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCt             >  4:139382/1‑149 (MQ=255)
                                                                                                                                     aTTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc          >  3:63989/1‑149 (MQ=255)
                                                                                                                                     aTTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc          >  1:162554/1‑149 (MQ=255)
                                                                                                                                           aTAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCt    >  3:399218/1‑149 (MQ=255)
                                                                                                                                             aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa  <  3:290036/149‑1 (MQ=255)
                                                                                                                                             aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa  >  2:66202/1‑149 (MQ=255)
                                                                                                                                             aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa  <  3:314573/149‑1 (MQ=255)
                                                                                                                                             aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa  >  4:290036/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTA  >  NC_000913/247242‑247531

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: