Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 247,393 G→T 6.9% A228S (GCC→TCC)  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,3930GT 6.9% 242.8 / 5.4 70A228S (GCC→TCC) yafLputative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (29/36);  new base T (2/3);  total (31/39)
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.98e-01

TGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGC  >  NC_000913/247247‑247537
                                                                                                                                                  |                                                                                                                                                
tgGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGcc                                                                                                                                                <  3:476334/149‑1 (MQ=255)
      acaacGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaatgg                                                                                                                                          <  1:174842/149‑4 (MQ=255)
         acGCCCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGACGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTc                                                                                                                                       <  3:614947/149‑1 (MQ=255)
          cGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCt                                                                                                                                      <  4:42029/149‑1 (MQ=255)
                      gggAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGc                                                                                                                                                 <  4:724580/126‑1 (MQ=255)
                      gggAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGc                                                                                                                                                 >  3:724580/1‑126 (MQ=255)
                      gggAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtt                                                                                                                          <  1:753248/149‑1 (MQ=255)
                       ggAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAttt                                                                                                                         >  1:443500/1‑149 (MQ=255)
                          gaTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtttttt                                                                                                                      <  1:814624/149‑1 (MQ=255)
                          gaTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtttttt                                                                                                                      <  2:352039/149‑1 (MQ=255)
                          gaTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtttttt                                                                                                                      <  2:443500/149‑1 (MQ=255)
                           aTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTg                                                                                                                     <  2:231939/149‑1 (MQ=255)
                            tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                                 >  2:348989/1‑136 (MQ=255)
                            tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                                 <  1:348989/136‑1 (MQ=255)
                            tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                                 <  3:437988/136‑1 (MQ=255)
                            tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                                 >  4:437988/1‑136 (MQ=255)
                            tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTgg                                                                                                                    <  4:670076/149‑1 (MQ=255)
                            tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTgg                                                                                                                    <  3:81317/149‑1 (MQ=255)
                                          atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                             <  2:114686/110‑1 (MQ=255)
                                          atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                             >  1:114686/1‑110 (MQ=255)
                                          atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                             >  3:296120/1‑110 (MQ=255)
                                          atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                             <  4:296120/110‑1 (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                                                                               <  3:116056/149‑1 (MQ=255)
                                                                 gATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTa                                                                               <  2:145194/149‑1 (MQ=255)
                                                                     atatGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGAc                                                                           <  4:787326/149‑1 (MQ=255)
                                                                         gggCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGa                                                                                                        <  2:219345/116‑1 (MQ=255)
                                                                         gggCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGa                                                                                                        >  1:219345/1‑116 (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                                                                                     >  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:139382/1‑149 (MQ=255)
                                                                                                                             aCCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCt                   >  4:580890/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:314573/149‑1 (MQ=255)
                                                                                                                                        aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa        >  4:290036/1‑149 (MQ=255)
                                                                                                                                        aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa        <  3:290036/149‑1 (MQ=255)
                                                                                                                                        aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa        >  2:66202/1‑149 (MQ=255)
                                                                                                                                           ccGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCa     >  1:667716/1‑149 (MQ=255)
                                                                                                                                              aTTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGc  <  4:93816/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGC  >  NC_000913/247247‑247537

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