BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,975,807	A→G	L372P (CTG→CCG)	lysA ←	diaminopimelate decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,975,807	0	A	G	98.4%	587.2 / ‑1.2	189	L372P (CTG→CCG)	lysA	diaminopimelate decarboxylase
Reads supporting (aligned to +/- strand): ref base A (0/3); new base G (98/88); total (98/91)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

Only 100 of 193 total aligned reads displayed.

GCCTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCA  >  NC_000913/2975717‑2975896

 ccTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTa                                                                                        >  2:3623844/1‑93 (MQ=255)

 ccTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTa                                                                                        >  2:2020851/1‑93 (MQ=255)

  cTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAc                                                                                       >  2:62690/1‑93 (MQ=255)

   tGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAcc                                                                                      <  1:1418366/93‑1 (MQ=255)

   tGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAcc                                                                                      <  2:57108/93‑1 (MQ=255)

   tGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAca                                                                                      >  1:432763/1‑92 (MQ=255)

    gACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCa                                                                                     >  1:1682000/1‑93 (MQ=255)

     aCCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAg                                                                                    <  2:1676228/93‑1 (MQ=255)

         ttATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATaa                                                                                <  1:3067816/93‑1 (MQ=255)

         ttATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATaa                                                                                >  1:3225713/1‑93 (MQ=255)

         ttATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATaa                                                                                <  1:4679975/93‑1 (MQ=255)

         ttATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATaa                                                                                <  1:997429/93‑1 (MQ=255)

          tATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAAt                                                                               >  2:29223/1‑93 (MQ=255)

          tATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAAt                                                                               >  2:3578315/1‑93 (MQ=255)

          tATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAAt                                                                               >  2:2551724/1‑93 (MQ=255)

           aTCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATc                                                                              >  1:3896971/1‑93 (MQ=255)

            tCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCa                                                                             >  1:889585/1‑93 (MQ=255)

             cAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAc                                                                            <  2:2260877/93‑1 (MQ=255)

             cAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAc                                                                            >  1:2354732/1‑93 (MQ=255)

             cAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAc                                                                            >  2:2515959/1‑93 (MQ=255)

             cAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAc                                                                            <  1:2870565/93‑1 (MQ=255)

              aaaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAcc                                                                           <  2:4133205/93‑1 (MQ=255)

              aaaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAcc                                                                           >  1:3094115/1‑93 (MQ=255)

                 cAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGc                                                                        >  1:2979761/1‑93 (MQ=255)

                  aGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGcc                                                                       >  1:1250591/1‑93 (MQ=255)

                   gAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCt                                                                      >  2:1324024/1‑93 (MQ=255)

                   gAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAGTACCAGATAATCACCTGCCt                                                                      <  1:1126710/93‑1 (MQ=255)

                      cTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCa                                                                   >  1:1901451/1‑93 (MQ=255)

                      attCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCa                                                                   <  1:3835400/92‑1 (MQ=255)

                       ttCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCAc                                                                  <  1:1770708/93‑1 (MQ=255)

                       ttCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCAc                                                                  >  1:160999/1‑93 (MQ=255)

                        tCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACt                                                                 <  1:3508113/93‑1 (MQ=255)

                         cTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACtt                                                                <  2:3970097/93‑1 (MQ=255)

                          tGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTc                                                               >  1:3036171/1‑93 (MQ=255)

                          tGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTc                                                               >  1:668088/1‑93 (MQ=255)

                           ggTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTcc                                                              <  1:4623237/93‑1 (MQ=255)

                              aaCAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGc                                                           >  1:1336715/1‑93 (MQ=255)

                                 aGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAg                                                        <  1:419652/93‑1 (MQ=255)

                                  gCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAgg                                                       >  2:1486845/1‑93 (MQ=255)

                                   cGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgc                                                      <  2:1622557/93‑1 (MQ=255)

                                   cGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgc                                                      <  1:1877629/93‑1 (MQ=255)

                                   cGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgc                                                      >  2:4386291/1‑93 (MQ=255)

                                     gACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgcgc                                                    <  2:737366/93‑1 (MQ=255)

                                     gACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgcgc                                                    <  2:2655821/93‑1 (MQ=255)

                                      gcGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGcgcg                                                   <  1:2577867/92‑1 (MQ=255)

                                       cGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGTCAAGGCGCgg                                                  <  1:2733983/93‑1 (MQ=255)

                                       cGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCgg                                                  >  2:4307352/1‑93 (MQ=255)

                                       cGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCgg                                                  >  1:2813201/1‑93 (MQ=255)

                                        ggCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCggg                                                 <  1:3111137/93‑1 (MQ=255)

                                         gCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTCCCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGt                                                <  2:2598918/93‑1 (MQ=255)

                                         gCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGt                                                <  2:3713713/93‑1 (MQ=255)

                                           tATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGttt                                              >  2:2760818/1‑93 (MQ=255)

                                            aTTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTc                                             <  1:3257490/93‑1 (MQ=255)

                                            aTTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTc                                             >  2:4330700/1‑93 (MQ=255)

                                             ttGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCa                                            >  2:3905984/1‑93 (MQ=255)

                                             ttGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCa                                            >  1:2558610/1‑93 (MQ=255)

                                                   ttGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACAttt                                      <  1:2660499/93‑1 (MQ=255)

                                                    tGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTc                                     <  1:967836/93‑1 (MQ=255)

                                                      gATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTccc                                   <  2:776336/93‑1 (MQ=255)

                                                      gATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTccc                                   >  1:2340802/1‑93 (MQ=255)

                                                        tGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCGCCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTccccc                                 <  1:378429/93‑1 (MQ=255)

                                                        tGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTccccc                                 <  1:2070922/93‑1 (MQ=255)

                                                        tGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTccccc                                 >  2:2818132/1‑93 (MQ=255)

                                                        tGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTccccc                                 <  2:4283707/93‑1 (MQ=255)

                                                         gACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCt                                <  2:4484487/93‑1 (MQ=255)

                                                         gACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCt                                >  2:3079103/1‑93 (MQ=255)

                                                         gACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCt                                <  1:2447187/93‑1 (MQ=255)

                                                          aCATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCtt                               >  2:1067026/1‑93 (MQ=255)

                                                            aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             <  1:435026/93‑1 (MQ=255)

                                                            aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             <  1:4120970/93‑1 (MQ=255)

                                                            aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             <  1:3306794/93‑1 (MQ=255)

                                                            aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             <  1:758834/93‑1 (MQ=255)

                                                            aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             <  1:934815/93‑1 (MQ=255)

                                                            aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             <  2:349836/93‑1 (MQ=255)

                                                             ttGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCct                            >  2:2970977/1‑93 (MQ=255)

                                                               gATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCctgc                          >  1:3548630/1‑93 (MQ=255)

                                                                 tGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCctgctg                        <  1:856098/93‑1 (MQ=255)

                                                                   cgcCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTggg                      <  1:2629487/93‑1 (MQ=255)

                                                                     ccATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTa                    >  2:4133484/1‑93 (MQ=255)

                                                                      cATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTaa                   >  1:1453102/1‑93 (MQ=255)

                                                                       atatGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTaaa                  >  2:808347/1‑93 (MQ=255)

                                                                       atatGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTaaa                  >  1:1315029/1‑93 (MQ=255)

                                                                           gCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACa              <  2:3859960/93‑1 (MQ=255)

                                                                           gCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACa              >  1:730369/1‑93 (MQ=255)

                                                                            ccccTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACAt             >  2:4018164/1‑93 (MQ=255)

                                                                            ccccTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACAt             >  2:3426292/1‑93 (MQ=255)

                                                                            ccccTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACAt             >  1:4285772/1‑93 (MQ=255)

                                                                             cccTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATc            >  2:2090428/1‑93 (MQ=255)

                                                                              ccTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCg           >  2:4173018/1‑93 (MQ=255)

                                                                              ccTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCg           >  1:3748566/1‑93 (MQ=255)

                                                                               cTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGc          >  1:2261503/1‑93 (MQ=255)

                                                                               cTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGCTTCAACATTTCCCACTTCCTGCTGGGTAAAGACATCGc          >  2:465108/1‑93 (MQ=255)

                                                                                 gTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGccc        <  1:2701221/93‑1 (MQ=255)

                                                                                  tATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCg       >  1:982614/1‑93 (MQ=255)

                                                                                  tATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCg       >  2:4609547/1‑93 (MQ=255)

                                                                                    tCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGAt     >  2:1670105/1‑93 (MQ=255)

                                                                                    tCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGAt     >  1:1357230/1‑93 (MQ=255)

                                                                                     cATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGAtt    <  1:49336/93‑1 (MQ=255)

                                                                                      aTGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTc   >  2:1923954/1‑93 (MQ=255)

                                                                                       tGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTca  >  2:2426633/1‑93 (MQ=255)

GCCTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCA  >  NC_000913/2975717‑2975896

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