Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,656,386 C→T 8.7% intergenic (‑28/+150) pepB ← / ← iscX aminopeptidase B/Fe(2+) donor and activity modulator for cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,656,3860CT 8.7% 211.4 / 9.1 69intergenic (‑28/+150)pepB/iscXaminopeptidase B/Fe(2+) donor and activity modulator for cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base C (36/27);  new base T (2/4);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

CTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGC  >  NC_000913/2656244‑2656530
                                                                                                                                              |                                                                                                                                                
cTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAAAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCTGCCaac                                                                                                                                            <  3:426381/149‑3 (MQ=255)
    cccGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGc                                                                                                                                        >  4:564223/1‑149 (MQ=255)
     ccGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGtt                                                                                                                                           <  4:314563/145‑1 (MQ=255)
     ccGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGtt                                                                                                                                           >  3:314563/1‑145 (MQ=255)
       gTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGc                                                                                                                                     >  1:119141/1‑149 (MQ=255)
         tCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCTGCCAACTGGCGGccc                                                                                                                                   <  3:440068/149‑1 (MQ=255)
         tCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCTGCCAACTGGCGGccc                                                                                                                                   <  1:451122/149‑1 (MQ=255)
               gCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAg                                                                                                                                >  2:382019/1‑146 (MQ=255)
               gCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAg                                                                                                                                <  1:382019/146‑1 (MQ=255)
                cAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGAtt                                                                                                                            <  3:140904/149‑1 (MQ=255)
                 aGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTTATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGAttt                                                                                                                           <  4:403903/149‑1 (MQ=255)
                           cGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTaaa                                                                                                                 <  1:128872/149‑1 (MQ=255)
                            gTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAAt                                                                                                                >  2:82055/1‑149 (MQ=255)
                            gTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAAt                                                                                                                >  2:591650/1‑149 (MQ=255)
                               antattAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAt                                                                                                             >  3:113192/1‑149 (MQ=255)
                                ttattaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAtt                                                                                                            <  1:136009/149‑1 (MQ=255)
                                      aTGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGtt                                                                                                      <  2:71013/149‑1 (MQ=255)
                                         cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                   <  1:229490/149‑1 (MQ=255)
                                         cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGAGTTGCa                                                                                                   <  2:422647/149‑1 (MQ=255)
                                               gTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAgg                                                                                             <  2:299373/149‑1 (MQ=255)
                                                ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                            >  4:6455/1‑149 (MQ=255)
                                                ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                            <  4:26154/149‑1 (MQ=255)
                                                ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                            >  2:104682/1‑149 (MQ=255)
                                                 tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcg                                                                                           <  1:553465/149‑1 (MQ=255)
                                                   cTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcggc                                                                                         >  3:339435/1‑149 (MQ=255)
                                                    tttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCa                                                                                        >  3:379543/1‑149 (MQ=255)
                                                    tttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCa                                                                                        >  4:425940/1‑149 (MQ=255)
                                                     ttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCaa                                                                                       <  3:126172/149‑1 (MQ=255)
                                                     ttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCaa                                                                                       <  4:32770/149‑1 (MQ=255)
                                                           ccccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgagt                                                                                 >  2:99188/1‑149 (MQ=255)
                                                            cccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgagtg                                                                                <  1:104682/149‑1 (MQ=255)
                                                               aGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAg                                                                             <  2:455919/149‑1 (MQ=255)
                                                                    gcgTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCa                                                                        <  4:293877/149‑1 (MQ=255)
                                                                     cgTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAg                                                                       <  1:275567/149‑1 (MQ=255)
                                                                       tCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGa                                                                     <  2:119141/149‑1 (MQ=255)
                                                                          gCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGGGCa                                                                  >  1:591083/1‑149 (MQ=255)
                                                                            aGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATa                                                                >  3:34599/1‑149 (MQ=255)
                                                                              gTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAt                                                                                                             <  3:258029/102‑1 (MQ=255)
                                                                              gTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAt                                                                                                             >  4:258029/1‑102 (MQ=255)
                                                                                 gggTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATaa                                                           <  1:266995/149‑1 (MQ=255)
                                                                                   gTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACt                                                         >  1:445952/1‑149 (MQ=255)
                                                                                           ggTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGccc                                                                                                                                   >  4:133476/1‑67 (MQ=255)
                                                                                           ggTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGccc                                                                                                                                   <  3:133476/67‑1 (MQ=255)
                                                                                           ggTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGAAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACt                                                 >  2:211513/1‑149 (MQ=255)
                                                                                            gTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTg                                                >  4:558628/1‑149 (MQ=255)
                                                                                             tAATCTTCATCGCTTCTGTCATCTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTgg                                               >  2:287591/1‑149 (MQ=255)
                                                                                              aaTCTTCATCGCTTCTGTCATTCTAGTTATCCTTCTTTGTAAACGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGAGACTggg                                              >  2:219223/1‑149 (MQ=255)
                                                                                              aaTCTTCATCGCTTCTGTCATCTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTggg                                              <  1:287591/149‑1 (MQ=255)
                                                                                                      tCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGa                                      >  2:57404/1‑149 (MQ=255)
                                                                                                      tCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGa                                      <  1:57404/149‑1 (MQ=255)
                                                                                                              gTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGcc                              >  2:309052/1‑149 (MQ=255)
                                                                                                              gTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACAATGTGACTGGGGTGAACGAGCGCAGcc                              <  3:391730/149‑1 (MQ=255)
                                                                                                                  ttttAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACg                          >  2:302023/1‑149 (MQ=255)
                                                                                                                  ttttAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACg                          >  1:503014/1‑149 (MQ=255)
                                                                                                                      aGTTAACCTGCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAAGTGACTGGGGTGAACGAGCGCAGCCAACGCATc                      <  1:99188/149‑1 (MQ=255)
                                                                                                                            ccTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTgcggc                >  3:313553/1‑149 (MQ=255)
                                                                                                                              ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt              >  1:169001/1‑149 (MQ=255)
                                                                                                                              ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt              >  1:215872/1‑149 (MQ=255)
                                                                                                                              ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt              >  2:571488/1‑149 (MQ=255)
                                                                                                                              ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt              >  3:509280/1‑149 (MQ=255)
                                                                                                                                 tttGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGaa           >  2:72/1‑149 (MQ=255)
                                                                                                                                  ttGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAg          >  2:349510/1‑149 (MQ=255)
                                                                                                                                    gTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCa        >  1:496994/1‑149 (MQ=255)
                                                                                                                                      aaaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGa                                      >  1:112236/1‑117 (MQ=255)
                                                                                                                                      aaaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGa                                      <  2:112236/117‑1 (MQ=255)
                                                                                                                                         ggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAcg   >  2:45203/2‑149 (MQ=255)
                                                                                                                                          gggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGaa           >  2:355809/1‑140 (MQ=255)
                                                                                                                                          gggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGaa           <  1:355809/140‑1 (MQ=37)
                                                                                                                                          gggCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAcgc  >  4:491486/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGC  >  NC_000913/2656244‑2656530

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: