BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,555	(G)_5→6	coding (55/618 nt)	b3646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,550	1	.	G	99.6%	712.7 / ‑3.9	222	coding (50/618 nt)	b3646	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (114/107); total (115/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

Only 100 of 255 total aligned reads displayed.

CACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGT  >  NC_000913/3816458‑3816637

   aaaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg                                                                                      <  1:5164823/93‑1 (MQ=255)

   aaaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg                                                                                      <  2:3602683/93‑1 (MQ=255)

     aTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggg                                                                                    >  1:4948164/1‑93 (MQ=255)

     aTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggg                                                                                    >  1:3883786/1‑93 (MQ=255)

     aTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggg                                                                                    <  2:6106380/93‑1 (MQ=255)

      tCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑gggggg                                                                                   >  1:3308096/1‑92 (MQ=255)

       cGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggggc                                                                                  <  2:4644492/93‑3 (MQ=255)

        gCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggggcg                                                                                 >  2:233704/1‑90 (MQ=255)

           aCCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTg                                                                               >  2:2954367/1‑93 (MQ=255)

            cccTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTgg                                                                              >  2:1437241/1‑93 (MQ=255)

             ccTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTggc                                                                             >  1:3145480/1‑93 (MQ=255)

             ccTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTggc                                                                             <  1:2325748/93‑1 (MQ=255)

              cTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTggcg                                                                            >  1:119258/1‑93 (MQ=255)

                 ccAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGcc                                                                         <  1:3653688/93‑1 (MQ=255)

                 ccAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGcc                                                                         <  1:4911025/93‑1 (MQ=255)

                  cagacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCg                                                                        <  1:2053263/93‑1 (MQ=255)

                  cagacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCg                                                                        >  2:5469521/1‑93 (MQ=255)

                  cagacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCg                                                                        >  1:2890754/1‑93 (MQ=255)

                   agacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCgg                                                                       >  2:1602506/1‑93 (MQ=255)

                    gacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGga                                                                      <  2:3177466/93‑1 (MQ=255)

                    gacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGga                                                                      <  1:5950116/93‑1 (MQ=255)

                    gacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGga                                                                      >  1:3016121/1‑93 (MQ=255)

                    gacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGga                                                                      <  1:2552664/93‑1 (MQ=255)

                    gacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGga                                                                      >  2:4198948/1‑93 (MQ=255)

                     acagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgac                                                                     >  1:1819545/1‑93 (MQ=255)

                      cagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacg                                                                    >  2:1203166/1‑93 (MQ=255)

                      cagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacg                                                                    <  1:465477/93‑1 (MQ=255)

                       agaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacga                                                                   <  1:2283091/93‑1 (MQ=255)

                       agaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacga                                                                   <  2:4924165/93‑1 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacgac                                                                  <  2:4063400/93‑1 (MQ=255)

                         aTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGacgacg                                                                 >  1:4884949/1‑93 (MQ=255)

                          tttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACgg                                                                <  1:5440818/93‑1 (MQ=255)

                          tttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACgg                                                                <  2:6192051/93‑1 (MQ=255)

                          tttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACgg                                                                <  2:639370/93‑1 (MQ=255)

                           ttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgc                                                               <  1:567960/93‑1 (MQ=255)

                           ttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgc                                                               <  2:4914260/93‑1 (MQ=255)

                           ttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgc                                                               >  1:3273627/1‑93 (MQ=255)

                               aGGTAGATAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCAt                                                           <  2:4756091/93‑1 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCag                                                           >  2:707735/1‑92 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCAt                                                           >  2:3309748/1‑93 (MQ=255)

                                 ggAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATgg                                                         >  2:390004/1‑93 (MQ=255)

                                   agagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGAt                                                       >  1:6231508/1‑93 (MQ=255)

                                    gagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATa                                                      >  1:2644826/1‑93 (MQ=255)

                                    gagaACCATGCTGTTACACATTTAGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATa                                                      >  2:5074840/1‑93 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    >  1:4739599/1‑93 (MQ=255)

                                       aaCCATGCTGTTCCACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACAt                                                   <  2:4407309/93‑1 (MQ=255)

                                        aCCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACAtt                                                  >  1:1069271/1‑93 (MQ=255)

                                        aCCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATCGATACAtt                                                  >  2:5635074/1‑93 (MQ=255)

                                          cATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTg                                                >  2:4315827/1‑93 (MQ=255)

                                           aTGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTgg                                               >  1:5515150/1‑93 (MQ=255)

                                           aTGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTgg                                               >  2:4944016/1‑93 (MQ=255)

                                           aTGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTgg                                               >  2:596936/1‑93 (MQ=255)

                                             gCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCCTGGATACATTTGGCg                                             >  2:486238/1‑93 (MQ=255)

                                             gCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGAGCTTGCGCCCGCGCGACGGCGCATGGATACATTTggtc                                             >  1:925071/1‑91 (MQ=255)

                                              cTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGt                                            >  1:1441162/1‑93 (MQ=255)

                                                gTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGACGGACGACGGCGCATGGATACATTTGGGGTaa                                          >  1:2225511/1‑93 (MQ=255)

                                                 ttACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAat                                         <  1:3100033/93‑1 (MQ=255)

                                                  tACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAatt                                        <  1:707931/93‑1 (MQ=255)

                                                   acacATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAatta                                       <  2:871228/93‑1 (MQ=255)

                                                   acacATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAatta                                       <  2:3534341/93‑1 (MQ=255)

                                                    cacaTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGAGCATGGATACATTTGGCGTAattat                                      >  1:4012266/1‑93 (MQ=255)

                                                     acaTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatt                                     <  1:1530976/93‑1 (MQ=255)

                                                      caTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatta                                    <  2:2377453/93‑1 (MQ=255)

                                                      caTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatta                                    <  1:4234072/93‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   <  1:2587713/93‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   >  2:2823383/1‑93 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   <  2:5576761/93‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   >  1:3660244/1‑93 (MQ=255)

                                                         tttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTg                                 <  2:5600965/93‑1 (MQ=255)

                                                         tttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGCAATTATTATTg                                 >  1:3597064/1‑93 (MQ=255)

                                                              aTTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAcc                            >  1:3946737/1‑93 (MQ=255)

                                                                ttggttggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGc                          >  1:4216252/1‑93 (MQ=255)

                                                                 tggttggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGcc                         >  2:2078977/1‑93 (MQ=255)

                                                                  ggttggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACCGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCa                        >  1:2501567/1‑93 (MQ=255)

                                                                   gttggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCAc                       <  1:1132945/93‑1 (MQ=255)

                                                                      ggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGc                    <  1:5017295/93‑1 (MQ=255)

                                                                      ggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGc                    >  2:1087507/1‑93 (MQ=255)

                                                                      ggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGc                    >  2:518379/1‑93 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  2:3440193/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  2:4777171/93‑1 (MQ=255)

                                                                          ttACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAtt                <  2:2468542/93‑1 (MQ=255)

                                                                           tACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTg               <  1:204129/93‑1 (MQ=255)

                                                                           tACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTg               >  1:2684971/1‑93 (MQ=255)

                                                                           tACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTg               <  1:2794482/93‑1 (MQ=255)

                                                                            aCTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTgg              <  1:271458/93‑1 (MQ=255)

                                                                            aCTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTgg              <  2:4667874/93‑1 (MQ=255)

                                                                              tGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCg            >  1:3780432/1‑93 (MQ=255)

                                                                              tGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCg            >  1:3833561/1‑93 (MQ=255)

                                                                              tGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCg            >  2:2573617/1‑93 (MQ=255)

                                                                               gCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCgg           <  1:1006604/93‑1 (MQ=255)

                                                                               gCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCgg           >  1:1833372/1‑93 (MQ=255)

                                                                                ccGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGa          >  2:4725351/1‑93 (MQ=255)

                                                                                ccGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGa          <  1:6210108/93‑1 (MQ=255)

                                                                                 cGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAg         <  1:823265/93‑1 (MQ=255)

                                                                                 cGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAg         >  2:5822988/1‑93 (MQ=255)

                                                                                   aaGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAggg       >  2:5303446/1‑93 (MQ=255)

                                                                                     gCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTc     <  2:1706542/93‑1 (MQ=255)

                                                                                     gCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTc     >  2:3757892/1‑93 (MQ=255)

                                                                                     gCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTc     >  2:1941228/1‑93 (MQ=255)

                                                                                       cATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGGGGGGCAg   >  1:4356791/1‑93 (MQ=255)

CACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGT  >  NC_000913/3816458‑3816637

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