BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,926,900	T→G	15.2%	D647A (GAC→GCC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,926,900	0	T	G	15.2%	336.2 / 52.2	132	D647A (GAC→GCC)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base T (59/53); new base G (7/13); total (66/66)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 132 total aligned reads displayed.

GTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAA  >  NC_000913/1926804‑1926997

gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                <  1:3107537/100‑1 (MQ=255)

gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                <  1:2062233/100‑1 (MQ=255)

gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                <  1:1930597/100‑1 (MQ=255)

gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCgg                                                                                                <  1:721617/100‑1 (MQ=38)

gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCgg                                                                                                <  1:2747399/100‑1 (MQ=38)

gTTCCCTGCGCCAGCGCGACCAGAAAAGCATAGTCCATCGCTAGGCCGACGTACGATTTAAAGAGACCAGATTTACCGCCATGGCCTAAGTTCATGTCgg                                                                                                >  2:211606/1‑100 (MQ=255)

 ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGACTCAGGCCATGGCGGt                                                                                               <  1:2903305/100‑1 (MQ=255)

 ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGt                                                                                               <  2:2390909/100‑1 (MQ=38)

 ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGt                                                                                               <  1:2121376/100‑1 (MQ=38)

 ttCACTGCACCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGt                                                                                               <  1:60164/100‑1 (MQ=255)

  tCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                              >  1:2434524/1‑100 (MQ=255)

   cccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTAc                                                                                             >  1:2074768/1‑100 (MQ=255)

   cccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaa                                                                                             <  2:618282/100‑2 (MQ=37)

    ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaan                                                                                            <  2:2128675/100‑3 (MQ=37)

    ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaa                                                                                            <  2:2815138/100‑3 (MQ=37)

        cgcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGc                                                                                        <  1:354508/100‑1 (MQ=255)

         gcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGc                                                                                        <  2:1552204/99‑1 (MQ=255)

          ccAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCCTGGCCTGAGTCCATGTCGGTACagag                                                                                         >  2:3345424/1‑97 (MQ=255)

            aGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATa                                                                                    <  1:4152174/100‑1 (MQ=255)

               gcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                                 <  2:2813409/100‑1 (MQ=255)

                cgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAg                                                                                >  2:2571656/1‑100 (MQ=255)

                  aCCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGAt                                                                              <  1:926087/100‑1 (MQ=255)

                  aCCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGAt                                                                              <  1:1365252/100‑1 (MQ=255)

                    cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                                                            >  1:3257809/1‑100 (MQ=255)

                    cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                                                            <  2:2633519/100‑1 (MQ=255)

                     ataaaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                                                            <  1:1006788/97‑1 (MQ=255)

                     aGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGt                                                                           <  1:3068995/100‑1 (MQ=255)

                      gAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTc                                                                          >  1:4192262/1‑100 (MQ=255)

                      gAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTc                                                                          <  2:2829527/100‑1 (MQ=255)

                       aaaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCa                                                                         <  2:3010840/100‑1 (MQ=255)

                             acatTCCTGCGCTTCTCCTTCGTACGATTTAAAGTGACCACATTTACCTCCATGTTCTTAGTTCATGTCGGTACAGATCAATa                                                                                    <  2:2418121/81‑1 (MQ=255)

                               atTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGtttt                                                                 >  1:4009191/1‑100 (MQ=255)

                                ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                                                >  1:2869231/1‑100 (MQ=255)

                                ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                                                >  1:3811053/1‑100 (MQ=255)

                                ttCCATCGCTACGCCTTCGTACGATGTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                                                <  2:1850008/100‑1 (MQ=255)

                                 tCCATGGCTAGGCCTTCGTACCATTTAAATCGACCAGATTTACCTCTATGGCCTGAGTCCGTGGCGGTACAGAGCAATAAAAGATGGTCATCGGtttt                                                                 <  2:3868575/98‑1 (MQ=255)

                                 tCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                                                <  1:2561744/99‑1 (MQ=255)

                                 tCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                                               <  1:1194268/100‑1 (MQ=255)

                                 tCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                                               >  1:112554/1‑100 (MQ=255)

                                 tCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGAAATGGGTTTTCa                                                               >  2:233314/1‑100 (MQ=255)

                                  ccATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTGTCAg                                                              >  2:906725/1‑100 (MQ=255)

                                   cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                                             >  1:3855116/1‑100 (MQ=255)

                                   cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                                             >  1:2906024/1‑100 (MQ=255)

                                    aTCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAg                                                              <  2:2794407/98‑1 (MQ=255)

                                       gCTACGCCTTCGTACGATTTAAAGCGAGCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgc                                                         >  2:3527254/1‑100 (MQ=255)

                                        cTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgcg                                                        <  1:2057325/100‑1 (MQ=255)

                                        cTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgcg                                                        >  2:1782047/1‑100 (MQ=255)

                                        cTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgcg                                                        <  2:1630518/100‑1 (MQ=255)

                                           cGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCaa                                                     >  2:2219750/1‑100 (MQ=255)

                                            gCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAt                                                    <  1:1034682/100‑1 (MQ=255)

                                              cTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTAAGCTCACGCAAatt                                                  >  2:2405366/1‑100 (MQ=255)

                                                tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAg                                                >  1:1074397/1‑100 (MQ=255)

                                                tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAg                                                <  1:2477575/100‑1 (MQ=255)

                                                tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAg                                                <  1:3629395/100‑1 (MQ=255)

                                                 cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGc                                               <  1:2021390/100‑1 (MQ=255)

                                                  gTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg                                              >  1:733144/1‑100 (MQ=255)

                                                     cGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGAcc                                           >  1:626223/1‑100 (MQ=255)

                                                        tttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAt                                        >  2:3600634/1‑100 (MQ=255)

                                                            aaGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTg                                    <  1:840544/100‑1 (MQ=255)

                                                             aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                                   >  2:3880288/1‑100 (MQ=255)

                                                             aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                                   >  1:3965224/1‑100 (MQ=255)

                                                               cGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCg                                 <  2:1411951/100‑1 (MQ=255)

                                                                gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                                >  1:2090514/1‑100 (MQ=255)

                                                                gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                                <  1:1385287/100‑1 (MQ=255)

                                                                gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                                <  2:3334692/100‑1 (MQ=255)

                                                                gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                                >  1:3379625/1‑100 (MQ=255)

                                                                 ttgctctgtACCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGt                               >  2:3926328/9‑100 (MQ=21)

                                                                  ccAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGtt                              >  1:3561730/1‑100 (MQ=255)

                                                                   cAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATTGTCATCGGGTTTCAGCGCGCGCAATTGAGCGACCCATTTTGCCGGTTc                             >  2:839447/1‑100 (MQ=255)

                                                                    aGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTcc                            <  2:511411/100‑1 (MQ=255)

                                                                     gATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                           >  1:1407748/1‑100 (MQ=255)

                                                                       tttACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCaa                         <  2:2538972/100‑1 (MQ=255)

                                                                         tACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTCTGCCGGTTCCCAata                       >  1:4121363/1‑100 (MQ=255)

                                                                          aCCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat                      <  1:2712190/100‑1 (MQ=255)

                                                                          aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat                      >  1:3861908/1‑100 (MQ=38)

                                                                           ccGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTTAGCTCGCGCAATTTAGAGACCCATTTTGCAGGTTCTCAat                        >  2:4020862/1‑97 (MQ=255)

                                                                           ccGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt                     >  1:2620595/1‑100 (MQ=255)

                                                                           ccGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt                     >  1:4148155/1‑100 (MQ=255)

                                                                           ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt                     >  1:563204/1‑100 (MQ=38)

                                                                           ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt                     >  1:729492/1‑100 (MQ=38)

                                                                           ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt                     >  1:720591/1‑100 (MQ=38)

                                                                            cGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTg                    <  1:2583515/100‑1 (MQ=255)

                                                                             gCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc                   <  1:3210176/100‑1 (MQ=255)

                                                                                 tGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCt               >  1:1502340/1‑100 (MQ=255)

                                                                                  ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTg              <  1:1567860/100‑1 (MQ=255)

                                                                                   gCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTga             >  1:906380/1‑100 (MQ=255)

                                                                                    ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTgag            <  1:2883975/100‑1 (MQ=255)

                                                                                      tGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGaa          >  2:2076174/1‑100 (MQ=255)

                                                                                       gAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAAt         <  1:3049251/100‑1 (MQ=255)

                                                                                        aGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATc        <  2:3026153/100‑1 (MQ=255)

                                                                                        aGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATc        >  2:548899/1‑100 (MQ=255)

                                                                                        aGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATc        >  2:549875/1‑100 (MQ=255)

                                                                                          tCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTGGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGt      >  2:2331472/1‑100 (MQ=255)

                                                                                          tCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGt      >  2:4165809/1‑100 (MQ=255)

                                                                                          tCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCACGCAATTTAGCGACCCATTTTGCAGGTGCCCAATATTGCAACTGAGAATGGt      >  2:3879309/1‑100 (MQ=255)

                                                                                           ccATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTg     <  2:3224029/100‑1 (MQ=255)

                                                                                            cATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGc    <  2:1410129/100‑1 (MQ=255)

                                                                                            cATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATATAGCGACCCATTTTGCCGGATCCCAATACTGCACCTGAGAATCGTGc    >  2:3413363/1‑100 (MQ=255)

                                                                                             aTGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCa   <  1:2571663/100‑1 (MQ=255)

                                                                                               gTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCaa  <  2:2914874/99‑1 (MQ=255)

GTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAA  >  NC_000913/1926804‑1926997

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