BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,039,349	T→C	P77P (CCA→CCG)	mobA ←	MobA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,039,349	0	T	C	81.5%	618.2 / 149.4	303	P77P (CCA→CCG)	mobA	MobA
Reads supporting (aligned to +/- strand): ref base T (30/26); new base C (137/110); total (167/136)
Fisher's exact test for biased strand distribution p-value = 8.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 303 total aligned reads displayed.

AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGA  >  NC_000913/4039259‑4039440

aGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCgg                                                                                           <  1:61629/93‑1 (MQ=255)

 gggggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCggg                                                                                          <  1:4456181/93‑1 (MQ=255)

  ggggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGt                                                                                         <  1:6750956/93‑1 (MQ=255)

  ggggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGt                                                                                         <  1:6452979/93‑1 (MQ=255)

  ggggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGt                                                                                         >  2:3058580/1‑93 (MQ=255)

   gggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTa                                                                                        >  1:2546233/1‑93 (MQ=255)

      aaTGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATc                                                                                     >  2:6771005/1‑93 (MQ=255)

       aTGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATcc                                                                                    >  1:4541403/1‑93 (MQ=255)

         gTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGc                                                                                  <  1:6665795/93‑1 (MQ=255)

         gTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGc                                                                                  >  2:2847170/1‑93 (MQ=255)

          tAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGcc                                                                                 >  1:1091494/1‑93 (MQ=255)

             ggCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGt                                                                              <  1:2791113/93‑1 (MQ=255)

               cGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGa                                                                            >  2:122640/1‑93 (MQ=255)

                gTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGaa                                                                           <  1:851216/93‑1 (MQ=255)

                gTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGaa                                                                           >  1:5960624/1‑93 (MQ=255)

                gTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGaa                                                                           >  2:4701705/1‑93 (MQ=255)

                gTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGaa                                                                           >  1:4279676/1‑93 (MQ=255)

                  aTCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATc                                                                         >  2:5060050/1‑93 (MQ=255)

                  aTCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATc                                                                         <  2:3692684/93‑1 (MQ=255)

                  aTCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATc                                                                         <  2:5876361/93‑1 (MQ=255)

                   tCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCt                                                                        <  2:4086060/93‑1 (MQ=255)

                   tCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCt                                                                        >  1:5215008/1‑93 (MQ=255)

                     gCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTc                                                                      >  1:5141523/1‑93 (MQ=255)

                        cGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAAt                                                                   >  1:3820527/1‑93 (MQ=255)

                         gggCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATc                                                                  >  1:2034935/1‑93 (MQ=255)

                          ggCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCa                                                                 <  1:2689739/93‑1 (MQ=255)

                           gCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCAc                                                                <  1:2321475/93‑1 (MQ=255)

                           gCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCAc                                                                <  1:3930490/93‑1 (MQ=255)

                             aaaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACtt                                                              >  2:7481058/1‑93 (MQ=255)

                              aaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACttt                                                             <  1:4886925/93‑1 (MQ=255)

                              aaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACttt                                                             <  1:6964660/93‑1 (MQ=255)

                              aaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACttt                                                             <  1:6704019/93‑1 (MQ=255)

                                acaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCa                                                           <  2:1103466/93‑1 (MQ=255)

                                acaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCa                                                           <  2:1939371/93‑1 (MQ=255)

                                   aaaaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGAc                                                        <  1:3141307/93‑1 (MQ=255)

                                   aaaaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGAc                                                        >  2:2617012/1‑93 (MQ=255)

                                    aaaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGAcc                                                       >  1:6547542/1‑93 (MQ=255)

                                      aCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGc                                                     >  1:2625643/1‑93 (MQ=255)

                                      aCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGc                                                     >  2:4910653/1‑93 (MQ=255)

                                       ccACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCt                                                    <  2:138460/93‑1 (MQ=255)

                                        cACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCtt                                                   >  2:3753904/1‑93 (MQ=255)

                                        cACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCtt                                                   >  2:5503696/1‑93 (MQ=255)

                                        cACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCtt                                                   <  2:7124607/93‑1 (MQ=255)

                                         aCTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTg                                                  <  1:958156/93‑1 (MQ=255)

                                          cTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGc                                                 >  1:143858/1‑93 (MQ=255)

                                          cTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGc                                                 <  1:5303832/93‑1 (MQ=255)

                                                aGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATa                                           <  1:5782349/93‑1 (MQ=255)

                                                aGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATa                                           <  2:6671910/93‑1 (MQ=255)

                                                 gCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAg                                          >  2:3861209/1‑93 (MQ=255)

                                                  cTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGa                                         >  1:3044824/1‑93 (MQ=255)

                                                    tCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGTGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGAtt                                       <  1:1168394/93‑1 (MQ=255)

                                                    tCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGAtt                                       >  2:2647573/1‑93 (MQ=255)

                                                      ctgctgCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTc                                     <  2:4223278/93‑1 (MQ=255)

                                                      ctgctgCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTc                                     >  1:6835567/1‑93 (MQ=255)

                                                         ctgcATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCtg                                  >  1:1608532/1‑93 (MQ=255)

                                                         ctgcATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCtg                                  >  2:6268468/1‑93 (MQ=255)

                                                          tgcATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCtga                                 <  2:6007857/93‑1 (MQ=255)

                                                          tgcATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCtga                                 >  2:7218303/1‑93 (MQ=255)

                                                              ttACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGAc                             >  2:1794384/1‑93 (MQ=255)

                                                              ttACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGAc                             >  1:180110/1‑93 (MQ=255)

                                                               tACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACg                            >  1:5465859/1‑93 (MQ=255)

                                                                aCTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGa                           <  2:2456498/93‑1 (MQ=255)

                                                                aCTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGa                           <  2:7393639/93‑1 (MQ=255)

                                                                 cTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGAt                          >  2:2310588/1‑93 (MQ=255)

                                                                  tGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGAtt                         >  2:3644729/1‑93 (MQ=255)

                                                                   gAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTa                        >  2:2091230/1‑93 (MQ=255)

                                                                   gAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTa                        >  1:4063124/1‑93 (MQ=255)

                                                                   gAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTa                        >  1:703716/1‑93 (MQ=255)

                                                                    aaaGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAg                       >  1:7498213/1‑93 (MQ=255)

                                                                    aaaGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAg                       >  2:5416703/1‑93 (MQ=255)

                                                                     aaGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGc                      >  1:1260238/1‑93 (MQ=255)

                                                                      aGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCa                     <  1:4635242/93‑1 (MQ=255)

                                                                       gCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAt                    <  2:2336232/93‑1 (MQ=255)

                                                                        cATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAtt                   >  1:2468456/1‑93 (MQ=255)

                                                                        cATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAtt                   <  2:7335754/93‑1 (MQ=255)

                                                                         aTTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTa                  >  2:1655783/1‑93 (MQ=255)

                                                                          ttCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTaa                 >  1:5191117/1‑93 (MQ=255)

                                                                          ttCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTaa                 >  2:4137230/1‑93 (MQ=255)

                                                                           tCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAAc                <  1:2902783/93‑1 (MQ=255)

                                                                             cTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCa              >  2:1515965/1‑93 (MQ=255)

                                                                             cTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCa              <  1:3754367/93‑1 (MQ=255)

                                                                             cTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCa              >  1:3494091/1‑93 (MQ=255)

                                                                              tGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCac             >  1:4718807/1‑93 (MQ=255)

                                                                              tGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCac             >  2:865842/1‑93 (MQ=255)

                                                                                ccAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacga           <  2:1468805/93‑1 (MQ=255)

                                                                                ccAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacga           <  2:3959126/93‑1 (MQ=255)

                                                                                 cAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgac          <  2:4849734/93‑1 (MQ=255)

                                                                                 cAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgac          <  1:5727192/93‑1 (MQ=255)

                                                                                 cAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgac          >  2:1711895/1‑93 (MQ=255)

                                                                                  agagGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgacg         >  2:6813781/1‑93 (MQ=255)

                                                                                  agagGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgacg         >  1:5628399/1‑93 (MQ=255)

                                                                                   gagGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGt        <  2:4263278/93‑1 (MQ=255)

                                                                                   gagGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGt        >  2:5280481/1‑93 (MQ=255)

                                                                                   gagGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGt        >  1:2171972/1‑93 (MQ=255)

                                                                                    agGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTg       <  1:7310904/93‑1 (MQ=255)

                                                                                    agGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTg       >  2:643835/1‑93 (MQ=255)

                                                                                     gggCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTga      <  2:3371405/93‑1 (MQ=255)

                                                                                     gggCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTga      >  2:3273666/1‑93 (MQ=255)

                                                                                      ggCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTgag     >  1:3579590/1‑93 (MQ=255)

                                                                                      ggCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTgag     <  2:2687972/93‑1 (MQ=255)

AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGA  >  NC_000913/4039259‑4039440

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