BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,483,908	A→G	S942G (AGC→GGC)	hrpA →	helicase, ATP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,483,908	0	A	G	88.5%	437.9 / 42.2	183	S942G (AGC→GGC)	hrpA	helicase, ATP‑dependent
Reads supporting (aligned to +/- strand): ref base A (12/9); new base G (87/75); total (99/84)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 187 total aligned reads displayed.

AAACCGGTACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCTGAAAATT  >  NC_000913/1483818‑1483997

aaaCCGGTACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGc                                                                                         >  1:6474251/1‑93 (MQ=255)

aaaCCGGTACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGc                                                                                         >  1:6474250/1‑93 (MQ=255)

aaaCCGGTACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGc                                                                                         <  2:7622119/93‑1 (MQ=255)

aaaCCGGTACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCCTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGc                                                                                         >  2:5870950/1‑93 (MQ=255)

 aaCCGGTACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGcc                                                                                        <  1:6617406/93‑1 (MQ=255)

  aCCGGTACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCt                                                                                       <  2:2777233/93‑1 (MQ=255)

   ccGGTACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCtt                                                                                      <  2:3535779/93‑1 (MQ=255)

   ccGGTACGCCGTAATTTTGTACCCGCGCCAAAATATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCtt                                                                                      <  1:9957526/93‑1 (MQ=255)

     ggTACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGa                                                                                    >  2:8709115/1‑93 (MQ=255)

       tACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAgc                                                                                  <  1:598970/93‑1 (MQ=255)

       tACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAgc                                                                                  >  1:3823848/1‑93 (MQ=255)

       tACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAgc                                                                                  >  1:4988156/1‑93 (MQ=255)

        aCGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGCCGGCCTTGAgcg                                                                                 >  2:3531001/1‑93 (MQ=255)

           ccGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGa                                                                              <  2:7711643/93‑1 (MQ=255)

              tAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGtt                                                                           <  2:8028924/93‑1 (MQ=255)

                  tttGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTAcgg                                                                       <  2:2894356/93‑1 (MQ=255)

                   ttGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTAcggg                                                                      >  1:9276585/1‑92 (MQ=255)

                   ttGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTAcggc                                                                      <  2:1896222/93‑1 (MQ=255)

                   ttGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTAcggc                                                                      <  2:4196694/93‑1 (MQ=255)

                    tGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTAcggcg                                                                     <  2:7035743/93‑1 (MQ=255)

                      tACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGa                                                                   <  2:5814624/93‑1 (MQ=255)

                      tACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGa                                                                   >  1:1236136/1‑93 (MQ=255)

                      tACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGa                                                                   >  2:9593312/1‑93 (MQ=255)

                         ccGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGa                                                                <  1:2272527/93‑1 (MQ=255)

                          cgcgCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGAc                                                               >  1:6294639/1‑93 (MQ=255)

                          cgcgCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGAc                                                               <  1:7291523/93‑1 (MQ=255)

                          cgcgCCAAAATATGCCGAAGCGTTTTTAGGCCGCCTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGAc                                                               <  1:7767420/93‑1 (MQ=255)

                           gcgcCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGAcc                                                              <  2:50910/93‑1 (MQ=255)

                               cAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCg                                                          >  2:9468833/1‑93 (MQ=255)

                               cAAACTATGCCGAAGCGTTGTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCg                                                          <  2:4781603/93‑1 (MQ=255)

                                  aCTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTa                                                       >  2:5349924/1‑93 (MQ=255)

                                   cTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTAc                                                      <  1:91915/93‑1 (MQ=255)

                                   cTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTAc                                                      <  2:8452438/93‑1 (MQ=255)

                                    tATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTAcc                                                     >  2:1970286/1‑93 (MQ=255)

                                      tGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGt                                                   >  1:3797493/1‑93 (MQ=255)

                                       gCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGtt                                                  <  2:4148401/93‑1 (MQ=255)

                                        ccGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTg                                                 <  1:8698768/93‑1 (MQ=255)

                                        ccGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTg                                                 >  2:4358891/1‑93 (MQ=255)

                                         cGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGa                                                <  1:2172884/93‑1 (MQ=255)

                                         cGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGa                                                >  1:9091108/1‑93 (MQ=255)

                                          gAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGAc                                               <  1:7027234/93‑1 (MQ=255)

                                          gAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGAc                                               >  1:1843105/1‑93 (MQ=255)

                                          gAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGAc                                               >  1:8482468/1‑93 (MQ=255)

                                           aaGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGAcc                                              >  1:298724/1‑93 (MQ=255)

                                           aaGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGAcc                                              <  1:6546631/93‑1 (MQ=255)

                                            aGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACcg                                             >  1:5547371/1‑93 (MQ=255)

                                             gCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACcgc                                            >  1:3244984/1‑93 (MQ=255)

                                              cGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACcgcg                                           >  1:1023819/1‑93 (MQ=255)

                                              cGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGGTACCGTGGAcgggg                                           >  1:1023820/1‑89 (MQ=255)

                                               gTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGa                                          >  1:4273671/1‑93 (MQ=255)

                                                tttttAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGaa                                         <  2:7046751/93‑1 (MQ=255)

                                                tttttAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGaa                                         >  2:6412194/1‑93 (MQ=255)

                                                tttttAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGaa                                         >  1:8298335/1‑93 (MQ=255)

                                                 ttttAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGAACGGCGTTACCGTTGACCGCGAAg                                        >  1:7201200/1‑93 (MQ=255)

                                                  tttAGGCCGCGTCAAACCGCTGTAGTTACCGTTGCTCCGTGACCTTGAGCGCGAGTTACGGCGGATGACAGGCGGTACCGTTGGCCGCGaaat                                       >  2:2658351/1‑91 (MQ=255)

                                                  tttAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGa                                       <  2:1540748/93‑1 (MQ=255)

                                                  tttAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGa                                       <  1:9849991/93‑1 (MQ=255)

                                                  tttAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGa                                       <  2:1761344/93‑1 (MQ=255)

                                                   ttAGTCCCCGTAAAACCGCTGGAGTTACCGTTGCTCTACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTCCCGTTGACCGCGAAGAc                                      <  2:5126851/93‑1 (MQ=255)

                                                   ttAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGAc                                      <  2:9965087/93‑1 (MQ=255)

                                                   ttAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGAc                                      <  2:9121417/93‑1 (MQ=255)

                                                    tAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACt                                     >  1:4371881/1‑93 (MQ=255)

                                                     aGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTg                                    <  1:2082849/93‑1 (MQ=255)

                                                     aGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTg                                    <  1:8136610/93‑1 (MQ=255)

                                                     aGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTg                                    <  2:2120197/93‑1 (MQ=255)

                                                      ggCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTgg                                   >  2:1042164/1‑93 (MQ=255)

                                                       gCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGc                                  <  2:376443/93‑1 (MQ=255)

                                                        ccGCGTCAAACCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCa                                 <  1:2554804/93‑1 (MQ=255)

                                                         cgcgTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCAc                                >  1:7096543/1‑93 (MQ=255)

                                                                agccGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATc                         <  2:7067578/91‑1 (MQ=255)

                                                                aaCCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATc                         >  1:7440479/1‑93 (MQ=255)

                                                                aaCCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCGCTGGCATc                         >  1:5276838/1‑93 (MQ=255)

                                                                 aCCGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCa                        <  2:1531515/93‑1 (MQ=255)

                                                                  ccGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAg                       >  1:9625061/1‑93 (MQ=255)

                                                                  ccGCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAg                       <  2:6710184/93‑1 (MQ=255)

                                                                   cGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAgg                      >  2:411447/1‑93 (MQ=255)

                                                                    gCTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGt                     >  1:4042752/1‑93 (MQ=255)

                                                                     cTGGAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTg                    >  2:7275575/1‑93 (MQ=255)

                                                                       ggAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGcc                  >  2:2695368/1‑93 (MQ=255)

                                                                       ggAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGcc                  <  1:3621580/93‑1 (MQ=255)

                                                                        gAGTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGccc                 <  2:2632969/93‑1 (MQ=255)

                                                                        gAGTTACCGTTGCTCGACGGCATTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGccc                 <  1:1119065/93‑1 (MQ=255)

                                                                        gAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGccc                 <  1:2204365/93‑1 (MQ=255)

                                                                          gTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGa               >  2:8539090/1‑93 (MQ=255)

                                                                          gTTACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGa               >  1:358097/1‑93 (MQ=255)

                                                                           ttACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGAt              <  2:2349270/93‑1 (MQ=255)

                                                                            tACCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCAATGGGATCAGGTGCCCGATc             <  2:4069833/93‑1 (MQ=255)

                                                                             aCCGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCa            >  2:6133787/1‑93 (MQ=255)

                                                                               cGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCAcc          >  2:8364096/1‑93 (MQ=255)

                                                                               cGTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCAcc          >  1:1625990/1‑93 (MQ=255)

                                                                                gTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCt         >  2:6606462/1‑93 (MQ=255)

                                                                                gTTGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCt         >  2:1533888/1‑93 (MQ=255)

                                                                                 ttGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGTCCCATCACCTg        >  1:1431907/1‑93 (MQ=255)

                                                                                 ttGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCTg        >  1:6719765/1‑93 (MQ=255)

                                                                                  tGCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCTGa       >  1:4759497/1‑93 (MQ=255)

                                                                                   gCTCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCTGaa      >  1:2699720/1‑93 (MQ=255)

                                                                                     tCGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCTGaaaa    <  1:4071373/93‑1 (MQ=255)

                                                                                      cGACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCTGAAAAt   >  1:6021241/1‑93 (MQ=255)

                                                                                       gACGGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCTGAAAAtt  >  2:2445963/1‑93 (MQ=255)

                                                                                       gACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCTGAAAAtt  <  2:9281505/93‑1 (MQ=255)

AAACCGGTACGCCGTAATTTTGTACCCGCGCCAAACTATGCCGAAGCGTTTTTAGGCCGCGTCAAACCGCTGGAGTTACCGTTGCTCGACAGCCTTGAGCGCGAGTTACGGCGGATGACCGGCGTTACCGTTGACCGCGAAGACTGGCACTGGGATCAGGTGCCCGATCACCTGAAAATT  >  NC_000913/1483818‑1483997

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 17 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`