BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	821,422	A→G	S100S (AGT→AGC)	ybhN ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	821,422	0	A	G	96.0%	441.7 / 4.2	149	S100S (AGT→AGC)	ybhN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/3); new base G (71/72); total (74/75)
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.98e-01

Only 100 of 153 total aligned reads displayed.

ACTGAGCGAGAAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCGGGCGAG  >  NC_000913/821332‑821511

aCTGAGCGAGAAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCt                                                                                         <  2:4213835/93‑1 (MQ=255)

aCTGAGCGAGAAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACt                                                                                         <  2:4474120/93‑1 (MQ=255)

 cTGAGCGAGAAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTg                                                                                        <  2:292949/93‑1 (MQ=255)

   gAGCGAGAAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAg                                                                                      >  2:4266612/1‑93 (MQ=255)

   gAGCGAGAAAATCAGCGTAATAGTGCTGCCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATTCCGCCGGCCCAGGTGCTGAg                                                                                      <  2:4223459/93‑1 (MQ=255)

    aGCGAGAAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGc                                                                                     >  1:942694/1‑93 (MQ=255)

     gCGAGAAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCg                                                                                    <  1:86741/93‑1 (MQ=255)

       gagaAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTc                                                                                  >  2:1568490/1‑93 (MQ=255)

        agaAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCa                                                                                 >  1:2241117/1‑93 (MQ=255)

        agaAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCa                                                                                 >  2:4019300/1‑93 (MQ=255)

        agaAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCa                                                                                 >  1:1482129/1‑93 (MQ=255)

         gaAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAg                                                                                <  2:2093950/93‑1 (MQ=255)

            aaTCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGtt                                                                             <  2:2834863/93‑1 (MQ=255)

            aaTCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGtt                                                                             >  2:3554394/1‑93 (MQ=255)

             aTCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTg                                                                            >  2:2405999/1‑93 (MQ=255)

             aTCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTg                                                                            >  2:1801646/1‑93 (MQ=255)

              tCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGa                                                                           <  1:3767461/93‑1 (MQ=255)

              tCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGa                                                                           <  1:4364190/93‑1 (MQ=255)

              tCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGa                                                                           <  1:123195/93‑1 (MQ=255)

              tCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCCCCGACCCAGGTGCTGAGCGTCAGGTTGa                                                                           >  1:2095544/1‑93 (MQ=255)

                cgcgTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAg                                                                         >  2:2607327/1‑93 (MQ=255)

                cgcgTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAg                                                                         <  2:3027193/93‑1 (MQ=255)

                cgcgTAATAGTGCTACCGCGCAGCTACAGACGAGAGTACAAACGATGGCGCATACTAATGCGGCCGACCCAGGGGCTGAGCGTCAGGTTGAAg                                                                         <  2:295316/93‑1 (MQ=255)

                    tAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGt                                                                     >  2:2177978/1‑93 (MQ=255)

                    tAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGt                                                                     <  2:565382/93‑1 (MQ=255)

                     aaTAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTa                                                                    <  2:626937/93‑1 (MQ=255)

                      aTAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAg                                                                   <  2:3283298/93‑1 (MQ=255)

                       tAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGCTGAAGGCGTAga                                                                  >  2:3253409/1‑92 (MQ=255)

                        aGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATCCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCa                                                                 <  2:1531466/93‑1 (MQ=255)

                         gTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAg                                                                >  1:1652221/1‑93 (MQ=255)

                            cTACCGGGCAGCCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCTCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATa                                                             <  1:2523289/93‑1 (MQ=255)

                            cTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATTCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATa                                                             <  1:514409/93‑1 (MQ=255)

                            cTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATa                                                             <  2:2020080/93‑1 (MQ=255)

                            cTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATa                                                             <  1:1422004/93‑1 (MQ=255)

                            cTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGCTGAAGGCGTAGCAGATa                                                             >  1:786791/1‑93 (MQ=255)

                             tgccGGGCAACGCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATaa                                                            <  1:677963/91‑1 (MQ=255)

                             tACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATaa                                                            >  2:1256273/1‑93 (MQ=255)

                             tACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATaa                                                            <  1:3926185/93‑1 (MQ=255)

                              aCCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATaaa                                                           <  1:4199739/93‑1 (MQ=255)

                              aCCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATaaa                                                           <  2:2456841/93‑1 (MQ=255)

                                 gggCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGa                                                        <  1:2690764/93‑1 (MQ=255)

                                 gggCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGa                                                        <  1:722862/93‑1 (MQ=255)

                                  ggCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGac                                                       <  1:2316102/93‑1 (MQ=255)

                                  ggCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGac                                                       >  1:151356/1‑93 (MQ=255)

                                  ggCAACCCCAGACGAGAGGACAAAAGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACgcc                                                       >  2:3702586/1‑91 (MQ=255)

                                   gCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGaca                                                      >  1:1447806/1‑93 (MQ=255)

                                   gCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGaca                                                      <  1:2253830/93‑1 (MQ=255)

                                   gCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGaca                                                      >  1:1314683/1‑93 (MQ=255)

                                    cAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGacac                                                     <  1:1589104/93‑1 (MQ=255)

                                     aaCCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACAcc                                                    <  2:446283/93‑1 (MQ=255)

                                     aaCCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACAcc                                                    <  1:4064688/93‑1 (MQ=255)

                                     aaCCCCAGACGAGAGTACAAACGATAGAGCATACCAATGCCGCCGACCCAGGCGCTGAGCGTCAGGGTGAAGGCGTAGCACAGAAACGACAcc                                                    >  1:283804/1‑93 (MQ=255)

                                      aCCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCa                                                   <  2:1790819/93‑1 (MQ=255)

                                      aCCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCa                                                   >  1:3145659/1‑93 (MQ=255)

                                       ccccAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAg                                                  >  1:3114830/1‑93 (MQ=255)

                                         ccAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCa                                                <  2:2048948/93‑1 (MQ=255)

                                         ccAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCa                                                >  1:3465424/1‑93 (MQ=255)

                                         ccAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCa                                                <  1:104162/93‑1 (MQ=255)

                                          cAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACGCcatcat                                               >  2:109849/1‑93 (MQ=255)

                                          cAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCAt                                               <  1:2776532/93‑1 (MQ=255)

                                            gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCa                                             >  2:1746508/1‑93 (MQ=255)

                                            gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCa                                             >  1:2249581/1‑93 (MQ=255)

                                             aCGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCAc                                            >  2:709512/1‑93 (MQ=255)

                                              cGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCAcc                                           >  1:3036089/1‑93 (MQ=255)

                                              cGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCAcc                                           <  1:3530245/93‑1 (MQ=255)

                                              cGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCAcc                                           <  2:3081901/93‑1 (MQ=255)

                                                agagTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTg                                         <  2:3947431/93‑1 (MQ=255)

                                                  agTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGgc                                       <  2:2348138/93‑1 (MQ=255)

                                                    tACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGgcgc                                     >  1:1126329/1‑93 (MQ=255)

                                                       aaaCGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTc                                  >  1:2244895/1‑93 (MQ=255)

                                                       aaaCGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTc                                  >  2:4322005/1‑93 (MQ=255)

                                                       aaaCGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTc                                  <  1:3130513/93‑1 (MQ=255)

                                                        aaCGATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCg                                 >  2:3068287/1‑93 (MQ=255)

                                                           gATAGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCa                              >  2:2260881/1‑93 (MQ=255)

                                                              aGCGCATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCTCCTGGCGCTTCGCCAGtt                           <  2:2395839/93‑1 (MQ=255)

                                                               gcgcATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGttt                          <  2:1458654/93‑1 (MQ=255)

                                                               gcgcATACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGttt                          >  1:3624531/1‑93 (MQ=255)

                                                                   aTACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGa                      >  1:2310523/1‑93 (MQ=255)

                                                                   aTACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGa                      <  2:3157706/93‑1 (MQ=255)

                                                                    tACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGAc                     >  2:2045090/1‑93 (MQ=255)

                                                                    tACCAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGAc                     >  2:2029958/1‑93 (MQ=255)

                                                                      ccAATGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCg                   >  1:3658132/1‑93 (MQ=255)

                                                                        aaTGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCa                 <  1:623252/93‑1 (MQ=255)

                                                                        aaTGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCa                 <  1:1828262/93‑1 (MQ=255)

                                                                        aaTGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCa                 >  1:2531374/1‑93 (MQ=255)

                                                                         aTGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAg                <  1:898205/93‑1 (MQ=255)

                                                                          tGCCGCCGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGt               <  2:3879463/93‑1 (MQ=255)

                                                                           gccgccGACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTa              <  1:2699240/93‑1 (MQ=255)

                                                                              gccgACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTaaaa           <  1:1061447/93‑1 (MQ=255)

                                                                              cccgACCCACGTGATGAGCGTCCGCTTGAACGCGTAGCAGATAACCGACCCCAGCATCGCCTGGCGGTTCGCCAGTTTGTGACCGCAGTaaaa           <  2:3638197/92‑1 (MQ=255)

                                                                               ccgACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAg          <  1:1075126/93‑1 (MQ=255)

                                                                                 gACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCg        <  2:3902622/93‑1 (MQ=255)

                                                                                 gACCCAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCg        >  1:1659003/1‑93 (MQ=255)

                                                                                  aCCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCgg       <  1:563398/93‑1 (MQ=255)

                                                                                   cccAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCggg      >  1:1846219/1‑93 (MQ=255)

                                                                                    ccAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCGGGc     >  1:1423420/1‑93 (MQ=255)

                                                                                     cAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCGGGCg    >  1:2775230/1‑93 (MQ=255)

                                                                                     cAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCGGGCg    >  2:4623893/1‑93 (MQ=255)

                                                                                     cAGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCGGGCg    >  1:4427828/1‑93 (MQ=255)

                                                                                      aGGTGCTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCGGGCGa   <  2:4661849/93‑1 (MQ=255)

ACTGAGCGAGAAAATCCGCGTAATAGTGCTACCGGGCAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGACCGCAGTAAAAGCGGGCGAG  >  NC_000913/821332‑821511

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