BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,519,341	A→G	noncoding (69/76 nt)	lysV →	tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,519,341	0	A	G	99.1%	268.5 / ‑4.5	111	noncoding (69/76 nt)	lysV	tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (45/65); total (46/65)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

Only 100 of 114 total aligned reads displayed.

CCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGC  >  NC_000913/2519252‑2519429

cccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGccc                                                                                       <  1:4405853/93‑1 (MQ=37)

 ccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCcca                                                                                      >  1:155142/1‑93 (MQ=37)

 ccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCcca                                                                                      >  1:3888210/1‑93 (MQ=37)

 ccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCcca                                                                                      >  2:31893/1‑93 (MQ=37)

     catcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCaa                                                                                  <  1:2635196/93‑1 (MQ=255)

      ctcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAAt                                                                                 <  2:3560916/92‑1 (MQ=255)

      atcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAAt                                                                                 <  2:1842765/93‑1 (MQ=255)

      atcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAAt                                                                                 <  1:431316/93‑1 (MQ=255)

      atcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAAt                                                                                 >  1:3926578/1‑93 (MQ=255)

       tcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATg                                                                                <  2:2766721/93‑1 (MQ=255)

             accaccTCGGGTCGTTAGCTCAGTTGTTAGAGCAGTTAACTTTTAATCATGTGTGCGCAGGTTCGAATCATGCACGGCCCACCAATGTaaaaa                                                                          <  2:3996417/93‑1 (MQ=11)

             accacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaaaa                                                                          >  2:3525739/1‑93 (MQ=255)

               cacTTCGGGTCGTTAGCTCAGTTGGTAGAGCATTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAg                                                                        <  1:2340697/93‑1 (MQ=255)

               cacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAg                                                                        <  1:3212485/93‑1 (MQ=255)

               cacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAg                                                                        <  1:1437415/93‑1 (MQ=255)

                acTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgc                                                                       >  1:885140/1‑93 (MQ=25)

                acTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgc                                                                       <  1:3069801/93‑1 (MQ=25)

                  ttCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgcgc                                                                     >  2:1116079/1‑93 (MQ=18)

                   tCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGTTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGcc                                                                    <  1:2934856/93‑1 (MQ=12)

                    cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGccc                                                                   <  2:3588441/93‑1 (MQ=14)

                    cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGccc                                                                   <  1:3791352/93‑1 (MQ=14)

                                    ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa                                                   <  1:1071029/93‑1 (MQ=14)

                                    ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa                                                   >  1:2889100/1‑93 (MQ=14)

                                    ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa                                                   >  2:1895623/1‑93 (MQ=14)

                                     tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTAc                                                  >  1:721168/1‑93 (MQ=17)

                                     tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTAc                                                  <  1:721107/93‑1 (MQ=17)

                                     tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTAc                                                  <  2:1624573/93‑1 (MQ=17)

                                       gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTa                                                >  2:2515932/1‑93 (MQ=18)

                                        tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTAt                                               >  1:3292919/1‑93 (MQ=18)

                                         agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATc                                              <  2:2810870/93‑1 (MQ=18)

                                         agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATc                                              <  1:4039514/93‑1 (MQ=18)

                                          gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCt                                             >  1:4158417/1‑93 (MQ=18)

                                          gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCt                                             <  1:2536539/93‑1 (MQ=21)

                                          gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCt                                             >  1:3990083/1‑93 (MQ=18)

                                           agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTg                                            <  2:4853/93‑1 (MQ=21)

                                           agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTg                                            >  1:4285454/1‑93 (MQ=18)

                                            gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGc                                           <  1:2841415/93‑1 (MQ=18)

                                            gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGc                                           >  1:3448427/1‑93 (MQ=18)

                                            gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGc                                           >  2:3860783/1‑93 (MQ=18)

                                             cAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCg                                          <  2:4593132/93‑1 (MQ=21)

                                              aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                         <  1:245232/93‑1 (MQ=21)

                                              aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                         <  1:2475037/93‑1 (MQ=18)

                                              aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                         >  1:2691866/1‑93 (MQ=18)

                                              aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                         >  2:2765496/1‑93 (MQ=33)

                                               gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGAt                                        >  1:2409414/1‑93 (MQ=18)

                                               gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGAt                                        <  1:4579450/93‑1 (MQ=18)

                                                ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATa                                       <  2:2661382/93‑1 (MQ=18)

                                                 tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAc                                      >  2:1779913/1‑93 (MQ=18)

                                                 tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAc                                      >  2:3064318/1‑93 (MQ=18)

                                                  taCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACt                                     <  1:4190909/92‑1 (MQ=18)

                                                   aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTc                                    <  1:3854606/93‑1 (MQ=21)

                                                   aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTc                                    <  2:155192/93‑1 (MQ=18)

                                                    cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCa                                   <  2:603215/93‑1 (MQ=21)

                                                    cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCa                                   >  1:834650/1‑93 (MQ=18)

                                                    cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCa                                   >  1:1451584/1‑93 (MQ=18)

                                                     ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaa                                  <  2:3596245/93‑1 (MQ=18)

                                                      tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                 >  1:691725/1‑93 (MQ=18)

                                                        tAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGa                               >  1:1217079/1‑93 (MQ=18)

                                                        tAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGa                               <  1:322023/93‑1 (MQ=21)

                                                        tAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGa                               <  1:1250262/93‑1 (MQ=21)

                                                         aaTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAt                              <  2:3811931/93‑1 (MQ=18)

                                                         aaTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAt                              <  2:3365139/93‑1 (MQ=18)

                                                          aTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                             <  1:1559442/93‑1 (MQ=18)

                                                          aTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                             >  1:2467711/1‑93 (MQ=18)

                                                          aTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                             <  1:302059/93‑1 (MQ=21)

                                                           tCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTc                            <  2:589052/93‑1 (MQ=21)

                                                           tCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTc                            <  1:2150238/93‑1 (MQ=21)

                                                           tCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTc                            <  1:3851631/93‑1 (MQ=21)

                                                           tCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTc                            >  2:2482785/1‑93 (MQ=18)

                                                            cAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCg                           <  2:773702/93‑1 (MQ=21)

                                                               ttGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAc                        <  1:2591798/93‑1 (MQ=18)

                                                                tGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAcc                       <  1:2675027/93‑1 (MQ=18)

                                                                 ggTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCt                      >  1:2765750/1‑93 (MQ=18)

                                                                   tCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgc                    >  1:4502606/1‑93 (MQ=18)

                                                                   tCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgc                    <  1:2762067/93‑1 (MQ=18)

                                                                       aGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcagca                <  1:1079857/93‑1 (MQ=18)

                                                                       aGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcagca                <  1:4666572/93‑1 (MQ=21)

                                                                        ggTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGcagcag               <  2:1820622/93‑1 (MQ=21)

                                                                        ggTTCGAATCCTGCACGGACCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGcagcag               >  2:2766760/1‑93 (MQ=16)

                                                                          ttCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGt             >  2:2600386/1‑93 (MQ=18)

                                                                           tCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt            <  2:1267651/93‑1 (MQ=18)

                                                                            cGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt           >  2:4541626/1‑93 (MQ=18)

                                                                             gAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttg          >  1:4316667/1‑93 (MQ=18)

                                                                              aaTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAGGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga         <  2:423175/93‑1 (MQ=16)

                                                                              aaTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga         <  2:2528378/93‑1 (MQ=18)

                                                                              aaTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga         <  2:68392/93‑1 (MQ=18)

                                                                              aaTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga         >  1:3334476/1‑93 (MQ=18)

                                                                              aaTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga         <  1:3473647/93‑1 (MQ=18)

                                                                              aaTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga         >  2:2668068/1‑93 (MQ=18)

                                                                               aTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag        <  1:2962891/93‑1 (MQ=18)

                                                                               aTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag        >  2:4119638/1‑93 (MQ=18)

                                                                               aTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag        <  2:1332447/93‑1 (MQ=21)

                                                                               aTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag        <  2:697035/93‑1 (MQ=21)

                                                                                tCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt       >  2:1441885/1‑93 (MQ=18)

                                                                                tCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt       >  2:1501866/1‑93 (MQ=18)

                                                                                tCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgaga       >  1:1917074/1‑92 (MQ=18)

                                                                                  cTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttg     >  1:4631917/1‑93 (MQ=18)

                                                                                  cTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttg     >  1:1989264/1‑93 (MQ=18)

                                                                                     cACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc  <  1:1939431/93‑1 (MQ=21)

                                                                                     cACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc  <  2:602973/93‑1 (MQ=18)

CCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGC  >  NC_000913/2519252‑2519429

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