BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,075,220	(T)_5→6	coding (256/357 nt)	yggP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,075,215	1	.	T	99.4%	656.4 / ‑3.6	176	coding (261/357 nt)	yggP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base T (86/89); total (86/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

Only 100 of 198 total aligned reads displayed.

ATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑TTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAA  >  NC_000913/3075123‑3075302

ataAGTGGCATTGCCGCCAAAATTTTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑cc                                                                                           <  1:4274992/93‑1 (MQ=255)

 taAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑t                                                                                         >  1:1371869/1‑93 (MQ=255)

 taAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑t                                                                                         >  2:3687899/1‑93 (MQ=255)

 taAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑t                                                                                         >  1:3262450/1‑93 (MQ=255)

  aaGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑tt                                                                                        <  2:5512307/93‑1 (MQ=255)

  aaGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑tt                                                                                        >  1:5610532/1‑93 (MQ=255)

   aGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑ttt                                                                                       <  1:4573758/93‑1 (MQ=255)

      ggCATTGCCGCCACAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑tttttt                                                                                    <  2:5461912/93‑2 (MQ=255)

      ggCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑tttttt                                                                                    >  1:951466/1‑92 (MQ=255)

      ggCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑tttttt                                                                                    >  1:4902701/1‑92 (MQ=255)

       gCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑tttttta                                                                                   <  1:1007180/93‑3 (MQ=255)

         aTTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑ttttttata                                                                                 <  1:3513502/93‑5 (MQ=255)

           tGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATAtt                                                                               <  2:459561/93‑1 (MQ=255)

           tGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATAtt                                                                               <  2:6046151/93‑1 (MQ=255)

            gccgccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATAttt                                                                              <  1:5719157/93‑1 (MQ=255)

            gccgccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATAttt                                                                              <  1:2623825/93‑1 (MQ=255)

             ccgccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTg                                                                             <  2:5532733/93‑1 (MQ=255)

              cgccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGc                                                                            <  2:4026479/93‑1 (MQ=255)

               gccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGcc                                                                           <  1:381330/93‑1 (MQ=255)

                ccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCa                                                                          <  2:2921939/93‑1 (MQ=255)

                ccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCa                                                                          <  1:883472/93‑1 (MQ=255)

                 cAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAg                                                                         >  2:1421094/1‑93 (MQ=255)

                 cAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAg                                                                         <  1:291024/93‑1 (MQ=255)

                  aaaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGt                                                                        <  2:964668/93‑1 (MQ=255)

                  aaaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGt                                                                        <  2:4400427/93‑1 (MQ=255)

                   aaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTg                                                                       <  1:4064486/93‑1 (MQ=255)

                    aaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTga                                                                      >  1:5900180/1‑93 (MQ=255)

                    aaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTga                                                                      <  2:3496032/93‑1 (MQ=255)

                    aaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTga                                                                      >  2:4438382/1‑93 (MQ=255)

                     atatTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTgag                                                                     <  2:2491774/93‑1 (MQ=255)

                     atatTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTgag                                                                     <  1:5780789/93‑1 (MQ=255)

                      tatTCGTAGCTGTAGCCCGCTGAATATCTGCTTGGGAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTgaga                                                                    >  1:2833830/1‑93 (MQ=255)

                      tatTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTgaga                                                                    <  2:661966/93‑1 (MQ=255)

                      tatTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTgaga                                                                    >  2:2284862/1‑93 (MQ=255)

                          cGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGAtttt                                                                >  2:3302416/1‑93 (MQ=255)

                           gTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGAttttt                                                               <  2:3867721/93‑1 (MQ=255)

                            tAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTa                                                              >  1:2423830/1‑93 (MQ=255)

                              gCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTAcc                                                            <  1:2365914/93‑1 (MQ=255)

                              gCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTAcc                                                            >  1:5804279/1‑93 (MQ=255)

                                tGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCa                                                          >  2:1372657/1‑93 (MQ=255)

                                 gTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCAc                                                         >  2:2594876/1‑93 (MQ=255)

                                    gCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTc                                                      >  1:5240722/1‑93 (MQ=255)

                                     cccGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCg                                                     >  2:534652/1‑93 (MQ=255)

                                      ccGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGa                                                    <  2:1018417/93‑1 (MQ=255)

                                        gCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACa                                                  <  1:4736249/93‑1 (MQ=255)

                                         cTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACaa                                                 >  2:2590189/1‑93 (MQ=255)

                                         cTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACaa                                                 <  2:5118444/93‑1 (MQ=255)

                                          tGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAAt                                                >  2:4714408/1‑93 (MQ=255)

                                           gAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATa                                               <  1:1875853/93‑1 (MQ=255)

                                            aaTATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATaa                                              >  2:434983/1‑93 (MQ=255)

                                            aaTATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATaa                                              <  1:555569/93‑1 (MQ=255)

                                             atatCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAc                                             >  1:1500640/1‑93 (MQ=255)

                                             atatCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAc                                             <  1:958022/93‑1 (MQ=255)

                                              tatCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcc                                            >  2:1472047/1‑93 (MQ=255)

                                                 ccGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCg                                         >  2:5874154/1‑93 (MQ=255)

                                                  cGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCgg                                        >  2:4844011/1‑93 (MQ=255)

                                                   gCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGc                                       >  1:2794070/1‑93 (MQ=255)

                                                   gCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCCGc                                       >  1:5101295/1‑93 (MQ=255)

                                                    cTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATCACCCCGGca                                      >  2:1105981/1‑93 (MQ=255)

                                                    cTTGGTAACCCCATCGCCGGTTGCAAGACAAAACGCTGG‑CC‑TTTTTAAGATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGca                                      <  1:3634103/93‑1 (MQ=255)

                                                     ttGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGcac                                     >  2:2358942/1‑93 (MQ=255)

                                                        gTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACAtt                                  >  2:1585844/1‑93 (MQ=255)

                                                         tAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTc                                 >  2:2160485/1‑93 (MQ=255)

                                                          atccccATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCa                                <  1:583812/91‑1 (MQ=255)

                                                          aaCCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCa                                <  2:3199131/93‑1 (MQ=255)

                                                           aCCCCATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCAt                               <  2:4855636/93‑1 (MQ=255)

                                                             cccATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCCCATTCATGc                             >  2:4175179/1‑93 (MQ=255)

                                                             cccATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGc                             >  1:1123767/1‑93 (MQ=255)

                                                             cccATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGc                             <  1:3742789/93‑1 (MQ=255)

                                                             cccATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGc                             >  2:109945/1‑93 (MQ=255)

                                                               cATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGccc                           <  2:2730889/93‑1 (MQ=255)

                                                               cATCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGccc                           <  2:4956521/93‑1 (MQ=255)

                                                                aTCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCg                          <  2:1212698/93‑1 (MQ=255)

                                                                 tCGCCGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCgg                         >  2:908265/1‑93 (MQ=255)

                                                                    ccGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATGACCCCGGCACATTCATGCCCGGTaa                      <  1:3655043/93‑1 (MQ=255)

                                                                     cgattGCAATACAAAACGCTGGTCT‑TTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAAt                     <  1:1006166/90‑1 (MQ=255)

                                                                     cGGTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAAt                     <  1:4297127/93‑1 (MQ=255)

                                                                      ggTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATg                    >  2:6096957/1‑93 (MQ=255)

                                                                      ggTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATg                    >  1:5012462/1‑93 (MQ=255)

                                                                       gTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa                   >  2:3661795/1‑93 (MQ=255)

                                                                       gTTGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa                   >  1:5502437/1‑93 (MQ=255)

                                                                        ttGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGAc                  >  2:5660944/1‑93 (MQ=255)

                                                                         tGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGAcc                 >  1:1760146/1‑93 (MQ=255)

                                                                         tGCAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGAcc                 <  1:2964978/93‑1 (MQ=255)

                                                                           cAATACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCgg               >  2:3900324/1‑93 (MQ=255)

                                                                             aTACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCTGCACATTCATGCCCGGTAATGACCGgat             >  1:773678/1‑93 (MQ=255)

                                                                             aTACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGgat             >  1:3120735/1‑93 (MQ=255)

                                                                             aTACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGgat             >  1:4186971/1‑93 (MQ=255)

                                                                              tACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGgatg            >  1:1661094/1‑93 (MQ=255)

                                                                              tACAAAACGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGgatg            <  1:3449817/93‑1 (MQ=255)

                                                                                 aaaaCGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGAtt         >  1:663193/1‑93 (MQ=255)

                                                                                   aaCGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGAtttt       >  2:3072160/1‑93 (MQ=255)

                                                                                   aaCGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGAtttt       <  2:2645444/93‑1 (MQ=255)

                                                                                    aCGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTc      >  1:634095/1‑93 (MQ=255)

                                                                                     cGCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCt     <  1:2495833/93‑1 (MQ=255)

                                                                                      gCTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTa    <  1:1653399/93‑1 (MQ=255)

                                                                                       cTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTaa   <  1:2309152/93‑1 (MQ=255)

                                                                                       cTGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTaa   >  1:4366437/1‑93 (MQ=255)

                                                                                        tGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTaaa  >  1:4201502/1‑93 (MQ=255)

                                                                                        tGG‑CCTTTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTaaa  <  2:1796249/93‑1 (MQ=255)

ATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CC‑TTTTT‑ATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAA  >  NC_000913/3075123‑3075302

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