Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,589 G→C 15.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,5890GC15.3% 355.1 / 56.7 139A97P (GCA→CCA) yafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (55/62);  new base C (12/9);  total (68/71)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-02

Only 100 of 153 total aligned reads displayed.

CAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCA  >  NC_000913/252491‑252688
                                                                                                  |                                                                                                   
cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaagc                                                                                                    >  2:2034005/1‑100 (MQ=255)
cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                     >  1:602190/1‑99 (MQ=255)
 aaaTGTGAATAATCCATTTCCGCCACAGTTGTGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGTGCGTTTGATGagcaagca                                                                                                   <  1:3799257/100‑1 (MQ=255)
 aaaTGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaagca                                                                                                   <  2:1184356/100‑1 (MQ=255)
 aaaTGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaagc                                                                                                    <  1:400591/99‑1 (MQ=255)
    tGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATgg                                                                                                >  2:437289/1‑100 (MQ=255)
         ataatCCATTTCCGCCACAGAAACGACAATTCAACAGAACGAATGACGAAGCGCATTTTGTATATTGACAGGAGGCGTTTGATGAGCACACAAGGTTGCt                                                                                           >  2:1047697/1‑100 (MQ=255)
            atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                        <  1:2329174/100‑1 (MQ=255)
            atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAt                                                                                        <  2:3219459/100‑1 (MQ=255)
            atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAt                                                                                        <  2:227009/100‑1 (MQ=255)
             tCCATTTCCGCCACAGTTGCGCCAATTCAGTAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAt                                                                                        <  2:1984705/99‑1 (MQ=255)
             tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                       <  1:2941973/100‑2 (MQ=255)
             tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                       <  1:93235/100‑2 (MQ=255)
              ccATTTCCGCCACAGTTGCGCCAATTCAGTAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTg                                                                                      <  1:84628/100‑1 (MQ=255)
              ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTg                                                                                      <  1:2026104/100‑1 (MQ=255)
               cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                     <  2:3389026/100‑1 (MQ=255)
               cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                     >  2:1418943/1‑100 (MQ=255)
               cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                     <  1:3488757/100‑1 (MQ=255)
                 tttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGTGCGTTTGATGAGCAACCATGCTTGCTCAtcaaac                                                                                   <  2:1836020/100‑6 (MQ=38)
                 tttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaac                                                                                   <  1:2578661/100‑6 (MQ=39)
                     cGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                               <  1:706289/100‑1 (MQ=255)
                      gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTg                                                                              >  1:1965151/1‑100 (MQ=255)
                        cacaGTTGCGCCAATTCAGCAGAACGAATTACGAAGCGTATTTGGTATATTGTCAGGGGGCGGTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaa                                                                            <  1:3036714/100‑1 (MQ=255)
                        cacaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaa                                                                            <  1:1112323/100‑1 (MQ=255)
                        aacaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaa                                                                            >  1:2258546/2‑100 (MQ=255)
                         acaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaaa                                                                           <  2:1466169/100‑1 (MQ=255)
                               gcgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                     <  2:3730698/100‑1 (MQ=255)
                                cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGa                                                                      >  2:2717112/1‑98 (MQ=255)
                                cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                    <  2:2550002/100‑1 (MQ=255)
                                 gcCAATTCATTTGAGCGAATTACGAAACGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCTTGGTTTCTCATTGCCATTCTGAAACCTGaa                                                                     <  2:2258340/98‑1 (MQ=255)
                                 gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                   <  1:1836946/100‑1 (MQ=255)
                                 gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                   <  2:3880735/100‑1 (MQ=255)
                                 gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTAAGGGGGCGTTTGATGAGAAAGCATGGTTGATCATTGCCATTCTGAAACCTGAAcc                                                                   >  2:524902/1‑100 (MQ=255)
                                   cAATTCAGCAGAACGAATGACGAAGAGCAATTGGTATATTGTCAAGGGGCGCTTGATGAGCAAGAATAGTTGCTCATTGCCATTGTGAAACCTGaaggtc                                                                 >  2:2813865/1‑96 (MQ=255)
                                       tCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaa                                                             <  2:1683870/100‑1 (MQ=255)
                                        cagcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                            >  1:1112533/1‑100 (MQ=255)
                                        cagcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                            >  1:1246501/1‑100 (MQ=255)
                                        cagcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                            >  2:3212207/1‑100 (MQ=255)
                                          gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATAGTTGCACATTGCCATTCTGAAACCTGAACCTCATAAACt                                                          >  1:745342/1‑100 (MQ=255)
                                            agAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                        <  1:1416777/100‑1 (MQ=255)
                                             gAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                       <  1:1355550/100‑1 (MQ=255)
                                             gAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                       <  2:3280091/100‑1 (MQ=255)
                                             gAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                       >  1:3302511/1‑100 (MQ=255)
                                               aCGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTc                                                     >  1:3074249/1‑100 (MQ=255)
                                                cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                    >  2:3772656/1‑100 (MQ=255)
                                                cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                    <  1:1599912/100‑1 (MQ=255)
                                                cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                    <  1:2343556/100‑1 (MQ=255)
                                                 gAATGACGAAGCGCATTTGGTATATTGTCAGGTGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                    <  2:3928951/99‑1 (MQ=255)
                                                 gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                    <  2:2206739/99‑1 (MQ=255)
                                                  aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                  <  1:2560451/100‑1 (MQ=255)
                                                  aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                  <  1:2780915/100‑1 (MQ=255)
                                                  aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTAATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                  >  1:351005/1‑100 (MQ=255)
                                                    tGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGAt                                                <  2:2156083/100‑1 (MQ=255)
                                                     gACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATa                                               <  2:227616/100‑1 (MQ=255)
                                                     gACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATa                                               >  2:3099111/1‑100 (MQ=255)
                                                       cGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaac                                             >  1:579886/1‑100 (MQ=255)
                                                         aaGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGACTCGAGATaacaa                                           >  1:838905/1‑100 (MQ=255)
                                                          aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaacaac                                          <  1:1066885/100‑1 (MQ=255)
                                                           gcgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAAcc                                         <  1:2881571/100‑1 (MQ=255)
                                                           gcgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAAcc                                         <  1:2881692/100‑1 (MQ=255)
                                                               aTTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAAt                                     >  1:102844/1‑100 (MQ=255)
                                                                  tGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCa                                  >  2:3781752/1‑100 (MQ=255)
                                                                    gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                                >  1:1800402/1‑100 (MQ=255)
                                                                    gTATATTGTCAGGGGGAGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAACTGAAACTCATAAACTGGCTCGAGATAAAAAACAAATgaa                                  >  2:948385/1‑96 (MQ=255)
                                                                     tataTTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGATATAACAACCAAATGCATaa                               >  2:3020904/1‑100 (MQ=255)
                                                                     tataTTGTCAGGAGGCGTATGATGAGCAAGCATGGTTGCTCATTACCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaa                               >  2:1252295/1‑100 (MQ=255)
                                                                        atTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAt                            >  2:3652390/1‑100 (MQ=255)
                                                                          tGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                          <  2:3432008/100‑1 (MQ=255)
                                                                           gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                         <  2:511116/100‑1 (MQ=255)
                                                                             cAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGa                       >  1:3419104/1‑100 (MQ=255)
                                                                               gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaa                     >  1:370357/1‑100 (MQ=255)
                                                                                ggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaaa                    >  1:1685329/1‑100 (MQ=255)
                                                                                 gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGACTCGAGATAACAACCAAATGCATAAAATTGGGaaa                     >  2:3918188/1‑98 (MQ=255)
                                                                                  ggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATg                  <  2:1995278/100‑1 (MQ=255)
                                                                                   gCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATgg                 <  2:3419065/100‑1 (MQ=255)
                                                                                     gTTTGATGAGCAAGCATGGTTGCTCATTGTCATTTTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa               <  1:16388/100‑1 (MQ=255)
                                                                                     gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa               <  2:799320/100‑1 (MQ=255)
                                                                                      ggcaatgaGCAACCATGCTTGCTCATTACCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTTGGAAAATGGCa               >  2:1212423/5‑99 (MQ=38)
                                                                                       gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa               >  2:1896235/4‑98 (MQ=255)
                                                                                       gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGa             >  2:109559/4‑100 (MQ=255)
                                                                                        tgatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAAAAACCAAATGCCTACCATTGGGAACTTGGCAgga            >  1:1687384/1‑98 (MQ=255)
                                                                                        caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  1:3318568/3‑100 (MQ=255)
                                                                                         gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg           >  1:990694/1‑100 (MQ=255)
                                                                                         gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg           <  1:3568405/100‑1 (MQ=255)
                                                                                         gatgaGCAAGCATGATTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGCAATAACAACCAAATGCATAAAATTGGGAAAATCGCAGaa            >  2:369353/1‑99 (MQ=255)
                                                                                         aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg           >  1:2246681/2‑100 (MQ=255)
                                                                                          atgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc          >  1:2121149/1‑100 (MQ=255)
                                                                                           tgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCg         <  2:1696623/100‑1 (MQ=255)
                                                                                           tgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCg         >  2:3668191/1‑100 (MQ=255)
                                                                                           tgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCg         >  1:2945631/1‑100 (MQ=255)
                                                                                           tgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCg         >  1:3685133/1‑100 (MQ=255)
                                                                                            gaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGCAAATGGCAGAAGCGt        >  2:3708538/1‑100 (MQ=255)
                                                                                            gaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGt        >  1:3625870/1‑100 (MQ=255)
                                                                                             agcacgcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAAAGGCAGAAGCGtt       >  2:3431212/1‑100 (MQ=255)
                                                                                             agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGtt       >  2:3293449/1‑100 (MQ=255)
                                                                                             agcaagcaTGGTTGCTCATTGCCATTATGAAACCTAAACCTCATAAACTGGCTCGGGATAACAACAACAGACATAAAATTTGAAAAATGGCAAAAGCggt       >  2:20513/1‑98 (MQ=255)
                                                                                               caagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTc     >  2:764811/1‑100 (MQ=255)
                                                                                                 agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc   >  2:2513075/1‑100 (MQ=255)
                                                                                                  gcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa  <  1:2783138/100‑1 (MQ=255)
                                                                                                  gcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa  <  1:1202425/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCA  >  NC_000913/252491‑252688

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: