BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	302,564	A→G	Y117C (TAT→TGT)	yagU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	302,564	0	A	G	99.4%	513.5 / ‑5.0	164	Y117C (TAT→TGT)	yagU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (97/66); total (98/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 175 total aligned reads displayed.

CCAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTT  >  NC_000913/302474‑302653

ccAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGgttgtt                                                                                         >  2:4219831/1‑93 (MQ=255)

ccAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGgttgtt                                                                                         >  1:2002015/1‑93 (MQ=255)

ccAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGgttgtt                                                                                         >  1:1898611/1‑93 (MQ=255)

 cAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGttgttg                                                                                        <  2:1823804/93‑1 (MQ=255)

 cAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGttgttg                                                                                        >  1:1521858/1‑93 (MQ=255)

  aaTGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTtgttgt                                                                                       <  2:2810139/93‑1 (MQ=255)

   aTGCTGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTtgtg                                                                                      >  2:4145022/1‑93 (MQ=255)

    tGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTtgtgt                                                                                     <  2:3251200/93‑1 (MQ=255)

     gCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTtgtgtg                                                                                    >  1:1022644/1‑93 (MQ=255)

        gCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTc                                                                                 <  2:55478/93‑1 (MQ=255)

         cTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCg                                                                                >  1:1145655/1‑93 (MQ=255)

          tGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGc                                                                               >  2:3412897/1‑93 (MQ=255)

           gTTTATACGTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCt                                                                              >  2:4581490/1‑93 (MQ=255)

           gTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCt                                                                              >  2:3665610/1‑93 (MQ=255)

           gTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCt                                                                              >  1:4588822/1‑93 (MQ=255)

             ttATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGa                                                                            <  2:3804577/93‑1 (MQ=255)

              tataCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGaa                                                                           <  1:1837141/93‑1 (MQ=255)

              tataCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGaa                                                                           <  2:2341841/93‑1 (MQ=255)

               ataCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAg                                                                          <  1:4363630/93‑1 (MQ=255)

                taCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGt                                                                         <  2:2217994/93‑1 (MQ=255)

                 tccTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTa                                                                        <  1:2254222/92‑1 (MQ=255)

                 aCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTa                                                                        <  1:4034778/93‑1 (MQ=255)

                 aCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTa                                                                        >  2:2931166/1‑93 (MQ=255)

                   cTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAtt                                                                      >  2:190356/1‑93 (MQ=255)

                   cTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAtt                                                                      <  1:1180267/93‑1 (MQ=255)

                   cTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAtt                                                                      <  2:515815/93‑1 (MQ=255)

                    tttGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAttt                                                                     >  1:4134394/1‑93 (MQ=255)

                    tttGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAttt                                                                     <  2:3338339/93‑1 (MQ=255)

                       gCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCa                                                                  <  2:1069653/93‑1 (MQ=255)

                        ccGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaa                                                                 <  2:1807091/93‑1 (MQ=255)

                        ccGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaa                                                                 >  1:1772008/1‑93 (MQ=255)

                          gggCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaaaa                                                               >  2:3106275/1‑93 (MQ=255)

                            gCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAAt                                                             <  1:589480/93‑1 (MQ=255)

                              aTGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTa                                                           >  2:80365/1‑93 (MQ=255)

                              aTGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTa                                                           <  2:140214/93‑1 (MQ=255)

                               tGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTaa                                                          >  1:266193/1‑93 (MQ=255)

                               tGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTACAAAAATTaa                                                          >  2:3943657/1‑93 (MQ=255)

                                gTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTaaa                                                         >  2:2797806/1‑93 (MQ=255)

                                 tCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAAc                                                        >  2:2142234/1‑93 (MQ=255)

                                 tCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAAc                                                        <  1:1320811/93‑1 (MQ=255)

                                 tCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAAc                                                        <  2:2972283/93‑1 (MQ=255)

                                   tttAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAActc                                                      >  1:787665/1‑93 (MQ=255)

                                   tttAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAActc                                                      >  2:2078540/1‑93 (MQ=255)

                                   tttAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAActc                                                      <  1:3175345/93‑1 (MQ=255)

                                    ttAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAActct                                                     <  1:4204806/93‑1 (MQ=255)

                                    ttAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAActct                                                     <  2:4512305/93‑1 (MQ=255)

                                      aaCTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGTTGTCGGTTGTTGTGTGGTCGCTGAAGTCTTTCCAAAAATTAAACTCTgg                                                   <  2:2983679/93‑1 (MQ=255)

                                      aaCTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTgg                                                   <  1:3477065/93‑1 (MQ=255)

                                       aCTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGc                                                  <  1:273062/93‑1 (MQ=255)

                                       aCTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGc                                                  <  1:2308694/93‑1 (MQ=255)

                                         tGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAg                                                >  2:924343/1‑93 (MQ=255)

                                         tGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAg                                                >  1:2904120/1‑93 (MQ=255)

                                            gctGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTAGCAGGGc                                             >  2:1637457/3‑93 (MQ=255)

                                            gTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTTGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGc                                             >  1:2243910/1‑93 (MQ=255)

                                            gTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGc                                             <  1:4420732/93‑1 (MQ=255)

                                            gTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGc                                             <  1:2891098/93‑1 (MQ=255)

                                            gTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGc                                             >  1:3786085/1‑93 (MQ=255)

                                            gTTGGTGTTACGCACATTATCTTTTCGATAGTGTGTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATGTCCAAAAACTAAACTCTGGCAGGGc                                             >  1:4679661/1‑93 (MQ=255)

                                             ttGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCt                                            >  2:433170/1‑93 (MQ=255)

                                                 tgtTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACt                                        >  1:3955756/1‑93 (MQ=255)

                                                  gtTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTg                                       >  1:314821/1‑93 (MQ=255)

                                                     aCGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCa                                    >  2:794271/1‑93 (MQ=255)

                                                     aCGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCa                                    <  1:1369258/93‑1 (MQ=255)

                                                      cGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAg                                   <  2:232295/93‑1 (MQ=255)

                                                      cGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAg                                   <  2:3887757/93‑1 (MQ=255)

                                                      cGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAg                                   <  2:26767/93‑1 (MQ=255)

                                                      cGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAg                                   >  1:801346/1‑93 (MQ=255)

                                                       gCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAgg                                  >  1:737296/1‑93 (MQ=255)

                                                       gCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAgg                                  >  1:907720/1‑93 (MQ=255)

                                                       gCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGACAGGGCTTACTGGCAgg                                  >  1:1610601/1‑93 (MQ=255)

                                                        cacaTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAggg                                 >  1:3822801/1‑92 (MQ=255)

                                                        cacaTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAggg                                 >  2:2057079/1‑92 (MQ=255)

                                                         acaTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTg                                >  1:1806581/1‑93 (MQ=255)

                                                          caTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGc                               >  2:1583488/1‑93 (MQ=255)

                                                            ttATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCtt                             >  1:1380581/1‑93 (MQ=255)

                                                            ttATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCtt                             <  2:276404/93‑1 (MQ=255)

                                                             tATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCttt                            <  1:3853126/93‑1 (MQ=255)

                                                              aTCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTa                           >  2:1477549/1‑93 (MQ=255)

                                                               tCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg                          >  1:4280182/1‑93 (MQ=255)

                                                                  tttCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGccc                       <  2:423125/93‑1 (MQ=255)

                                                                   ttCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCa                      >  2:3782819/1‑93 (MQ=255)

                                                                      gATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACt                   >  2:1842214/1‑93 (MQ=255)

                                                                      gATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACt                   >  1:282556/1‑93 (MQ=255)

                                                                       aTAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACtt                  <  1:567441/93‑1 (MQ=255)

                                                                         aGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACtttt                >  2:2781201/1‑93 (MQ=255)

                                                                          gtgtTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACttttt               <  2:2548275/93‑1 (MQ=255)

                                                                          gtgtTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAAAttttt               >  1:1359772/1‑93 (MQ=255)

                                                                           tgtTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACtttttt              >  2:4098649/1‑92 (MQ=255)

                                                                           tgtTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTg              >  1:4474387/1‑93 (MQ=255)

                                                                            gtTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTGACCCAACTTTTTGt             >  1:125301/1‑93 (MQ=255)

                                                                             tttGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGtt            >  1:3014646/1‑93 (MQ=255)

                                                                             tttGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGtt            >  2:705318/1‑93 (MQ=255)

                                                                             tttGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTGtt            >  1:400215/1‑93 (MQ=255)

                                                                               tGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCa          <  1:2492110/93‑1 (MQ=255)

                                                                                gCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTGTTCat         >  2:4368178/1‑93 (MQ=255)

                                                                                 cTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTtata        >  2:1092802/1‑93 (MQ=255)

                                                                                  tGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTATCCCCACTTTTTGTTCatat       >  1:1675549/1‑93 (MQ=255)

                                                                                   gTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATg      <  2:3496831/93‑1 (MQ=255)

                                                                                   gTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACCTTTTTTTCATATg      >  1:4030521/1‑93 (MQ=255)

                                                                                      ggTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtt   >  2:2939166/1‑93 (MQ=255)

CCAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTT  >  NC_000913/302474‑302653

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