BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	380,496	(C)_5→6	coding (16/444 nt)	b0359 ←	putative transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	380,491	1	.	C	98.6%	655.0 / 0.9	216	coding (21/444 nt)	b0359	putative transferase
Reads supporting (aligned to +/- strand): ref base . (2/1); new base C (91/122); total (93/123)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 241 total aligned reads displayed.

ATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑C‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAC  >  NC_000913/380399‑380572

aTAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTgg                                                                                     >  1:4332411/1‑93 (MQ=255)

aTAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTgg                                                                                     <  2:5931021/93‑1 (MQ=255)

aTAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTgg                                                                                     <  2:5146171/93‑1 (MQ=255)

 tAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑c                                                                                   >  1:6043231/1‑93 (MQ=255)

 tAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑c                                                                                   >  2:2182550/1‑93 (MQ=255)

 tAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑c                                                                                   <  2:41002/93‑1 (MQ=255)

   aCCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑c‑cc                                                                                >  1:2528915/1‑93 (MQ=255)

    ccGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑c‑ccc                                                                               <  2:3043307/93‑1 (MQ=255)

     cGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑c‑cccc                                                                              >  2:1253865/1‑93 (MQ=255)

           tttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAAt                                                                         >  1:106314/1‑93 (MQ=255)

                tctcGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCa                                                                    <  1:2253136/93‑1 (MQ=255)

                tctcGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCa                                                                    >  2:386966/1‑93 (MQ=255)

                  tcGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGa                                                                  >  2:625943/1‑93 (MQ=255)

                  tcGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGa                                                                  <  2:1834153/93‑1 (MQ=255)

                  tcGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGa                                                                  >  2:4980746/1‑93 (MQ=255)

                     gTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACaa                                                               >  1:467291/1‑93 (MQ=255)

                      tGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACCAc                                                              >  2:3251211/1‑93 (MQ=255)

                      tGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑CGCCCCTGAATCTCCAGACAAc                                                              <  2:6164854/93‑1 (MQ=255)

                       ggCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAAcc                                                             >  1:5009497/1‑93 (MQ=255)

                         cAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCaa                                                           <  1:2708033/93‑1 (MQ=255)

                         cAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCaa                                                           <  1:1889365/93‑1 (MQ=255)

                          aaaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑CACCCCTGAATCTCCAGACAACCAat                                                          >  1:5855227/1‑93 (MQ=255)

                           aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAata                                                         >  2:4748424/1‑93 (MQ=255)

                           aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAata                                                         >  2:2658850/1‑93 (MQ=255)

                           aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAata                                                         <  1:2680285/93‑1 (MQ=255)

                            aaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAatat                                                        <  2:3139475/93‑1 (MQ=255)

                            aaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAatat                                                        >  2:2780221/1‑93 (MQ=255)

                            aaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAatat                                                        <  2:6213695/93‑1 (MQ=255)

                            aaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAatat                                                        <  2:687622/93‑1 (MQ=255)

                             aCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATc                                                       >  1:5365294/1‑93 (MQ=255)

                             aCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATc                                                       <  2:561674/93‑1 (MQ=255)

                             aCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATc                                                       >  1:1534350/1‑93 (MQ=255)

                              ccAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCa                                                      <  2:1033456/93‑1 (MQ=255)

                              ccAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCa                                                      >  2:4199540/1‑93 (MQ=255)

                              ccAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCa                                                      >  2:89585/1‑93 (MQ=255)

                               cAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCAc                                                     >  2:2326280/1‑93 (MQ=255)

                                aaTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACt                                                    >  1:6202249/1‑93 (MQ=255)

                                 aTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACtt                                                   >  1:1199572/1‑93 (MQ=255)

                                 aTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACtt                                                   >  1:1744336/1‑93 (MQ=255)

                                 aTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACtt                                                   <  2:3228728/93‑1 (MQ=255)

                                  tttttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTa                                                  <  1:559163/93‑1 (MQ=255)

                                  tttttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTa                                                  >  2:4491999/1‑93 (MQ=255)

                                  tttttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTa                                                  >  2:2518110/1‑93 (MQ=255)

                                    tttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTaaa                                                <  1:6012420/93‑1 (MQ=255)

                                    tttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTaaa                                                <  1:925954/93‑1 (MQ=255)

                                       aCGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATaa                                             <  1:2755008/93‑1 (MQ=255)

                                       aCGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATaa                                             <  1:336958/93‑1 (MQ=255)

                                         gCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGt                                           >  1:5249723/1‑93 (MQ=255)

                                         gCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGt                                           <  1:746443/93‑1 (MQ=255)

                                          ccATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTg                                          >  2:3530547/1‑93 (MQ=255)

                                          ccATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTg                                          <  1:2382308/93‑1 (MQ=255)

                                           cATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGa                                         <  1:2940748/93‑1 (MQ=255)

                                            aTTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGAt                                        <  2:5839620/93‑1 (MQ=255)

                                             ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                       >  2:1375932/1‑93 (MQ=255)

                                             ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                       <  1:5084425/93‑1 (MQ=255)

                                             ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                       >  1:652172/1‑93 (MQ=255)

                                              tGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAg                                      >  1:6262584/1‑93 (MQ=255)

                                                cTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTc                                    >  1:4254891/1‑93 (MQ=255)

                                                cTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTc                                    <  1:3080359/93‑1 (MQ=255)

                                                 tGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCt                                   >  2:843704/1‑93 (MQ=255)

                                                  gATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                  >  1:3916031/1‑93 (MQ=255)

                                                  gATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                  >  1:4610797/1‑93 (MQ=255)

                                                    ttattGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                <  2:5365295/93‑1 (MQ=255)

                                                     tattGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                               <  1:5374110/93‑1 (MQ=255)

                                                     tattGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                               <  2:3778557/93‑1 (MQ=255)

                                                       ttGTGTCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAc                             <  1:4022137/93‑1 (MQ=255)

                                                       ttGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAc                             <  2:5632567/93‑1 (MQ=255)

                                                         gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                           <  1:4168038/93‑1 (MQ=255)

                                                           gCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGt                         >  2:1960859/1‑93 (MQ=255)

                                                           gCCAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGt                         >  2:3729762/1‑93 (MQ=255)

                                                            ccAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                        <  2:641601/93‑1 (MQ=255)

                                                             cAGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                       <  1:3006051/93‑1 (MQ=255)

                                                              aGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                      >  1:2916556/1‑93 (MQ=255)

                                                              aGGACGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                      >  1:3948696/1‑93 (MQ=255)

                                                                  cGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGa                  <  2:529941/93‑1 (MQ=255)

                                                                  cGAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGa                  >  2:2298395/1‑93 (MQ=255)

                                                                   gAATAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGAc                 <  1:1709774/93‑1 (MQ=255)

                                                                    aaTAAACGCATAATTACTAGTCTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACt                <  2:4784857/93‑1 (MQ=255)

                                                                     aTAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTa               <  1:1340172/93‑1 (MQ=255)

                                                                     aTAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTa               <  1:102894/93‑1 (MQ=255)

                                                                     aTAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTa               <  2:1680300/93‑1 (MQ=255)

                                                                     aTAAACGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTa               >  2:1759005/1‑93 (MQ=255)

                                                                       aaaCGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGt             <  2:4882158/93‑1 (MQ=255)

                                                                       aaaCGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGt             <  1:3859575/93‑1 (MQ=255)

                                                                        aaCGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTc            >  1:4711136/1‑93 (MQ=255)

                                                                          cGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAt          <  1:3736362/93‑1 (MQ=255)

                                                                          cGCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAt          <  1:6117600/93‑1 (MQ=255)

                                                                           gCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGCTAGTCTTAATCCTCGTTTTTAGACTCGTCAtt         <  2:5177810/93‑1 (MQ=14)

                                                                           gCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAtt         <  1:3590114/93‑1 (MQ=255)

                                                                           gCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAtt         >  2:2944845/1‑93 (MQ=255)

                                                                           gCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAtt         >  2:2100002/1‑93 (MQ=255)

                                                                           gCATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAtt         <  2:1129953/93‑1 (MQ=255)

                                                                            cATAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg        <  2:5978532/93‑1 (MQ=255)

                                                                             aTAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg       >  2:5351299/1‑93 (MQ=255)

                                                                             aTAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg       <  2:3977361/93‑1 (MQ=255)

                                                                              tAATTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGga      >  1:5201704/1‑93 (MQ=255)

                                                                               aaTTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag     <  1:1364048/93‑1 (MQ=255)

                                                                                aTTACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgaga    >  2:4271350/1‑93 (MQ=16)

                                                                                 ttACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa   <  1:2496599/93‑1 (MQ=12)

                                                                                  tACTAGACTGGCC‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc  <  2:3436004/93‑1 (MQ=9)

ATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑C‑CCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAC  >  NC_000913/380399‑380572

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