Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,589 G→C 7.3% A97P (GCA→CCA)  yafO → mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913252,5890GC 7.3% 480.6 / 24.5 137A97P (GCA→CCA) yafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (59/68);  new base C (6/4);  total (65/72)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-02

Only 100 of 143 total aligned reads displayed.

CTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATC  >  NC_000913/252441‑252736
                                                                                                                                                    |                                                                                                                                                   
cTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                                                                     >  3:761164/1‑149 (MQ=255)
   tACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCAt                                                                                                                                                  <  3:149908/149‑1 (MQ=255)
   tACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCAt                                                                                                                                                  <  1:963833/149‑1 (MQ=255)
    aCCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTGGATGAGCAAGCATg                                                                                                                                                 <  2:608306/149‑1 (MQ=255)
     ccTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATgg                                                                                                                                                <  4:198446/149‑1 (MQ=255)
           cATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTc                                                                                                                                          <  4:14182/149‑1 (MQ=255)
            aTTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                                                                         >  1:275728/1‑149 (MQ=255)
            aTTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                                                                         <  1:379802/149‑1 (MQ=255)
              tAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAtt                                                                                                                                       <  3:383468/149‑1 (MQ=255)
              tAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAtt                                                                                                                                       <  4:632851/149‑1 (MQ=255)
               aaTCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTg                                                                                                                                      >  4:690689/1‑149 (MQ=255)
                 tCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                    >  4:343290/1‑149 (MQ=255)
                 tCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                    >  3:1049840/1‑149 (MQ=255)
                 tCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                    >  4:607941/1‑149 (MQ=255)
                  cAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCa                                                                                                                                   <  1:557162/149‑1 (MQ=255)
                   aaaTTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCACTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAt                                                                                                                                  >  3:480338/1‑149 (MQ=255)
                   aaaTTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAt                                                                                                                                  <  3:247455/149‑1 (MQ=255)
                   aaaTTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaacg                                                                                                                                  <  4:214582/149‑7 (MQ=255)
                      tttGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCTCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                                                               >  2:692375/1‑149 (MQ=255)
                      tttGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCGTGGTTGCTCATTGCCATTCt                                                                                                                               >  3:887669/1‑149 (MQ=255)
                      tttGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                                                               >  2:896999/1‑149 (MQ=255)
                          aGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaaa                                                                                                                           <  4:692962/149‑1 (MQ=255)
                           gCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAc                                                                                                                          <  1:372100/149‑1 (MQ=255)
                                ttGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                                                                     <  3:316909/149‑1 (MQ=255)
                                  gCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                                                                   >  1:106242/1‑149 (MQ=255)
                                   cTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCt                                                                                                                  >  1:906133/1‑149 (MQ=255)
                                   cTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCt                                                                                                                  >  4:215754/1‑149 (MQ=255)
                                    tCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                                                                                 >  3:119875/1‑149 (MQ=255)
                                    tCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                                                                                 <  2:494394/149‑1 (MQ=255)
                                    tCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCa                                                                                                                                         <  4:748281/125‑1 (MQ=255)
                                    tCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCa                                                                                                                                         >  3:748281/1‑125 (MQ=255)
                                     cATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                                                                >  2:234954/1‑149 (MQ=255)
                                          tCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAAc                                                                                                           <  3:874863/149‑1 (MQ=255)
                                           cATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                                                                          <  1:502068/149‑1 (MQ=255)
                                           cATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                                                                          <  2:87815/149‑1 (MQ=255)
                                            aTCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                                                                         <  4:1009278/149‑1 (MQ=255)
                                            aTCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                                                                         <  1:1044474/149‑1 (MQ=255)
                                            aTCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                                                                         <  1:613720/149‑1 (MQ=255)
                                             tCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                                                                        <  1:896999/149‑1 (MQ=255)
                                              cTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                                                                       <  2:1056959/149‑1 (MQ=255)
                                               tGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCt                                                                                                      >  4:819972/1‑149 (MQ=255)
                                                ggCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTc                                                                                                     <  4:155763/149‑1 (MQ=255)
                                                  cAAATGAGAATCATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCga                                                                                                   <  1:275859/149‑1 (MQ=255)
                                                     aTGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGAt                                                                                                <  4:661176/149‑1 (MQ=255)
                                                      tGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATa                                                                                               >  2:963870/1‑149 (MQ=255)
                                                         gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCa                                                                                                                                         >  2:635717/1‑104 (MQ=255)
                                                         gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCa                                                                                                                                         <  1:635717/104‑1 (MQ=255)
                                                          aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaacaa                                                                                           >  3:562685/1‑149 (MQ=255)
                                                              atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                                                                                <  3:726564/108‑1 (MQ=255)
                                                              atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaa                                                                                       >  4:872554/1‑149 (MQ=255)
                                                               tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                                                               <  3:3270/108‑1 (MQ=255)
                                                               tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                                                               >  4:3270/1‑108 (MQ=255)
                                                               tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaaa                                                                                      >  2:107688/1‑149 (MQ=255)
                                                                  aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                                                                   >  4:599651/1‑149 (MQ=255)
                                                                   tttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCa                                                                                  >  4:351142/1‑149 (MQ=255)
                                                                        gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaaaa                                                                             >  2:855581/1‑149 (MQ=255)
                                                                         ccACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATCACAACCAAATGCATAAAAt                                                                            >  4:736432/1‑149 (MQ=255)
                                                                           acaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                                                                          <  4:558502/149‑1 (MQ=255)
                                                                                 gcgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaaa                                                                    <  3:215754/149‑1 (MQ=255)
                                                                                  cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAAt                                                                   <  3:447468/149‑1 (MQ=255)
                                                                                   gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATg                                                                  >  1:687788/1‑149 (MQ=255)
                                                                                    ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATgg                                                                 <  3:599651/149‑1 (MQ=255)
                                                                                    ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATgg                                                                 <  4:887669/149‑1 (MQ=255)
                                                                                     cAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                                                                <  3:448655/149‑1 (MQ=255)
                                                                                     cAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                                                                <  4:971530/149‑1 (MQ=255)
                                                                                            gcagAACGATTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                                                                                 >  4:930396/1‑93 (MQ=255)
                                                                                            gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                                                                                 <  3:930396/93‑1 (MQ=255)
                                                                                             cagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGt                                                        <  4:562685/149‑1 (MQ=255)
                                                                                                   gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa                                                  >  2:137549/1‑149 (MQ=255)
                                                                                                    aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCAt                                                 >  3:261271/1‑149 (MQ=255)
                                                                                                      tGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGc                                               <  4:1049840/149‑1 (MQ=255)
                                                                                                             gcgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                                                                       >  1:764460/1‑86 (MQ=255)
                                                                                                             gcgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                                                                       <  2:764460/86‑1 (MQ=255)
                                                                                                               gcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGaa                                      <  3:343290/149‑1 (MQ=255)
                                                                                                                      gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTatat                               >  2:937358/1‑149 (MQ=255)
                                                                                                                       tataTTGTCAGGGGACGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtatt                              >  4:237314/1‑149 (MQ=255)
                                                                                                                             gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAAt                        <  4:43692/149‑1 (MQ=255)
                                                                                                                              tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                                                                        >  4:1042777/1‑100 (MQ=255)
                                                                                                                              tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                                                                        <  3:1042777/100‑1 (MQ=255)
                                                                                                                              tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATa                       >  2:642340/1‑149 (MQ=255)
                                                                                                                               cAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATaa                      <  3:690689/149‑1 (MQ=255)
                                                                                                                                 gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGaa                                      >  2:611129/1‑131 (MQ=255)
                                                                                                                                  ggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACAt                   >  3:491307/1‑149 (MQ=255)
                                                                                                                                   gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATa                  <  4:976507/149‑1 (MQ=255)
                                                                                                                                   gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATa                  <  3:872554/149‑1 (MQ=255)
                                                                                                                                   gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATa                  <  2:259603/149‑1 (MQ=255)
                                                                                                                                    ggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATAc                 >  1:972633/1‑149 (MQ=255)
                                                                                                                                       tggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaaa              >  1:76981/6‑149 (MQ=255)
                                                                                                                                       gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaaa              >  2:523800/1‑149 (MQ=255)
                                                                                                                                           gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg          >  4:577037/1‑149 (MQ=255)
                                                                                                                                           gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg          >  1:642301/1‑149 (MQ=255)
                                                                                                                                           aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaaaa                                                                             >  4:418715/2‑82 (MQ=255)
                                                                                                                                           aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaaaa                                                                             <  3:418715/81‑1 (MQ=255)
                                                                                                                                           aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg          >  4:23596/2‑149 (MQ=255)
                                                                                                                                              gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGaaa       >  1:606079/1‑149 (MQ=255)
                                                                                                                                                gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACt     >  3:565887/1‑149 (MQ=255)
                                                                                                                                                gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACt     >  1:1053583/1‑149 (MQ=255)
                                                                                                                                                  aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTAt   >  2:614613/1‑149 (MQ=255)
                                                                                                                                                  aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTAt   >  1:996942/1‑149 (MQ=255)
                                                                                                                                                   agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATc  <  4:731561/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATC  >  NC_000913/252441‑252736

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