BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,758,577	G→T	N8K (AAC→AAA)	yfjH ←	CP4‑57 prophage; protein YfjH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,758,577	0	G	T	99.6%	2056.8 / 0.1	479	N8K (AAC→AAA)	yfjH	CP4‑57 prophage; protein YfjH
Reads supporting (aligned to +/- strand): ref base G (1/1); new base T (274/203); total (275/204)
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 491 total aligned reads displayed.

TTTCTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAATACTAC  >  NC_000913/2758479‑2758674

 ttCTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGtttt                                                                                                >  2:5011889/1‑101 (MQ=255)

  tCTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                               <  2:5011583/101‑1 (MQ=255)

  tCTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                               >  2:5011895/1‑101 (MQ=255)

  tCTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                               >  1:5011896/1‑101 (MQ=255)

    tACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttttt                                                                                             <  1:5011225/101‑1 (MQ=255)

        aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCt                                                                                         <  2:5011340/101‑1 (MQ=255)

        aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCt                                                                                         >  1:5011906/1‑101 (MQ=255)

         aTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCtt                                                                                        >  2:5011907/1‑101 (MQ=255)

            tGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTaaa                                                                                     >  2:5011914/1‑101 (MQ=255)

             gTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                    >  1:5011918/1‑101 (MQ=255)

             gTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                    <  2:5011487/101‑1 (MQ=255)

              ttCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGa                                                                                   <  1:5011157/101‑1 (MQ=255)

              ttCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGa                                                                                   >  2:5011924/1‑101 (MQ=255)

               tCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGAt                                                                                  >  1:5011928/1‑101 (MQ=255)

                cATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATa                                                                                 <  1:5011389/101‑1 (MQ=255)

                 aTAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATa                                                                                 <  1:5011406/100‑1 (MQ=255)

                 aTAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATAc                                                                                >  1:5011933/1‑101 (MQ=255)

                  tAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATa                                                                                 >  1:5011941/1‑99 (MQ=255)

                    aaaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCa                                                                             >  2:5011947/1‑101 (MQ=255)

                     aaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCAt                                                                            >  1:5011949/1‑101 (MQ=255)

                      aaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCAtt                                                                           <  2:5011435/101‑1 (MQ=255)

                       aTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                               <  1:5011321/80‑1 (MQ=255)

                       aTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATa                                                                                 >  1:5011954/1‑94 (MQ=255)

                          aTTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtc                                                                       >  2:5011961/1‑101 (MQ=255)

                           ttCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtct                                                                      >  1:5011964/1‑101 (MQ=255)

                             cAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctct                                                                    <  2:5011752/101‑1 (MQ=255)

                             cAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctct                                                                    >  2:5011969/1‑101 (MQ=255)

                               gagaAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATctctct                                                                  <  1:5010802/101‑1 (MQ=255)

                               gagaAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATctctct                                                                  <  1:5011559/101‑1 (MQ=255)

                                agaAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTg                                                                 >  2:5011976/1‑101 (MQ=255)

                                   aaTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCaa                                                              >  1:5011982/1‑101 (MQ=255)

                                    aTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAc                                                             <  2:5011024/101‑1 (MQ=255)

                                     tgaTAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                            >  1:5011987/3‑101 (MQ=255)

                                      cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCa                                                           >  2:5011988/1‑101 (MQ=255)

                                       aTAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCaa                                                          >  1:5011993/1‑101 (MQ=255)

                                        tAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCaaa                                                         <  1:5011350/101‑1 (MQ=255)

                                          aaaCCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAg                                                        >  2:5012005/1‑100 (MQ=255)

                                          aaaCCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGt                                                       >  1:5012003/1‑101 (MQ=255)

                                          aaaCCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGt                                                       <  1:5011649/101‑1 (MQ=255)

                                           aaCCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTg                                                      <  2:5011615/101‑1 (MQ=255)

                                            accaccATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGa                                                     >  2:5012011/1‑101 (MQ=255)

                                             ccaccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGaa                                                    <  2:5010922/101‑1 (MQ=255)

                                             ccaccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGaa                                                    <  1:5011028/101‑1 (MQ=255)

                                              caccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAc                                                   >  1:5012017/1‑101 (MQ=255)

                                              caccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAc                                                   <  2:5011254/101‑1 (MQ=255)

                                               accaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAcc                                                  >  2:5012020/1‑101 (MQ=255)

                                                ccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCa                                                 >  2:5012025/1‑101 (MQ=255)

                                                 caTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCaa                                                <  2:5011175/101‑1 (MQ=255)

                                                 caTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCaa                                                <  1:5011366/101‑1 (MQ=255)

                                                  aTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAAt                                               <  2:5011725/101‑1 (MQ=255)

                                                  aTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAAt                                               <  2:5011326/101‑1 (MQ=255)

                                                    aaaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATga                                             <  2:5011513/101‑1 (MQ=255)

                                                    aaaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATga                                             >  1:5012033/1‑101 (MQ=255)

                                                     aaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACAT‑GGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATgaga                                           >  2:5012034/1‑101 (MQ=255)

                                                      aTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATgaga                                           <  2:5011933/101‑1 (MQ=255)

                                                       tCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGagag                                          <  2:5011550/101‑1 (MQ=255)

                                                        cAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAgg                                         >  1:5012044/1‑101 (MQ=255)

                                                        cAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAgg                                         >  2:5012042/1‑101 (MQ=255)

                                                         aaTGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGc                                        <  2:5011457/101‑1 (MQ=255)

                                                           tGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGcaa                                      >  2:5012052/1‑101 (MQ=255)

                                                           tGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGcaa                                      >  1:5012053/1‑101 (MQ=255)

                                                             tATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGcaaca                                    <  2:5011601/101‑1 (MQ=255)

                                                              aTTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGcaacaa                                   >  2:5012063/1‑101 (MQ=255)

                                                              aTTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGcaacaa                                   >  1:5012061/1‑101 (MQ=255)

                                                               ttAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAg                                  >  1:5012068/1‑101 (MQ=255)

                                                                tAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGa                                 >  2:5012072/1‑101 (MQ=255)

                                                                 agagGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGaa                                <  1:5011598/101‑1 (MQ=255)

                                                                  gagGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAAt                               >  1:5012080/1‑101 (MQ=255)

                                                                   agGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATg                              >  1:5012087/1‑101 (MQ=255)

                                                                   agGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATg                              >  2:5012086/1‑101 (MQ=255)

                                                                    ggTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGa                             <  1:5011673/101‑1 (MQ=255)

                                                                     gTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGaa                            >  2:5012092/1‑101 (MQ=255)

                                                                         tAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAat                        <  1:5011479/101‑1 (MQ=255)

                                                                         tAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAat                        <  1:5011464/101‑1 (MQ=255)

                                                                         tAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAat                        <  2:5011363/101‑1 (MQ=255)

                                                                          aaaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatc                       >  2:5012101/1‑101 (MQ=255)

                                                                          aaaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatc                       >  1:5012102/1‑101 (MQ=255)

                                                                          aaaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatc                       <  2:5011573/101‑1 (MQ=255)

                                                                           aaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcg                      >  2:5012109/1‑101 (MQ=255)

                                                                           aaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcg                      <  1:5011503/101‑1 (MQ=255)

                                                                            aaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcga                     >  1:5012112/1‑101 (MQ=255)

                                                                             aaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcgat                    <  2:5011532/101‑1 (MQ=255)

                                                                             aaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcgat                    <  2:5011280/101‑1 (MQ=255)

                                                                              aTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcgatc                   >  1:5012120/1‑101 (MQ=255)

                                                                               ttGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcgatcg                  >  2:5012126/1‑101 (MQ=255)

                                                                               ttGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcgatcg                  <  2:5011402/101‑1 (MQ=255)

                                                                                 gTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTaa                                                                         >  1:5012131/1‑44 (MQ=255)

                                                                                  tATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGcc               <  1:5011546/101‑1 (MQ=255)

                                                                                   aTTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCg              >  2:5012134/1‑101 (MQ=255)

                                                                                     ttttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGtt            >  1:5012139/1‑101 (MQ=255)

                                                                                      tttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTa           <  1:5011444/101‑1 (MQ=255)

                                                                                      tttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTa           <  2:5011375/101‑1 (MQ=255)

                                                                                      tttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTa           >  1:5012145/1‑101 (MQ=255)

                                                                                       ttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAg          >  2:5012147/1‑101 (MQ=255)

                                                                                        tGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGa         >  1:5012154/1‑101 (MQ=255)

                                                                                         ggTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGaa        >  2:5012156/1‑101 (MQ=255)

                                                                                           tACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAta      <  1:5011294/101‑1 (MQ=255)

                                                                                           tACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAta      <  2:5011230/101‑1 (MQ=255)

                                                                                           tACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAta      >  1:5012163/1‑101 (MQ=255)

                                                                                            aCATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAtac     >  2:5012166/1‑101 (MQ=255)

TTTCTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAATACTAC  >  NC_000913/2758479‑2758674

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