Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 247,389 A→C 6.7% R226R (CGA→CGC) ‡ 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,3890AC 6.7% 250.9 / 10.7 75R226R (CGA→CGC) ‡yafLputative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base A (33/37);  new base C (2/3);  total (35/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 = 4.32e-01

GATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGC  >  NC_000913/247242‑247537
                                                                                                                                                   |                                                                                                                                                    
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)
              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                                                                                                                                 <  1:348989/136‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                                 <  3:437988/136‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                                 >  4:437988/1‑136 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                                 >  2:348989/1‑136 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTgg                                                                                                                    <  3:81317/149‑1 (MQ=255)
                                 tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTgg                                                                                                                    <  4:670076/149‑1 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                             <  4:296120/110‑1 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                             >  3:296120/1‑110 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                             >  1:114686/1‑110 (MQ=255)
                                               atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                             <  2:114686/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)
                                                                             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: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        >  4:290036/1‑149 (MQ=255)
                                                                                                                                             aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa        <  3:314573/149‑1 (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)
                                                                                                                                                   |                                                                                                                                                    
GATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGC  >  NC_000913/247242‑247537

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: