Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 247,381 A→C 6.7% I224L (ATA→CTA)  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,3810AC 6.7% 255.5 / 9.6 74I224L (ATA→CTA) yafLputative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base A (30/39);  new base C (2/3);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

CCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCT  >  NC_000913/247233‑247529
                                                                                                                                                    |                                                                                                                                                    
ccGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCAGTCGGa                                                                                                                                                      <  2:32176/149‑1 (MQ=255)
 cGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGAt                                                                                                                                                     <  2:624284/149‑1 (MQ=255)
  gggCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATa                                                                                                                                                    <  4:166755/149‑1 (MQ=255)
   ggCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATaa                                                                                                                                                   <  3:442154/149‑1 (MQ=255)
         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                                                                                                              <  1:814624/149‑1 (MQ=255)
                                        gaTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtttttt                                                                                                              <  2:443500/149‑1 (MQ=255)
                                        gaTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCAtttttt                                                                                                              <  2:352039/149‑1 (MQ=255)
                                         aTTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTg                                                                                                             <  2:231939/149‑1 (MQ=255)
                                          tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                         >  2:348989/1‑136 (MQ=255)
                                          tttGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAg                                                                                                                         >  4:437988/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                                                                                                                                     >  1:114686/1‑110 (MQ=255)
                                                        atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                     >  3:296120/1‑110 (MQ=255)
                                                        atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                     <  2:114686/110‑1 (MQ=255)
                                                        atatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGCTAATCGGCTTATCCGaa                                                                                                                                     <  4:296120/110‑1 (MQ=255)
                                                         tatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAAc                                                                                                                                    >  3:234369/1‑110 (MQ=255)
                                                         tatCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAAc                                                                                                                                    <  4:234369/110‑1 (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                                                                       <  2:145194/149‑1 (MQ=255)
                                                                               gATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTa                                                                       <  3:116056/149‑1 (MQ=255)
                                                                                   atatGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGAc                                                                   <  4:787326/149‑1 (MQ=255)
                                                                                      tGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAtt                                                                                                                                            <  3:765243/73‑1 (MQ=255)
                                                                                      tGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAtt                                                                                                                                            >  4:765243/1‑73 (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: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)
                                                                                                                                                    |                                                                                                                                                    
CCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCACAGCCGCGAGATAGCCGATCATATGGGCGTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCT  >  NC_000913/247233‑247529

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: