BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,004	A→G	Y326C (TAT→TGT)	fruR →	FruR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,004	0	A	G	86.4%	630.2 / 78.1	236	Y326C (TAT→TGT)	fruR	FruR
Reads supporting (aligned to +/- strand): ref base A (12/20); new base G (116/88); total (128/108)
Fisher's exact test for biased strand distribution p-value = 5.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 249 total aligned reads displayed.

GTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACAT  >  NC_000913/88913‑89096

gTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTAt                                                                                             <  1:5679379/93‑1 (MQ=255)

 tCGCAGAGCGTGTGCTGTAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTc                                                                                            <  1:83671/93‑1 (MQ=255)

 tCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTc                                                                                            >  2:15325/1‑93 (MQ=255)

 tCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTc                                                                                            <  1:633760/93‑1 (MQ=255)

 tCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATc                                                                                            >  2:3188316/1‑93 (MQ=255)

 tCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATc                                                                                            <  2:118529/93‑1 (MQ=255)

   gCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTcgc                                                                                          <  2:3381068/93‑1 (MQ=255)

   cgagagCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTcgc                                                                                          <  1:5392011/91‑1 (MQ=255)

     agagCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTcgccg                                                                                        >  1:6097657/1‑93 (MQ=255)

      gagCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTcgccgc                                                                                       >  1:5229553/1‑93 (MQ=255)

      gagCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTcgccgc                                                                                       >  1:718878/1‑93 (MQ=255)

       agCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCcgcg                                                                                      >  1:4042928/1‑93 (MQ=255)

        gCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCgcgg                                                                                     <  1:1924730/93‑1 (MQ=255)

        gCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCgcgg                                                                                     <  1:2946414/93‑1 (MQ=255)

          gtgtGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCgcggcg                                                                                   >  1:6374615/1‑93 (MQ=255)

           tgtgCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGt                                                                                  >  1:3552162/1‑93 (MQ=255)

                tGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCa                                                                             >  2:1627184/1‑93 (MQ=255)

                tGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCa                                                                             >  2:1323560/1‑93 (MQ=255)

                 ggAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAg                                                                            >  1:5964802/1‑93 (MQ=255)

                 ggAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAg                                                                            <  2:4399409/93‑1 (MQ=255)

                 ggAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAg                                                                            >  2:2815062/1‑93 (MQ=255)

                  gagaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGc                                                                           <  1:1603905/93‑1 (MQ=255)

                  gagaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGc                                                                           <  1:1687170/93‑1 (MQ=255)

                    gaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCg                                                                         <  2:277750/93‑1 (MQ=255)

                     aTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGt                                                                        >  1:6605815/1‑93 (MQ=255)

                     aTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGt                                                                        >  1:4617683/1‑93 (MQ=255)

                     aTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGt                                                                        >  2:5388699/1‑93 (MQ=255)

                       tGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAg                                                                      <  1:4299808/93‑1 (MQ=255)

                       tGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAg                                                                      <  2:4106591/93‑1 (MQ=255)

                        gTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGc                                                                     <  2:6060619/93‑1 (MQ=255)

                         tCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCt                                                                    >  1:3425968/1‑93 (MQ=255)

                           cTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTaa                                                                  >  1:3214730/1‑93 (MQ=255)

                            tGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAg                                                                 >  2:1020984/1‑93 (MQ=255)

                               cAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCcg                                                              <  1:601785/93‑1 (MQ=255)

                                aaGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCcgc                                                             <  2:5058320/93‑1 (MQ=255)

                                aaGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCcgc                                                             <  2:5271466/93‑1 (MQ=255)

                                 aGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCcgcg                                                            >  1:4327083/1‑93 (MQ=255)

                                   ccTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGaa                                                          <  1:2104975/93‑1 (MQ=255)

                                       gACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACaaa                                                      <  2:4616238/93‑1 (MQ=255)

                                        aCGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACaaaa                                                     <  2:3968471/93‑1 (MQ=255)

                                         cGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACaaaaa                                                    >  1:3008675/1‑93 (MQ=255)

                                          gAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAAt                                                   <  2:2428045/93‑1 (MQ=255)

                                          gAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAAt                                                   <  2:6469319/93‑1 (MQ=255)

                                           aaCCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCAGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAACATa                                                  >  2:3019357/1‑93 (MQ=255)

                                            aCCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATAc                                                 >  1:5537248/1‑93 (MQ=255)

                                            aCCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATAc                                                 >  2:2023324/1‑93 (MQ=255)

                                             ccGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATAcg                                                >  2:6583371/1‑93 (MQ=255)

                                              cgcgTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATAcgc                                               <  1:4813190/93‑1 (MQ=255)

                                               gcgTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATAcgcg                                              >  2:4085470/1‑93 (MQ=255)

                                                cgTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACgcgc                                             >  2:2125011/1‑93 (MQ=255)

                                                  tAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCa                                           >  1:1786116/1‑93 (MQ=255)

                                                    aGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAgg                                         >  1:3933926/1‑93 (MQ=255)

                                                    aGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAgg                                         <  2:1833066/93‑1 (MQ=255)

                                                    aGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCACGAACAAAAATACGCGCCAgg                                         <  2:768160/93‑1 (MQ=255)

                                                     gCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGt                                        <  2:3430592/93‑1 (MQ=255)

                                                       cAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGa                                      >  2:4579366/1‑93 (MQ=255)

                                                         aaaCCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAAt                                    <  1:4974336/93‑1 (MQ=255)

                                                          aaCCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAAtt                                   >  1:1137736/1‑93 (MQ=255)

                                                           aCCTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAAttt                                  >  2:737953/1‑93 (MQ=255)

                                                             cTGGTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTcc                                >  2:1155659/1‑93 (MQ=255)

                                                                gTTTAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCctc                             >  1:6557141/1‑93 (MQ=255)

                                                                 tttAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCctct                            <  1:108288/93‑1 (MQ=255)

                                                                 tttAACGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCctct                            >  1:4469406/1‑93 (MQ=255)

                                                                   tAACGAGAATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTgg                          <  2:6506489/93‑1 (MQ=255)

                                                                    aaCGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGgc                         >  2:364180/1‑93 (MQ=255)

                                                                    aaCGCGCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGgc                         >  1:4398550/1‑93 (MQ=255)

                                                                     aCGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGgcg                        <  1:1331733/93‑1 (MQ=255)

                                                                      cgcgCATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGgcgc                       >  2:3550660/1‑93 (MQ=255)

                                                                       gcgcATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCACTCTGGcgcg                      >  2:5986848/1‑93 (MQ=255)

                                                                        cgcATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGt                     <  1:1314404/93‑1 (MQ=255)

                                                                        cgcATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGt                     <  1:2115949/93‑1 (MQ=255)

                                                                         gcATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGGGAATTTTCCTCTGGCGCGTa                    >  1:3675263/1‑93 (MQ=255)

                                                                          cATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTag                   <  2:1733311/93‑1 (MQ=255)

                                                                          cATTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTag                   >  2:5463864/1‑93 (MQ=255)

                                                                          cATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTag                   <  2:3857008/93‑1 (MQ=255)

                                                                           ataaaaCGTAATCTCTGCCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTaga                  <  1:102368/90‑1 (MQ=255)

                                                                           aTTAAACGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTaga                  >  1:2282401/1‑93 (MQ=255)

                                                                              aaaCGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTa               <  2:107891/93‑1 (MQ=255)

                                                                              aaaCGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTa               >  2:6345514/1‑93 (MQ=255)

                                                                              aaaCGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTa               <  2:493784/93‑1 (MQ=255)

                                                                              aaaCGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTa               <  1:3789146/93‑1 (MQ=255)

                                                                               aaCGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTAc              <  2:6414095/93‑1 (MQ=255)

                                                                               aaCGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTAc              >  1:6165153/1‑93 (MQ=255)

                                                                               aaCGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTAc              >  1:1390025/1‑93 (MQ=255)

                                                                                aCGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACg             >  2:420887/1‑93 (MQ=255)

                                                                                aCGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACg             >  2:2548087/1‑93 (MQ=255)

                                                                                 cGTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACgg            <  1:2601552/93‑1 (MQ=255)

                                                                                  gTAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACggg           <  2:4815634/93‑1 (MQ=255)

                                                                                   tAATCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGa          <  2:1215339/93‑1 (MQ=255)

                                                                                    aaTCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAc         <  1:4084591/93‑1 (MQ=255)

                                                                                     aTCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACt        <  1:589773/93‑1 (MQ=255)

                                                                                     aTCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACt        >  1:2486789/1‑93 (MQ=255)

                                                                                     aTCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACt        >  1:2713036/1‑93 (MQ=255)

                                                                                     aTCTCTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACt        >  1:6538011/1‑93 (MQ=255)

                                                                                      tctcTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTg       <  1:3296350/93‑1 (MQ=255)

                                                                                      tctcTGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTg       >  1:5786235/1‑93 (MQ=255)

                                                                                      tctcTGTCGCCGCGGCGTGCTCAGCCATAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTg       <  1:1528450/93‑1 (MQ=255)

                                                                                       ctctGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTgg      <  1:3169546/93‑1 (MQ=255)

                                                                                        tctGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGa     <  2:5877524/93‑1 (MQ=255)

                                                                                        tctGTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGa     <  1:5051091/93‑1 (MQ=255)

GTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACAT  >  NC_000913/88913‑89096

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