Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,250,107 T→A 5.8% E204V (GAG→GTG)  nth ← DNA glycosylase/apyrimidinic (AP) lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,250,1070TA 5.8% 241.2 / 2.8 79E204V (GAG→GTG) nthDNA glycosylase/apyrimidinic (AP) lyase
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base A (2/3);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTATCGGTACTCATGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCTTTTCCACTTCCAGCACCGT  >  NC_002947/1249959‑1250255
                                                                                                                                                    |                                                                                                                                                    
cTTTATCGGTACTCATGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACt                                                                                                                                                      >  7:260184/1‑149 (MQ=255)
     tCGGTCCTCATGAGTCGCTCCAAGCCACGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGca                                                                                                                                                  <  8:227662/148‑1 (MQ=255)
      cGGTACTCATGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGcac                                                                                                                                                 >  8:148058/1‑148 (MQ=255)
        gTACTCATGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACaaa                                                                                                                                              >  3:214405/1‑149 (MQ=255)
           cTCATGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCGAGTACAAGCACaaaac                                                                                                                                            <  8:224822/148‑3 (MQ=255)
            tCATGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCtt                                                                                                                                          >  6:185517/1‑149 (MQ=255)
            tCATGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCGAGTACAAGCACaaaacct                                                                                                                                          <  8:53985/149‑5 (MQ=255)
            gcATGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATATCCAGTAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCGAGTACAAGCACaaaagct                                                                                                                                          <  8:359268/148‑5 (MQ=255)
              aTGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTc                                                                                                                                         >  3:334654/1‑148 (MQ=255)
              aTGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTc                                                                                                                                         >  8:325461/1‑148 (MQ=255)
              aTGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTc                                                                                                                                         <  7:325461/148‑1 (MQ=255)
                 aGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATc                                                                                                                                     <  7:123602/149‑1 (MQ=255)
                          aagccaagGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTg                                                                                                                            >  3:10500/1‑149 (MQ=255)
                          aagccaagGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTg                                                                                                                            >  6:132025/1‑149 (MQ=255)
                           agccaagGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGc                                                                                                                           <  1:199191/149‑1 (MQ=255)
                            gccaagGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTg                                                                                                                                  <  4:273590/141‑1 (MQ=255)
                            gccaagGGGCCATTGGCCGTCCATATAACCAGAAAGCAGAAGCTCGTAAAAAAGATTTTTTCAATCGATGAGAGTGGTTTTTGGCTATCGCCTCAATCGTCGCAGGTTTTGTGCTTGTACTCGCACAAATCTTAGATCCTg                                                                                                                                  >  3:273590/1‑141 (MQ=255)
                             ccaagGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCg                                                                                                                         >  1:225748/1‑149 (MQ=255)
                             ccaagGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAAGCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCg                                                                                                                         >  5:298002/1‑149 (MQ=255)
                              caagGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGc                                                                                                                        <  4:214405/149‑1 (MQ=255)
                                 ggggCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCAcc                                                                                                                     >  8:316849/1‑149 (MQ=255)
                                    gCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTg                                                                                                                  <  7:148058/149‑1 (MQ=255)
                                       aTTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCt                                                                                                               >  3:13675/1‑149 (MQ=255)
                                           gCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTCGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTgcg                                                                                                           >  2:257666/1‑149 (MQ=255)
                                                    ataACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTgg                                                                                                    <  2:225748/147‑1 (MQ=255)
                                                      aaCCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGc                                                                                                                                                   >  7:360479/1‑98 (MQ=255)
                                                      aaCCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGc                                                                                                                                                   <  8:360479/98‑1 (MQ=255)
                                                      aaCCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGcaca                                                                                                >  7:270636/1‑149 (MQ=255)
                                                       aCCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCacac                                                                                               >  7:7686/1‑149 (MQ=255)
                                                        ccAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACg                                                                                              >  4:313857/1‑149 (MQ=255)
                                                        ccAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACg                                                                                              >  4:319250/1‑149 (MQ=255)
                                                                cAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGccc                                                                                      <  5:185517/149‑1 (MQ=255)
                                                                     gCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCCTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCa                                                                                 <  6:279025/149‑1 (MQ=255)
                                                                      cTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAg                                                                                <  5:132025/149‑1 (MQ=255)
                                                                             aaaaGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGc                                                                         <  4:10500/149‑1 (MQ=255)
                                                                                aGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAg                                                                      >  6:135294/1‑149 (MQ=255)
                                                                                 gATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTATCGCCCGTGCAGGATCAGCCAgt                                                                     <  8:270636/149‑1 (MQ=255)
                                                                                  aTTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGAACCTTGGCTTGCGCGCCTGGCACACGTATCGACCGTGCATGATAAGCCAgtg                                                                    >  3:5334/1‑149 (MQ=255)
                                                                                   ttttttCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAgtg                                                                    <  4:13675/148‑1 (MQ=255)
                                                                                          aaTCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGt                                                            >  7:213208/1‑149 (MQ=255)
                                                                                           ctCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCa                                                                                 <  2:220335/126‑1 (MQ=255)
                                                                                           aTCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAg                                                                                <  6:100355/128‑1 (MQ=255)
                                                                                           aTCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAg                                                                                >  5:100355/1‑128 (MQ=255)
                                                                                            tCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCg                                                          <  8:260184/149‑1 (MQ=255)
                                                                                               aTGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAg                                                                                >  8:319899/1‑124 (MQ=255)
                                                                                               aTGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAg                                                                                <  7:319899/124‑1 (MQ=255)
                                                                                                tttcaGTTGTTTTTTGCTATCCGCTCAATCGTCGCAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGc                                                       <  4:5334/144‑1 (MQ=255)
                                                                                                  tgAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAg                                                     <  3:313857/148‑1 (MQ=255)
                                                                                                  tgAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAg                                                     <  3:319250/148‑1 (MQ=255)
                                                                                                   gAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGATGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAgg                                                    >  3:227583/1‑148 (MQ=255)
                                                                                                    aGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTa                                                  <  1:257666/149‑1 (MQ=255)
                                                                                                    aGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACACGCACAAATCTTCGATCCTGCCGCTGCCGAACCTTGGCTTGCGAGCCTGGCACACCTAGCGCGCGTGCAGGATCATCCACTGGGGCGCGTCGAGctgc                                                    >  8:255296/1‑144 (MQ=255)
                                                                                                     gTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAg                                                 >  3:137258/1‑149 (MQ=255)
                                                                                                      tGGTTTTTGGCTATCGGCTCAATCATCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGt                                                >  3:307464/1‑149 (MQ=255)
                                                                                                             ttgCTCTCGGCTAAATCGTCGGAGGTGTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGTCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATAAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTgg                                          <  2:80630/146‑1 (MQ=255)
                                                                                                             tGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTgg                                          <  6:298002/148‑1 (MQ=255)
                                                                                                             tGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTgg                                          <  4:334654/148‑1 (MQ=255)
                                                                                                             tGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTgg                                          <  7:316849/148‑1 (MQ=255)
                                                                                                              ggCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTggg                                         >  3:93073/1‑148 (MQ=255)
                                                                                                                        tCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGca                                                                                                  <  8:87147/81‑1 (MQ=255)
                                                                                                                        tCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGca                                                                                                  >  7:87147/1‑81 (MQ=255)
                                                                                                                          aaTCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGa                              <  5:135294/147‑1 (MQ=255)
                                                                                                                          aaTCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAg                             >  7:3302/1‑148 (MQ=255)
                                                                                                                           aTCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTgg                                          >  7:358622/1‑134 (MQ=255)
                                                                                                                           aTCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTgg                                          >  7:33908/1‑134 (MQ=255)
                                                                                                                           aTCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTgg                                          <  8:358622/134‑1 (MQ=255)
                                                                                                                           aTCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTgg                                          <  8:33908/134‑1 (MQ=255)
                                                                                                                               tcgGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGttt                       <  8:7686/149‑1 (MQ=255)
                                                                                                                               tcgGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGttt                       <  8:213208/149‑1 (MQ=255)
                                                                                                                                 ggAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTc                      <  8:3302/148‑1 (MQ=255)
                                                                                                                                  gaagatttGTGCGAGTACAAGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCtt                    >  8:250928/6‑149 (MQ=255)
                                                                                                                                   aGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCttt                   >  3:93339/1‑149 (MQ=255)
                                                                                                                                    ggTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCtttt                  >  6:336157/1‑149 (MQ=255)
                                                                                                                                     gtgttGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTTGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCTTTTc                 <  7:255296/146‑1 (MQ=255)
                                                                                                                                     gTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCtttt                  <  4:227583/148‑1 (MQ=255)
                                                                                                                                       tttGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACAAGTAGCGCACGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACATGAGCAGTTTCTTTTcga               >  3:97376/1‑147 (MQ=255)
                                                                                                                                           tgCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCTTTTCCACTTc           <  8:39836/149‑1 (MQ=255)
                                                                                                                                             cTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCTTTTCCACTTCCa         >  3:57008/1‑149 (MQ=255)
                                                                                                                                             cTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCTTTTCCACTTCCa         <  4:93073/149‑1 (MQ=255)
                                                                                                                                                   cTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCTTTTCCACTTCCAGCACCg   >  1:59130/1‑149 (MQ=255)
                                                                                                                                                    tCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCTTTTCCACTTCCAGCACCGt  >  2:268951/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTTATCGGTACTCATGAGTCGCTCCAAGCCAAGGGGCCATTGGCCGTCCTTATAACCAGAAAGCAGACGCTGGTAAAAAAGATTTTTTCAATCGATGTGAGTGGTTTTTGGCTATCGGCTCAATCGTCGGAGGTTTTGTGCTTGTACTCGCACAAATCTTCGATCCTGCAGCTGCCGCACCTTGGCTTGCGCGCCTGGCACACGTAGCGCCCGTGCAGGATCAGCCAGTGGTGGGCGTCGAGCAGGTAGTCCTTGGGTACAAACTTGAGCAGTTTCTTTTCCACTTCCAGCACCGT  >  NC_002947/1249959‑1250255

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑