BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,458,272	(G)_5→4	coding (218/1311 nt)	b2342 ←	putative acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,458,268	0	G	.	99.3%	613.7 / ‑3.4	136	coding (222/1311 nt)	b2342	putative acyltransferase
Reads supporting (aligned to +/- strand): ref base G (1/0); new base . (74/61); total (75/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 150 total aligned reads displayed.

CCTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTA  >  NC_000913/2458177‑2458356

 cTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggg                                                                                        <  2:2196450/93‑1 (MQ=255)

 cTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggg                                                                                        <  2:4077189/93‑1 (MQ=255)

 cTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggg                                                                                        >  2:3639786/1‑93 (MQ=255)

 cTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggg                                                                                        >  1:4609624/1‑93 (MQ=255)

  tGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTgggg                                                                                       >  1:3500734/1‑93 (MQ=255)

  tGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTgggg                                                                                       >  1:2030354/1‑93 (MQ=255)

      aaCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggcttc                                                                                   <  1:2343755/93‑5 (MQ=255)

         tGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAgg                                                                               >  2:4161164/1‑93 (MQ=255)

          ggTAGCGCAAGCGCGGCTGCCGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGc                                                                              <  2:4317172/93‑1 (MQ=255)

             aGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCAtt                                                                           >  2:4665998/1‑93 (MQ=255)

              gcgcAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCAttt                                                                          <  1:487953/93‑1 (MQ=255)

              gcgcAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCAttt                                                                          <  1:3197601/93‑1 (MQ=255)

               cgcAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTg                                                                         <  2:4012982/93‑1 (MQ=255)

               cgcAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTg                                                                         <  2:1094815/93‑1 (MQ=255)

                gcAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGt                                                                        <  2:758700/93‑1 (MQ=255)

                  aaGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTAc                                                                      <  1:2173075/93‑1 (MQ=255)

                   aGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACg                                                                     >  2:1300566/1‑93 (MQ=255)

                   aGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACg                                                                     >  2:1024798/1‑93 (MQ=255)

                      gcgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGAcc                                                                  <  2:3316923/93‑1 (MQ=255)

                      gcgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGAcc                                                                  <  1:3412778/93‑1 (MQ=255)

                      gcgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGAcc                                                                  <  1:2959306/93‑1 (MQ=255)

                       cgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGAccg                                                                 >  2:853742/1‑92 (MQ=255)

                       cgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCt                                                                 >  1:1558328/1‑93 (MQ=255)

                        ggCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGTCATTTGTACGACCTg                                                                >  2:1050158/1‑93 (MQ=255)

                          cTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGAc                                                              >  2:2899901/1‑93 (MQ=255)

                            gACGCTGTAAGCATCGGTATGCACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGGT‑CGTGCTTCAGGCATTTGTACGACCTGACCa                                                            >  1:2895397/1‑93 (MQ=255)

                              cGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCaaa                                                          >  2:3460077/1‑93 (MQ=255)

                               ccTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAg                                                         <  2:1452765/92‑1 (MQ=255)

                                  gTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGAcc                                                      >  1:2865519/1‑93 (MQ=255)

                                   tAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCa                                                     >  1:4381494/1‑93 (MQ=255)

                                    aaGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAg                                                    <  2:1168127/93‑1 (MQ=255)

                                     aGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgt                                                   <  2:2758373/93‑1 (MQ=255)

                                     aGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgt                                                   <  2:2758374/93‑1 (MQ=255)

                                      gCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTGCACGCGCAATGTT‑CGGGCTTCAGGCATTTGTACGACCTGACAAAAGACCAgtt                                                  <  1:2408112/93‑1 (MQ=255)

                                       cATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgttg                                                 >  1:4608208/1‑93 (MQ=255)

                                       cATCGGTATGTACATTCATTCCCGTACCGAGAACAATGTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgttg                                                 <  1:3368457/93‑1 (MQ=255)

                                        aTCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgttgt                                                >  2:988604/1‑93 (MQ=255)

                                        aTCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgttgt                                                <  1:1393376/93‑1 (MQ=255)

                                         tCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑TGTGCTTCAGGCATTTGTACGACCTGACCATACACCAgttgtt                                               >  2:373481/1‑93 (MQ=255)

                                          cGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTc                                              >  2:2732329/1‑93 (MQ=255)

                                            gTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCaa                                            >  1:2548731/1‑93 (MQ=255)

                                            gTATGTACACTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCaa                                            >  1:3441786/1‑93 (MQ=255)

                                              aTGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATc                                          <  1:2375269/93‑1 (MQ=255)

                                               tGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCa                                         >  2:204525/1‑93 (MQ=255)

                                               tGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCa                                         >  1:1791933/1‑93 (MQ=255)

                                                gTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCAc                                        <  2:3071487/93‑1 (MQ=255)

                                                 tACAGTCATTCGCGTGCCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACt                                       <  1:44628/93‑1 (MQ=255)

                                                  aCATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGAATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACtt                                      <  1:1481188/93‑1 (MQ=255)

                                                   tattcattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTc                                     <  1:620519/92‑1 (MQ=255)

                                                   cattcattCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTc                                     <  1:3794231/93‑1 (MQ=255)

                                                     ttcattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcgg                                   >  1:4651734/1‑93 (MQ=255)

                                                     ttcattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcgg                                   >  2:380448/1‑93 (MQ=255)

                                                     ttcattcCCGTAACGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcgg                                   <  2:4291071/93‑1 (MQ=255)

                                                       cattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcggcg                                 >  1:1912275/1‑93 (MQ=255)

                                                        attcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcggcgg                                <  1:3264797/93‑1 (MQ=255)

                                                         ttcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCggg                               <  2:3326820/93‑1 (MQ=255)

                                                         ttcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCggg                               >  2:4528698/1‑93 (MQ=255)

                                                         ttcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCggg                               >  2:2258373/1‑93 (MQ=255)

                                                         ttcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCggg                               >  1:3152898/1‑93 (MQ=255)

                                                          tcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCgggg                              <  1:226269/93‑1 (MQ=255)

                                                          tcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCgggg                              <  1:4740045/93‑1 (MQ=255)

                                                          tcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCgggg                              >  2:2379742/1‑93 (MQ=255)

                                                           cccGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGa                             >  2:1383331/1‑93 (MQ=255)

                                                            ccGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAt                            >  1:2042246/1‑93 (MQ=255)

                                                              gTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtct                          <  1:1096297/93‑1 (MQ=255)

                                                               tACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtctc                         <  1:1253754/93‑1 (MQ=255)

                                                               tACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtctc                         <  2:2061120/93‑1 (MQ=255)

                                                               tACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtctc                         >  1:2920250/1‑93 (MQ=255)

                                                               tACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtctc                         <  2:3459795/93‑1 (MQ=255)

                                                                 ccGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGc                       <  1:260573/93‑1 (MQ=255)

                                                                  cGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCCAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGcc                      >  2:3709735/1‑92 (MQ=255)

                                                                  cGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGCTCTCGCt                      >  1:907760/1‑93 (MQ=255)

                                                                   gagaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGGTCTCGCTg                     >  1:4020790/1‑93 (MQ=255)

                                                                   gagaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTg                     >  2:1965398/1‑93 (MQ=255)

                                                                   gagaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTg                     >  2:4092439/1‑93 (MQ=255)

                                                                   gagaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGAGGGGATCTCGCTg                     >  2:1975004/1‑93 (MQ=255)

                                                                    agaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGc                    >  2:4531265/1‑93 (MQ=255)

                                                                    agaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGc                    >  2:1820575/1‑93 (MQ=255)

                                                                    agaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGc                    <  2:3509405/93‑1 (MQ=255)

                                                                     gaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCg                   <  2:2401293/93‑1 (MQ=255)

                                                                     gaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCg                   >  1:3763175/1‑93 (MQ=255)

                                                                      aaCAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGt                  <  1:541084/93‑1 (MQ=255)

                                                                      aaCAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGt                  <  1:3281978/93‑1 (MQ=255)

                                                                      aaCAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGt                  >  1:981403/1‑93 (MQ=255)

                                                                       aCAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTg                 >  2:3768117/1‑93 (MQ=255)

                                                                       aCAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTg                 >  1:1420964/1‑93 (MQ=255)

                                                                       aCAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTg                 >  2:3061341/1‑93 (MQ=255)

                                                                         aaTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGcc               >  1:2188546/1‑93 (MQ=255)

                                                                         aaTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGcc               <  2:515709/93‑1 (MQ=255)

                                                                          aTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCca              >  2:181275/1‑93 (MQ=255)

                                                                             tCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagca           <  2:2156125/93‑1 (MQ=255)

                                                                             tCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagca           <  2:2156126/93‑1 (MQ=255)

                                                                             tCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagca           <  2:2652000/93‑1 (MQ=255)

                                                                              cACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagcag          <  2:4679984/93‑1 (MQ=255)

                                                                                cgcgCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGtt        <  2:3120611/93‑1 (MQ=255)

                                                                                     aaTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCt   >  1:801760/1‑93 (MQ=255)

                                                                                     aaTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCt   >  1:933326/1‑93 (MQ=255)

                                                                                     aaTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCt   >  1:3728076/1‑93 (MQ=255)

                                                                                     aaTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCt   >  2:1248387/1‑93 (MQ=255)

                                                                                      aTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTa  >  2:1476793/1‑93 (MQ=255)

CCTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTA  >  NC_000913/2458177‑2458356

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