BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,929	A→G	E399G (GAG→GGG)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,929	0	A	G	99.3%	866.7 / ‑3.5	275	E399G (GAG→GGG)	thrC	threonine synthase
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (131/142); total (132/143)
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.46e-01

Only 100 of 279 total aligned reads displayed.

AATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGT  >  NC_000913/4839‑5018

aaTATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGc                                                                                         >  1:3855589/1‑93 (MQ=255)

aaTATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGc                                                                                         <  2:6731406/93‑1 (MQ=255)

aaTATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGc                                                                                         <  2:3724391/93‑1 (MQ=255)

 atatGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCt                                                                                        >  1:5043575/1‑93 (MQ=255)

 atatGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCt                                                                                        <  1:5814439/93‑1 (MQ=255)

   atGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTgg                                                                                      <  2:5650200/93‑1 (MQ=255)

      gCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAg                                                                                   >  2:4739743/1‑93 (MQ=255)

      gCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAg                                                                                   <  2:4073713/93‑1 (MQ=255)

        ttGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGaa                                                                                 >  2:3010303/1‑93 (MQ=255)

          gTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACg                                                                               <  1:6269885/93‑1 (MQ=255)

              cTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCt                                                                           >  1:5658342/1‑93 (MQ=255)

                cGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGa                                                                         <  2:5512175/93‑1 (MQ=255)

                  gCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGAtt                                                                       <  2:2292368/93‑1 (MQ=255)

                  gCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGAtt                                                                       <  1:2757/93‑1 (MQ=255)

                  gCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGAtt                                                                       <  2:6536700/93‑1 (MQ=255)

                  gCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGAtt                                                                       >  1:708476/1‑93 (MQ=255)

                   cACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGAttt                                                                      >  1:181653/1‑93 (MQ=255)

                   cACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGAttt                                                                      <  1:1226137/93‑1 (MQ=255)

                    aCCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTa                                                                     >  1:3298139/1‑93 (MQ=255)

                       gcgcATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTAccc                                                                  <  2:1748866/93‑1 (MQ=255)

                       gcgcATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTAccc                                                                  <  1:4901161/93‑1 (MQ=255)

                        cgcATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCt                                                                 >  1:3667689/1‑93 (MQ=255)

                         gcATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCtt                                                                <  1:6406457/93‑1 (MQ=255)

                         gcATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCtt                                                                >  1:3417690/1‑93 (MQ=255)

                          cATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTg                                                               >  1:6464594/1‑93 (MQ=255)

                          cATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTg                                                               >  1:1473416/1‑93 (MQ=255)

                           aTCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGc                                                              >  2:2486143/1‑93 (MQ=255)

                            tCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCt                                                             >  1:1389104/1‑93 (MQ=255)

                            tCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCt                                                             <  1:5216341/93‑1 (MQ=255)

                            tCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCt                                                             <  1:3936248/93‑1 (MQ=255)

                             ccGGTGAAATTTAAAGAGAGTGTGCAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCtt                                                            <  1:76397/93‑1 (MQ=255)

                                gCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTca                                                         <  2:438407/93‑1 (MQ=255)

                                  gAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTcaca                                                       >  2:3669520/1‑93 (MQ=255)

                                  gAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTcaca                                                       <  2:3133831/93‑1 (MQ=255)

                                   aaaTTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACAt                                                      <  2:806901/93‑1 (MQ=255)

                                    aaTTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATa                                                     >  1:1975992/1‑93 (MQ=255)

                                     aTTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATaa                                                    >  2:5413141/1‑93 (MQ=255)

                                      tttAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAAt                                                   <  2:6568427/93‑1 (MQ=255)

                                       ttAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATc                                                  <  2:4583495/93‑1 (MQ=255)

                                       ttAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTCGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATc                                                  <  2:2711631/93‑1 (MQ=255)

                                          aaGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGc                                               <  1:5755253/93‑1 (MQ=255)

                                          aaGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGc                                               >  2:5717234/1‑93 (MQ=255)

                                          aaGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGc                                               >  2:5022800/1‑93 (MQ=255)

                                          aaGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGc                                               >  1:6077490/1‑93 (MQ=255)

                                           agagagCGTGGAACCGATCCTCGGTGAACCGTAGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTccc                                              <  1:4592806/93‑1 (MQ=255)

                                            gagagCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGccc                                             >  1:3636930/1‑93 (MQ=255)

                                            gagagCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGccc                                             <  1:4187012/93‑1 (MQ=255)

                                            gagagCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGccc                                             >  1:4242121/1‑93 (MQ=255)

                                             agagCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCccg                                            <  2:3343519/93‑1 (MQ=255)

                                             agagCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCccg                                            <  2:1361306/93‑1 (MQ=255)

                                                gCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCccgccg                                         >  1:2047123/1‑93 (MQ=255)

                                                  gTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGAt                                       <  1:2958160/93‑1 (MQ=255)

                                                   tGGAAGCGATTCTCGGTGAAAGGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGAtt                                      <  2:2195004/93‑1 (MQ=255)

                                                   tGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGAtt                                      <  1:1003374/93‑1 (MQ=255)

                                                   tGGAAGCGATTCTCGGTGAAACGTTGGATCGGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGAtt                                      <  2:1158752/93‑1 (MQ=255)

                                                    ggAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGAttt                                     <  2:6025309/93‑1 (MQ=255)

                                                    ggAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGAttt                                     <  2:2604290/93‑1 (MQ=255)

                                                    ggAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGAttt                                     <  2:1952708/93‑1 (MQ=255)

                                                     gAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGAtttt                                    >  1:4506983/1‑93 (MQ=255)

                                                     gAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGAtttt                                    >  1:589611/1‑93 (MQ=255)

                                                       aGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTtgc                                  >  2:2893974/1‑93 (MQ=255)

                                                       aGCGATTCTCGGTGAAACGTGGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCGCATAATCTGCCCGCCGATTTtgc                                  <  1:5264229/93‑1 (MQ=255)

                                                          gATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTtgctgc                               >  2:5434444/1‑93 (MQ=255)

                                                          gATTCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTtgctgc                               >  1:1043073/1‑93 (MQ=255)

                                                            ttCTCGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCtgcgt                             >  1:5424307/1‑93 (MQ=255)

                                                              ctcGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCagcgttg                           <  1:5555252/93‑1 (MQ=255)

                                                                cGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCtgcgttgcg                         <  1:3509800/93‑1 (MQ=255)

                                                                cGGTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCtgcgttgcg                         <  1:5931327/93‑1 (MQ=255)

                                                                  gTGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTa                       >  1:4725861/1‑93 (MQ=255)

                                                                   tGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTaa                      <  2:6249121/93‑1 (MQ=255)

                                                                   tGAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTaa                      >  2:465157/1‑93 (MQ=255)

                                                                    gAAACGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTCCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTaaa                     <  1:2467703/93‑1 (MQ=255)

                                                                     aaaCGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAAt                    <  1:86001/93‑1 (MQ=255)

                                                                     aaaCGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGCTGCGTAAAt                    >  2:942244/1‑93 (MQ=255)

                                                                      aaCGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAAtt                   >  1:6169808/1‑93 (MQ=255)

                                                                        cGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtga                 <  2:1424712/93‑1 (MQ=255)

                                                                        cGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtga                 <  1:3089800/93‑1 (MQ=255)

                                                                        cGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtga                 >  1:4860434/1‑93 (MQ=255)

                                                                        cGTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtga                 <  2:3485625/93‑1 (MQ=255)

                                                                         gTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgat                <  1:5331291/93‑1 (MQ=255)

                                                                         gTTGGATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgat                <  2:5284444/93‑1 (MQ=255)

                                                                            ggATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgatgat             <  2:6588798/93‑1 (MQ=255)

                                                                             gATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTCCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgatgatg            <  1:3009773/93‑1 (MQ=255)

                                                                             gATCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgatgatg            >  1:6660138/1‑93 (MQ=255)

                                                                              aTCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgatgatga           <  2:5881473/93‑1 (MQ=255)

                                                                               tCTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGaa          <  1:538268/93‑1 (MQ=255)

                                                                                cTGCCAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAat         >  1:2285490/1‑93 (MQ=255)

                                                                                 tGCCAAAAGGGCTGGCAGAACGTGCTGATTTCCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatc        <  2:2861541/93‑1 (MQ=255)

                                                                                   ccAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatcat      >  2:6139425/1‑93 (MQ=255)

                                                                                   ccAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatcat      <  1:3528835/93‑1 (MQ=255)

                                                                                    cAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatcatc     >  2:5204291/1‑93 (MQ=255)

                                                                                    cAAAAGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatcatc     <  2:3818482/93‑1 (MQ=255)

                                                                                     aaaaGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAAtcatca    >  2:1981502/1‑93 (MQ=255)

                                                                                      aaaGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAg   >  1:2576803/1‑93 (MQ=255)

                                                                                      aaaGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAg   <  2:531159/93‑1 (MQ=255)

                                                                                      aaaGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAg   >  2:4543507/1‑93 (MQ=255)

                                                                                      aaaGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAg   >  1:980905/1‑93 (MQ=255)

                                                                                       aaGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGt  >  1:2275523/1‑93 (MQ=255)

                                                                                       aaGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGt  >  1:1518908/1‑93 (MQ=255)

                                                                                       aaGGGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGt  <  2:1726134/93‑1 (MQ=255)

AATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGT  >  NC_000913/4839‑5018

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