BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,403,945	T→C	V131A (GTA→GCA)	accC →	acetyl CoA carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,403,945	0	T	C	99.6%	694.5 / ‑5.0	229	V131A (GTA→GCA)	accC	acetyl CoA carboxylase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (107/121); total (107/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

Only 100 of 246 total aligned reads displayed.

TCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCT  >  NC_000913/3403855‑3404034

tCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCAc                                                                                         >  2:1308219/1‑93 (MQ=255)

tCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCAc                                                                                         >  2:3804656/1‑93 (MQ=255)

  ttCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCg                                                                                       >  1:4226373/1‑93 (MQ=255)

  ttCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCg                                                                                       >  1:220786/1‑93 (MQ=255)

   tCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCgg                                                                                      <  1:1025809/93‑1 (MQ=255)

   tCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCgg                                                                                      >  2:1620412/1‑93 (MQ=255)

     aTTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGt                                                                                    >  1:3845198/1‑93 (MQ=255)

     aTTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGt                                                                                    <  2:720748/93‑1 (MQ=255)

     aTTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGt                                                                                    <  2:623871/93‑1 (MQ=255)

      ttGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGtt                                                                                   >  1:2449409/1‑93 (MQ=255)

      ttGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGtt                                                                                   >  1:2362739/1‑93 (MQ=255)

        ggCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCt                                                                                 <  1:3655145/93‑1 (MQ=255)

          cccGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGa                                                                               <  1:4105196/93‑1 (MQ=255)

           ccGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGAc                                                                              <  1:3744526/93‑1 (MQ=255)

            cGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACg                                                                             >  1:1950189/1‑93 (MQ=255)

            cGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACg                                                                             >  2:3377712/1‑93 (MQ=255)

             gAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACgg                                                                            >  1:3954336/1‑93 (MQ=255)

               aaGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGcc                                                                          >  1:2303995/1‑93 (MQ=255)

               aaGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGcc                                                                          >  1:716455/1‑93 (MQ=255)

               aaGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGcc                                                                          >  1:4553535/1‑93 (MQ=255)

               aaGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGcc                                                                          <  1:1387050/93‑1 (MQ=255)

                 gCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCg                                                                        >  2:1416582/1‑93 (MQ=255)

                 gCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCg                                                                        >  1:840580/1‑93 (MQ=255)

                  cAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGc                                                                       <  2:3770159/93‑1 (MQ=255)

                    gAAACCATTCGCCTGATGGGCGACAAAGTATGCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTg                                                                     >  2:786301/1‑93 (MQ=255)

                     aaaCCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTgg                                                                    >  1:1471466/1‑93 (MQ=255)

                     aaaCCATTCGCCTGATGGGCGACAAAGTATCCGCAAGCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTgg                                                                    <  2:2037636/93‑1 (MQ=255)

                      aaCCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTggg                                                                   <  2:2330493/93‑1 (MQ=255)

                      aaCCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTggg                                                                   <  2:293298/93‑1 (MQ=255)

                      aaCCATTCACCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTggg                                                                   >  2:4571878/1‑93 (MQ=255)

                       aCCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGc                                                                  >  2:3275407/1‑93 (MQ=255)

                        ccATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGcg                                                                 <  2:106389/93‑1 (MQ=255)

                         cATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGcga                                                                <  1:4278316/93‑1 (MQ=255)

                          aTTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGcgac                                                               <  2:4146811/93‑1 (MQ=255)

                           ttCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGcgacg                                                              >  2:86708/1‑93 (MQ=255)

                             cGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGat                                                            <  2:2685312/93‑1 (MQ=255)

                              gCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGata                                                           <  1:75202/93‑1 (MQ=255)

                              gCCTGATGGGCGACAAAGTATCCGCAAGCGCGGCGATGAAAAAAGCGGGCGTCCCTGGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGata                                                           <  2:3673170/93‑1 (MQ=255)

                               ccTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGatat                                                          <  2:2572054/93‑1 (MQ=255)

                               ccTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGatat                                                          <  2:1705195/93‑1 (MQ=255)

                                cTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATg                                                         >  1:790866/1‑93 (MQ=255)

                                cTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATg                                                         <  1:3246199/93‑1 (MQ=255)

                                cTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACAATATg                                                         <  2:574574/93‑1 (MQ=255)

                                 tGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATgg                                                        >  2:1997155/1‑93 (MQ=255)

                                  gATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGa                                                       <  1:1787683/93‑1 (MQ=255)

                                  gATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGa                                                       <  2:1026268/93‑1 (MQ=255)

                                  gATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGa                                                       <  1:2012036/93‑1 (MQ=255)

                                  gATGGGCGACAAAGTAGCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGa                                                       <  1:4122454/93‑1 (MQ=255)

                                    tGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATa                                                     >  1:4668437/1‑93 (MQ=255)

                                    tGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATa                                                     <  1:3186727/93‑1 (MQ=255)

                                     gggCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATaa                                                    <  2:1173315/93‑1 (MQ=255)

                                     gggCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATaa                                                    >  2:3930839/1‑93 (MQ=255)

                                      ggCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATaaa                                                   <  2:2308799/93‑1 (MQ=255)

                                       gCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATaaaa                                                  <  2:3636381/93‑1 (MQ=255)

                                        cGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATaaaaa                                                 <  1:1127075/93‑1 (MQ=255)

                                        cGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATaaaaa                                                 <  1:1137667/93‑1 (MQ=255)

                                         gACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCCTGCGCACCGGGGTCTGACCGCCCGCTGGGCGACGATATGGATAAAAAc                                                >  2:1502532/1‑93 (MQ=255)

                                          aCAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAAcc                                               <  1:592853/93‑1 (MQ=255)

                                           cAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCg                                              <  2:2470675/93‑1 (MQ=255)

                                            aaaGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGt                                             <  2:1490071/93‑1 (MQ=255)

                                            aaaGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGt                                             >  2:4072873/1‑93 (MQ=255)

                                            aaaGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGt                                             >  2:2849133/1‑93 (MQ=255)

                                              aGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGACCCTTGCGCACCGGGGTCTGACGGCCCGCTCGGCGTCGATACGGACAAGAACCGTGc                                           >  2:2451040/1‑93 (MQ=255)

                                                  tCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCAtt                                       <  2:421151/93‑1 (MQ=255)

                                                  tCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCAtt                                       >  1:304391/1‑93 (MQ=255)

                                                   ccGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTg                                      >  2:970402/1‑93 (MQ=255)

                                                   ccGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTg                                      >  1:3612501/1‑93 (MQ=255)

                                                     gCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCt                                    >  1:2169609/1‑93 (MQ=255)

                                                       aaTCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTaa                                  >  1:444348/1‑93 (MQ=255)

                                                       aaTCGCGGCGATGAAAAAAGCGGGCGCCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTaa                                  <  2:1561983/93‑1 (MQ=255)

                                                        aTCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTaaa                                 <  1:407907/93‑1 (MQ=255)

                                                        aTCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTaaa                                 <  1:2198005/93‑1 (MQ=255)

                                                        aTCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTaaa                                 >  1:2592165/1‑93 (MQ=255)

                                                         tCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAAc                                <  1:1885201/93‑1 (MQ=255)

                                                         tCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAAc                                <  1:3860039/93‑1 (MQ=255)

                                                             ggcgATGAAAAAAGCGGGCGTCCCTCGCGCGCCGGGTTCTGACGCCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCAt                            <  2:2884606/93‑1 (MQ=255)

                                                               cgATGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTg                          >  2:4012920/1‑93 (MQ=255)

                                                                 aTGAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGt                        <  2:3891927/93‑1 (MQ=255)

                                                                   gAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGGCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTa                      <  2:3055569/93‑1 (MQ=255)

                                                                   gAAAAAAGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTa                      <  1:2726623/93‑1 (MQ=255)

                                                                   aaaaaaaGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTa                      <  2:1814956/92‑1 (MQ=255)

                                                                      aaaaGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATcc                   <  1:706033/93‑1 (MQ=255)

                                                                       aaaTCGGGCGTCCCTTGCGCGCGGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCg                  <  1:3555139/93‑1 (MQ=255)

                                                                       aaaGCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCg                  >  1:4350498/1‑93 (MQ=255)

                                                                          gCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTg               >  1:524176/1‑93 (MQ=255)

                                                                          gCGGGCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTg               >  2:436438/1‑93 (MQ=255)

                                                                            gggCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGAt             >  1:3568771/1‑93 (MQ=255)

                                                                             ggCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGAtt            >  2:453860/1‑93 (MQ=255)

                                                                             ggCGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGAtt            <  2:2190953/93‑1 (MQ=255)

                                                                              gCGTCCCTTGCGCACCGGGTTCTGACGGCTCGGTGGGCGAAGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTa           <  1:3420298/93‑1 (MQ=255)

                                                                               cGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTAt          <  1:4602356/93‑1 (MQ=255)

                                                                               cGTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTAt          >  1:2231443/1‑93 (MQ=255)

                                                                                gTCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATc         <  2:2549927/93‑1 (MQ=255)

                                                                                 tCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCa        >  2:1965338/1‑93 (MQ=255)

                                                                                 tCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCa        >  2:2086566/1‑93 (MQ=255)

                                                                                 tCCCTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCa        <  1:376712/93‑1 (MQ=255)

                                                                                  cccTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCaa       <  2:27618/93‑1 (MQ=255)

                                                                                    cTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAg     <  1:3989070/93‑1 (MQ=255)

                                                                                    cTTGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAg     <  2:3347446/93‑1 (MQ=255)

                                                                                     ttGCGCACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGc    <  1:287156/93‑1 (MQ=255)

TCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCT  >  NC_000913/3403855‑3404034

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