BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,959,281	G→A	E330E (GAG→GAA)	rep →	rep helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,959,281	0	G	A	100.0%	1008.8 / NA	264	E330E (GAG→GAA)	rep	rep helicase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (125/139); total (125/139)

Only 100 of 270 total aligned reads displayed.

GCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAA  >  NC_000913/3959191‑3959370

  gTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGa                                                                                       >  1:3483114/1‑93 (MQ=255)

  gTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGa                                                                                       <  2:4944818/93‑1 (MQ=255)

  gTCTGTTCTCCGAACT‑GGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAt                                                                                      >  1:5621969/1‑93 (MQ=255)

    cTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGATCGCGTTACTGGCGAACTGATc                                                                                     <  1:1382250/93‑1 (MQ=255)

     tGTTCTCCGAACTGGGTTATGGAGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCg                                                                                    <  2:2365329/93‑1 (MQ=255)

      gTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGc                                                                                   >  1:2301003/1‑93 (MQ=255)

       ttCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGcc                                                                                  <  2:1546960/93‑1 (MQ=255)

        tctcCGAACTCGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGccc                                                                                 >  1:2556689/1‑93 (MQ=255)

         ctcCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCa                                                                                >  2:3377308/1‑93 (MQ=255)

           ccGTACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATc                                                                              <  2:3836521/93‑1 (MQ=255)

           ccGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGGGCGTTACTGGCGAACTGATCGCCCATc                                                                              >  1:1699661/1‑93 (MQ=255)

           ccGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATc                                                                              <  2:362819/93‑1 (MQ=255)

              aaCTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACt                                                                           >  1:732272/1‑93 (MQ=255)

                cTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGATCGCGTTACTGGCGAACTGATCGCCCATCACTTc                                                                         <  1:976966/93‑1 (MQ=255)

                cTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTc                                                                         <  1:1113879/93‑1 (MQ=255)

                cTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTc                                                                         <  2:525440/93‑1 (MQ=255)

                  gggTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGt                                                                       <  1:5222293/93‑1 (MQ=255)

                    gTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCa                                                                     >  1:338996/1‑93 (MQ=255)

                     ttATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCaa                                                                    <  2:623591/93‑1 (MQ=255)

                       aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATa                                                                  <  2:2136308/93‑1 (MQ=255)

                       aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATa                                                                  >  2:3244256/1‑93 (MQ=255)

                         ggCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATaaa                                                                >  2:1256435/1‑93 (MQ=255)

                           cgcgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAAc                                                              <  1:5003244/93‑1 (MQ=255)

                             cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGc                                                            <  1:5827308/93‑1 (MQ=255)

                             cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGc                                                            <  2:1905321/93‑1 (MQ=255)

                             cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGc                                                            >  2:6502323/1‑93 (MQ=255)

                              ggAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCa                                                           <  2:6215583/93‑1 (MQ=255)

                              ggAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCa                                                           <  2:2275838/93‑1 (MQ=255)

                              ggAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCa                                                           >  1:2819533/1‑93 (MQ=255)

                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          <  2:1817768/93‑1 (MQ=255)

                                aGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGt                                                         >  1:6019600/1‑93 (MQ=255)

                                aGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGt                                                         <  2:16877/93‑1 (MQ=255)

                                aGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGt                                                         >  2:943026/1‑93 (MQ=255)

                                 gCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTa                                                        <  2:2565421/93‑1 (MQ=255)

                                   tAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACa                                                      >  2:6331655/1‑93 (MQ=255)

                                   tAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACa                                                      >  2:6018880/1‑93 (MQ=255)

                                   tAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACa                                                      >  1:1293273/1‑93 (MQ=255)

                                    aaaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACaa                                                     <  1:4357060/93‑1 (MQ=255)

                                    aaaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACaa                                                     >  2:4299110/1‑93 (MQ=255)

                                      aaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAg                                                   <  1:3191294/93‑1 (MQ=255)

                                      aaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAg                                                   >  2:737684/1‑93 (MQ=255)

                                       aGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGa                                                  <  1:6354241/93‑1 (MQ=255)

                                        gTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAt                                                 >  1:1201219/1‑93 (MQ=255)

                                        gTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAt                                                 <  2:4178647/93‑1 (MQ=255)

                                        gTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAt                                                 >  2:5650834/1‑93 (MQ=255)

                                        gTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAt                                                 <  1:2069068/93‑1 (MQ=255)

                                         tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAtt                                                >  1:6283492/1‑93 (MQ=255)

                                         tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAtt                                                <  1:863495/93‑1 (MQ=255)

                                          aTTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTa                                               >  2:4489131/1‑93 (MQ=255)

                                           ttAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTAc                                              <  1:1569520/93‑1 (MQ=255)

                                            tAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGTCCATCACTTAGTCAATAAAACGCAGTACAAAGATTACg                                             >  2:2864056/1‑93 (MQ=255)

                                            tAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACg                                             >  1:6087997/1‑93 (MQ=255)

                                               gcgcGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCa                                          >  1:2992064/1‑93 (MQ=255)

                                               gcgcGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCa                                          >  1:3741561/1‑93 (MQ=255)

                                                 gcgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCAtt                                        <  1:6222907/93‑1 (MQ=255)

                                                 gcgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCAtt                                        <  2:2201535/93‑1 (MQ=255)

                                                 ccgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCAtt                                        <  1:5799978/92‑1 (MQ=255)

                                                  cgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTc                                       <  1:1422038/93‑1 (MQ=255)

                                                  cgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTc                                       <  2:3535362/93‑1 (MQ=255)

                                                  cgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTc                                       <  2:3181692/93‑1 (MQ=255)

                                                   gAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCt                                      >  1:2694092/1‑93 (MQ=255)

                                                   gAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCt                                      <  2:1600470/93‑1 (MQ=255)

                                                    aaTAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCtt                                     <  1:2516127/93‑1 (MQ=255)

                                                    aaTAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCtt                                     <  2:2998844/93‑1 (MQ=255)

                                                    aaTAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCtt                                     >  2:3293338/1‑93 (MQ=255)

                                                      tAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTa                                   <  2:5415447/93‑1 (MQ=255)

                                                      tAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTa                                   >  1:2949637/1‑93 (MQ=255)

                                                       aaCGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTAt                                  >  2:6486905/1‑93 (MQ=255)

                                                       aaCGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTAt                                  >  2:1051032/1‑93 (MQ=255)

                                                         cGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcg                                <  2:5134096/93‑1 (MQ=255)

                                                          gaagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgc                               <  1:6546901/93‑1 (MQ=255)

                                                          gaagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgc                               <  2:4086592/93‑1 (MQ=255)

                                                           aagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                              >  2:4000156/1‑93 (MQ=255)

                                                             gaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGt                            >  1:1667608/1‑93 (MQ=255)

                                                                cATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAc                         <  1:5715935/93‑1 (MQ=255)

                                                                 aTGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAcc                        >  2:3064637/1‑93 (MQ=255)

                                                                  tGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCa                       >  1:4147884/1‑93 (MQ=255)

                                                                    aGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATc                     >  1:4027489/1‑93 (MQ=255)

                                                                    aGGCTGAGCGCGTTACTGGCGAACTGATAGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATc                     >  1:2246034/1‑93 (MQ=255)

                                                                     ggCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCa                    >  1:4814400/1‑93 (MQ=255)

                                                                     ggCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCa                    <  1:1985338/93‑1 (MQ=255)

                                                                       cTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGt                  <  1:503497/93‑1 (MQ=255)

                                                                       cTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGt                  <  1:5306815/93‑1 (MQ=255)

                                                                        tGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTc                 >  1:5169028/1‑93 (MQ=255)

                                                                         gAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcg                >  1:3681876/1‑93 (MQ=255)

                                                                          aGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcgc               >  1:3908920/1‑93 (MQ=255)

                                                                            cgcgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCgg             <  2:4767039/93‑1 (MQ=255)

                                                                            cgcgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCgg             <  2:1298003/93‑1 (MQ=255)

                                                                            cgcgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCgg             <  1:2043030/93‑1 (MQ=255)

                                                                             gcgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCggg            <  2:4353990/93‑1 (MQ=255)

                                                                             gcgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCggg            <  1:3593010/93‑1 (MQ=255)

                                                                              cgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgt           >  2:4860640/1‑93 (MQ=255)

                                                                              cgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgt           >  1:3203424/1‑93 (MQ=255)

                                                                               gTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATGAGTCGCGGatg          >  1:5492106/1‑93 (MQ=255)

                                                                               gTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgtg          >  2:5816164/1‑93 (MQ=255)

                                                                               gTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgtg          >  1:2023931/1‑93 (MQ=255)

                                                                                 tACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgtggt        >  2:3898716/1‑91 (MQ=255)

                                                                                  aCTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGttt       <  2:4781180/93‑1 (MQ=255)

                                                                                    tGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGa     >  2:2293504/1‑93 (MQ=255)

                                                                                    tGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGa     >  1:1913052/1‑93 (MQ=255)

GCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAA  >  NC_000913/3959191‑3959370

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`