Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,594 G→C 17.3% W98C (TGG→TGC 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,5940GC17.3% 338.5 / 59.9 133W98C (TGG→TGCyafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (59/50);  new base C (10/13);  total (69/64)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 136 total aligned reads displayed.

TGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGT  >  NC_000913/252495‑252693
                                                                                                   |                                                                                                   
tGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATgg                                                                                                     >  1:2787586/1‑100 (MQ=255)
  agaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGtt                                                                                                   >  2:3129659/1‑100 (MQ=255)
   gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                  >  1:1516830/1‑100 (MQ=255)
   gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                  >  2:3145727/1‑100 (MQ=255)
    aataatCCATTTCCGCCACAGTTGCGCAAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGc                                                                                                 <  2:2063488/100‑1 (MQ=255)
     ataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCCAGCATGGTTGCt                                                                                                <  1:2392637/100‑1 (MQ=255)
     ataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCt                                                                                                <  2:955856/100‑1 (MQ=255)
     ataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCt                                                                                                <  1:2538868/100‑1 (MQ=255)
        atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                             >  1:1446018/1‑100 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAtt                                                                                            >  1:3386919/1‑100 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  1:527331/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  1:2349235/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  1:2347579/100‑2 (MQ=255)
          ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtca                                                                                           <  1:2208090/100‑3 (MQ=255)
          ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTCTCAGGGGGCGTTTGATGAGCAACCTTGCTTGCTCAtca                                                                                           <  2:2244480/100‑3 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                          <  1:1703815/100‑1 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaa                                                                                          <  1:3204621/100‑4 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaa                                                                                          <  1:3069751/100‑4 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTCTCAGCCAGCGTGTGATGAGCAAGCATGGTTGCTCATTg                                                                                           >  2:2560326/1‑99 (MQ=255)
            aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaa                                                                                         <  1:905992/100‑5 (MQ=39)
             tttCCGCCACAGTTGCGTCAATTCAGCCGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaac                                                                                        <  2:2510353/100‑6 (MQ=38)
             tttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaac                                                                                        <  1:3519290/100‑6 (MQ=39)
              ttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTAATGAGCAACCATGCTTGCTCAtcaaacg                                                                                       <  1:1535392/100‑7 (MQ=38)
                ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTTATGAGCAACCATGCTTGCTCAtcaaa                                                                                         <  2:1354287/96‑5 (MQ=39)
                ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                                     >  1:2408617/1‑100 (MQ=255)
                ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                                     >  1:1933792/1‑100 (MQ=255)
                 cGCCACAGTTGCGCGAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                    <  2:1567115/100‑1 (MQ=255)
                 cGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                    >  1:2329575/1‑100 (MQ=255)
                 cGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                    >  1:523307/1‑100 (MQ=255)
                  gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATNGTTGCTCATTGCCATTCTg                                                                                   <  1:2024998/100‑1 (MQ=255)
                  gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                                     >  1:270055/1‑98 (MQ=255)
                  gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTg                                                                                   >  2:2366023/1‑100 (MQ=255)
                    cacaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaa                                                                                 >  1:2822639/1‑100 (MQ=255)
                        gTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCt                                                                             >  2:3369455/1‑100 (MQ=255)
                         ttGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTg                                                                            <  2:564652/100‑1 (MQ=255)
                           gcgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                          >  1:3086024/1‑100 (MQ=255)
                           gcgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                          >  1:1648933/1‑100 (MQ=255)
                           gcgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGa                                                                           <  2:2054836/99‑1 (MQ=255)
                             gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATAAACGAGGATGGTCTCTACTTGCCAATCGGACACCTTAAc                                                                         >  2:243159/1‑99 (MQ=255)
                              ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCt                                                                       >  2:1163102/1‑100 (MQ=255)
                                aaTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                     <  2:1105841/100‑1 (MQ=255)
                                aaTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                     >  1:1805815/1‑100 (MQ=255)
                                      gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                               <  2:2205240/100‑1 (MQ=255)
                                         gAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                            <  1:3181892/100‑1 (MQ=255)
                                          aaCGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCt                                                           <  1:1271189/100‑1 (MQ=255)
                                            cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                         <  2:3484232/100‑1 (MQ=255)
                                             gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCga                                                        >  2:2199152/1‑100 (MQ=255)
                                              aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                       <  1:2288624/100‑1 (MQ=255)
                                              aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTCGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                       >  2:1491931/1‑100 (MQ=255)
                                              aaTGACGAAGCGCATTTGGTATATGGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                       <  2:602591/100‑1 (MQ=255)
                                               aTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgaga                                                      <  1:3487729/100‑1 (MQ=255)
                                                 gACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATa                                                    <  2:1425421/100‑1 (MQ=255)
                                                 gACGAAGCGCATTTGCTATATTGTCAGGGGACGATTGATGAGCAAGCATGCTTGCTCATTGACATTCTGAAACCTGAACCTCATAAATTGACTCGAGAt                                                     >  2:307583/1‑99 (MQ=255)
                                                  aCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgagaaga                                                   >  2:2024183/1‑97 (MQ=255)
                                                   cGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaac                                                  <  1:768707/100‑1 (MQ=255)
                                                    gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaaca                                                 <  1:3163259/100‑1 (MQ=255)
                                                    gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaaca                                                 <  1:851934/100‑1 (MQ=255)
                                                      aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCCTTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaac                                                  <  2:2252237/97‑1 (MQ=255)
                                                      aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaacaac                                               >  1:1865307/1‑100 (MQ=255)
                                                      aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGACTCGAGATaacaac                                               <  1:113229/100‑1 (MQ=255)
                                                          cATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaaa                                           <  2:2240569/100‑1 (MQ=255)
                                                          cATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaaa                                           <  2:2239871/100‑1 (MQ=255)
                                                            tttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATg                                         <  2:3348607/100‑1 (MQ=255)
                                                             ttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                        <  1:793390/100‑1 (MQ=255)
                                                             ttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                        <  2:1648925/100‑1 (MQ=255)
                                                             ttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                        >  2:3618188/1‑100 (MQ=255)
                                                                ttataTTGTCAGGTGTCGTTGTATGGGCATGCATGGTTTCTCATTGCCAGGCTGTAAGTTGAACCTCATAAACGGGCTCGAGATAACAACCAAAt                                          >  2:3039690/2‑95 (MQ=255)
                                                                gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                                     >  2:2806248/1‑100 (MQ=255)
                                                                    atTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATCCATAAAAt                                 >  2:3469482/1‑100 (MQ=255)
                                                                       gTCAGGGGGTGTTTGATGAGCAAGCAGGTGTGCTCATTGCCATTCTCAAACCTGAACATCATAAACTGTCTCGAGATAACAACCAAATGCATAAAATTgg                              <  2:376113/100‑1 (MQ=255)
                                                                       gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                              >  2:3584485/1‑100 (MQ=255)
                                                                        tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                             <  1:2206338/100‑1 (MQ=255)
                                                                        tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                             <  1:12388/100‑1 (MQ=255)
                                                                           gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                             <  1:2158048/97‑1 (MQ=255)
                                                                           gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaa                          >  1:588290/1‑100 (MQ=255)
                                                                             gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGCGATAACAACCAAATGCATAAAATTGGGTAAAt                        >  1:3438750/1‑100 (MQ=255)
                                                                               aatggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCGCGAGTTAACAACCAAATGCATAAAATTGGGAAAATgg                      >  2:2697038/8‑100 (MQ=38)
                                                                                atggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                     >  1:754508/7‑100 (MQ=38)
                                                                                 gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa                    >  2:410332/1‑100 (MQ=255)
                                                                                  tttGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAg                   <  2:73301/100‑1 (MQ=255)
                                                                                  ggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAg                   >  1:3472581/5‑100 (MQ=39)
                                                                                    tgatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTAGGAAAATGGCAGaa                 >  2:2094392/1‑100 (MQ=255)
                                                                                    tgatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAAACAAATGCATAAAATTGGGAAAATGGCAGaa                 >  2:2683882/1‑100 (MQ=255)
                                                                                    caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa                 >  1:2770981/3‑100 (MQ=255)
                                                                                     gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg                >  1:494517/1‑100 (MQ=255)
                                                                                     aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg                >  1:782987/2‑100 (MQ=255)
                                                                                     aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg                >  2:214427/2‑100 (MQ=255)
                                                                                     aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg                >  2:2880937/2‑100 (MQ=255)
                                                                                      atgaGCAACCATGCTTGCTCATTGGCATGCTGAAACCTGAACCTCATAAACTGGCGCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc               >  2:1209812/1‑100 (MQ=255)
                                                                                       tgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCg              >  2:880221/1‑100 (MQ=255)
                                                                                        gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGt             >  2:1219943/1‑100 (MQ=255)
                                                                                          gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt           >  2:747147/1‑100 (MQ=255)
                                                                                           caagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTc          <  1:241042/100‑1 (MQ=255)
                                                                                             agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc        >  2:1144416/1‑100 (MQ=255)
                                                                                             agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc        >  1:970879/1‑100 (MQ=255)
                                                                                             agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCACATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGTAGAAGCGTTTCGc        >  1:2565700/1‑100 (MQ=255)
                                                                                              gcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa       <  2:2745063/100‑1 (MQ=255)
                                                                                               caTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCAt      <  1:2016138/100‑1 (MQ=255)
                                                                                                aTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATg     <  2:3330828/100‑1 (MQ=255)
                                                                                                  ggTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCg   >  2:614135/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGT  >  NC_000913/252495‑252693

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