Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 247,387 C→G 6.8% R226G (CGA→GGA) ‡ 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,3870CG 6.8% 259.8 / 9.3 74R226G (CGA→GGA) ‡yafLputative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base G (2/3);  total (34/40)
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.06e-01

GATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCA  >  NC_000913/247242‑247534
                                                                                                                                                 |                                                                                                                                                   
gATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAt                                                                                                                                                  >  2:737602/1‑149 (MQ=255)
 aTTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATCTGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAtt                                                                                                                                                 <  4:694552/149‑1 (MQ=255)
 aTTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAtt                                                                                                                                                 >  4:614947/1‑149 (MQ=255)
 aTTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAtt                                                                                                                                                 >  4:116056/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                                                                                                                   <  2:352039/149‑1 (MQ=255)
                               gaTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtttttt                                                                                                                   <  1:814624/149‑1 (MQ=255)
                                aTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTg                                                                                                                  <  2:231939/149‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                              <  3:437988/136‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                              >  2:348989/1‑136 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                              >  4:437988/1‑136 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                              <  1:348989/136‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTgg                                                                                                                 <  4:670076/149‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTgg                                                                                                                 <  3:81317/149‑1 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                          <  2:114686/110‑1 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                          <  4:296120/110‑1 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                          >  1:114686/1‑110 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                          >  3:296120/1‑110 (MQ=255)
                                                tatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAAc                                                                                                                                         <  4:234369/110‑1 (MQ=255)
                                                tatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAAc                                                                                                                                         >  3:234369/1‑110 (MQ=255)
                                                 atCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGAtttt                                                                                                 <  3:399548/149‑1 (MQ=255)
                                                 atCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGAtttt                                                                                                 <  1:145493/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                                                                            <  2:145194/149‑1 (MQ=255)
                                                                      gATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTa                                                                            <  3:116056/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                                                                                          >  4:699374/1‑90 (MQ=255)
                                                                                                                   ccACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGaa                                                                                          <  3:699374/90‑1 (MQ=255)
                                                                                                                    cACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGAt                                                                                  <  2:550206/97‑1 (MQ=255)
                                                                                                                    cACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGAt                                                                                  >  1:550206/1‑97 (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     >  4:290036/1‑149 (MQ=255)
                                                                                                                                             aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa     >  2:66202/1‑149 (MQ=255)
                                                                                                                                             aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa     <  3:290036/149‑1 (MQ=255)
                                                                                                                                             aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa     <  3:314573/149‑1 (MQ=255)
                                                                                                                                                ccGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCa  >  1:667716/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCA  >  NC_000913/247242‑247534

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: