BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,335,654	T→A	15.7%	intergenic (‑64/+40)	eptA ← / ← adiC	lipid A phosphoethanolamine transferase/arginine:agmatine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,335,654	0	T	A	15.7%	450.2 / 67.9	173	intergenic (‑64/+40)	eptA/adiC	lipid A phosphoethanolamine transferase/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand): ref base T (54/90); new base A (14/13); total (70/103)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

Only 100 of 175 total aligned reads displayed.

ATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAT  >  NC_000913/4335555‑4335752

   cAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTc                                                                                                     >  2:1100946/1‑99 (MQ=255)

     aaGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTa                                                                                                  <  2:3016881/100‑1 (MQ=255)

      agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa                                                                                                 <  2:1315717/100‑1 (MQ=255)

      agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCAAACAAAGTGCGCATATGCTCGAAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa                                                                                                 >  2:665323/1‑100 (MQ=255)

         gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGt                                                                                              >  2:1739319/1‑100 (MQ=255)

          ggTCTTTTTAGTAGGCGCTTCAACCTCACGGTGTTTCCATCGAACACAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGtt                                                                                             <  1:2508958/100‑1 (MQ=255)

           gTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTg                                                                                            <  1:3545186/100‑1 (MQ=255)

           gTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTg                                                                                            <  2:2853499/100‑1 (MQ=255)

             cTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGc                                                                                          <  2:2594602/100‑1 (MQ=255)

             cTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTg                                                                                            <  2:3439428/98‑1 (MQ=255)

              tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCt                                                                                         <  2:3639181/100‑1 (MQ=255)

              tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGCt                                                                                         <  2:2600546/100‑1 (MQ=255)

              ttgttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGCt                                                                                         <  2:232929/97‑1 (MQ=255)

               ttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCt                                                                                         >  1:2180959/1‑99 (MQ=255)

               ttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGCtt                                                                                        <  2:3257507/100‑1 (MQ=255)

                ttcggtacGCGTTTCAATATCACGGTGTGTCCATACAACAAAGTTCGTATATTTTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGCtta                                                                                       <  1:3789233/96‑2 (MQ=16)

                  tagtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTa                                                                                   <  2:2134466/100‑1 (MQ=255)

                    gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATg                                                                                 >  1:2684873/1‑100 (MQ=255)

                    gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATg                                                                                 <  2:272387/100‑1 (MQ=255)

                    gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATg                                                                                 >  2:3297458/1‑100 (MQ=255)

                     tagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑t                                                                               >  1:430407/1‑100 (MQ=255)

                      agGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑tt                                                                              <  2:3228936/100‑1 (MQ=255)

                       ggCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑ttt                                                                             <  2:1360727/100‑1 (MQ=255)

                         cgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGc                                                                           <  1:3407938/100‑1 (MQ=255)

                         cgcTTCAACATCACGGTGTTTCCATAGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑ttt                                                                             >  1:2985452/1‑98 (MQ=255)

                             tCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGa                                                                       <  1:741105/100‑1 (MQ=255)

                              cAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTgg                                                                        >  2:229837/1‑98 (MQ=255)

                              cAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGAt                                                                      <  2:373701/100‑1 (MQ=255)

                               aaCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGAtt                                                                     >  2:454832/1‑100 (MQ=255)

                                aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTa                                                                    <  1:3111282/100‑1 (MQ=255)

                                aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAACTTAAGGTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTa                                                                    >  1:2721759/1‑100 (MQ=255)

                                  aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAg                                                                  <  2:1739259/100‑1 (MQ=255)

                                    cACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAggg                                                                <  1:90544/100‑1 (MQ=255)

                                     aCGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAgggg                                                               <  2:2202461/100‑1 (MQ=255)

                                        gtgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAAt                                                            <  1:3661879/100‑1 (MQ=255)

                                         tgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAAtt                                                           <  2:705393/100‑1 (MQ=255)

                                          gtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTa                                                          <  1:1547844/100‑1 (MQ=255)

                                          gtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTa                                                          >  2:192304/1‑100 (MQ=255)

                                             tCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATc                                                       <  2:3142149/100‑1 (MQ=255)

                                               cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCtt                                                     >  1:3301364/1‑100 (MQ=255)

                                               cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCtt                                                     <  2:29796/100‑1 (MQ=255)

                                                   gAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCt                                                 <  1:1692255/100‑1 (MQ=255)

                                                       aaaGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtt                                             <  2:2684893/100‑1 (MQ=255)

                                                       aaaGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGAAATTAATCTTTGTTTAtt                                             >  1:1040004/1‑100 (MQ=255)

                                                        aaGTGCGCATATGCTCGCAAAGCGAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                            <  1:505314/100‑1 (MQ=255)

                                                         aGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTgg                                           >  1:410829/1‑100 (MQ=255)

                                                         aGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                            <  1:332200/99‑1 (MQ=255)

                                                           tGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                         <  1:1716397/100‑1 (MQ=255)

                                                            gcgcATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGc                                        <  2:593600/100‑1 (MQ=255)

                                                            gcgcATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGc                                        <  1:3433332/100‑1 (MQ=255)

                                                             cgcATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCa                                       <  2:1898950/100‑1 (MQ=255)

                                                             cgcATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCa                                       <  2:1762301/100‑1 (MQ=255)

                                                               cATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATc                                     >  1:1346477/1‑100 (MQ=255)

                                                                atatGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAATGTTAGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCt                                    <  1:1079382/100‑1 (MQ=255)

                                                                 tatGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTa                                   <  1:1257440/100‑1 (MQ=255)

                                                                    gCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAgg                                <  2:3208234/100‑1 (MQ=255)

                                                                    gCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAgg                                <  2:2735211/100‑1 (MQ=255)

                                                                     cTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAggg                               <  2:1551155/100‑1 (MQ=255)

                                                                      tCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGa                              >  2:1373155/1‑100 (MQ=255)

                                                                       cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGAt                             >  2:1445472/1‑100 (MQ=255)

                                                                       cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGAt                             >  2:3315164/1‑100 (MQ=255)

                                                                          aaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACg                          <  2:385363/100‑1 (MQ=255)

                                                                           aagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACgg                         >  2:1224087/1‑100 (MQ=255)

                                                                            agcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACggg                        <  1:21240/100‑1 (MQ=255)

                                                                            agcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACggg                        <  2:3095816/100‑1 (MQ=255)

                                                                              caaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGtt                      >  2:3678985/1‑100 (MQ=255)

                                                                                 ataaaaGCCAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttttt                   >  2:936910/5‑100 (MQ=38)

                                                                                 aaTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttttt                   >  2:1069026/1‑100 (MQ=255)

                                                                                 aaTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttttt                   <  1:1551271/100‑1 (MQ=255)

                                                                                  aTTAAGCCAGCCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATGGGTGCATCTAAGGGATACGGGTTTTTa                  <  1:1099159/100‑1 (MQ=255)

                                                                                  aTTAAGCCAACCTTAAGTTGTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                  <  1:3696255/100‑1 (MQ=255)

                                                                                  aTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                  <  2:2023993/100‑1 (MQ=255)

                                                                                  aTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                  <  2:1047243/100‑1 (MQ=255)

                                                                                  aTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                  <  1:2572805/100‑1 (MQ=255)

                                                                                   ttAAGCCAACCTTAAGTTCTTAAGTTTGTCTT‑‑TTATT‑TTTGCTGGATTAAGGAGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                  <  1:2746585/99‑1 (MQ=255)

                                                                                   ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTGTG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                 <  2:2783865/100‑1 (MQ=255)

                                                                                   ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                 <  2:2843893/100‑1 (MQ=255)

                                                                                     aaGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGTGGAATTTATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                 >  1:2771015/1‑98 (MQ=255)

                                                                                     aaGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGt               <  2:1839405/100‑1 (MQ=255)

                                                                                      aGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTa              <  1:988117/100‑1 (MQ=255)

                                                                                      aGCCAACCTTAAGAACTTAAGGTTGGCTTAATTTTGCTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                 >  2:2798455/1‑100 (MQ=21)

                                                                                       gCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAg             >  1:2947033/1‑100 (MQ=255)

                                                                                       gCCAACCTTAAGAACTTAAGGTTGGCTTAATTTTGCTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATg                >  2:507642/1‑100 (MQ=21)

                                                                                       gCCAACCTTAAGAACTTAAGGTTGGC‑T‑‑TAAT‑‑TTTGCTGGATTAAGGGGAATTCATCTTTGCTTATTGGTGCATCTAAGGGAAACGGGTTTTTATGTAGcc           >  2:2449387/1‑100 (MQ=21)

                                                                                        ccAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAAGCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGc            <  1:276430/100‑1 (MQ=255)

                                                                                        ccAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGGTTATTGGTGCATCTAAGGGATACGGGttttttttt               >  1:3570940/1‑93 (MQ=38)

                                                                                        ccAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGc            >  2:3797057/1‑100 (MQ=255)

                                                                                        ccAACCTTAAGAACTTAAGGTTGGC‑T‑‑TAAT‑‑TTTGCTGGATTAAGGTCAATTAATCTTTGCGTATTGGTGCATCTAAGGGATACGGGGTTTTa                  >  1:3612990/1‑92 (MQ=18)

                                                                                         cAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGcc           >  1:3866018/1‑100 (MQ=255)

                                                                                         cAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGTGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGCTTTTATGTAGc            >  1:662317/1‑99 (MQ=38)

                                                                                         cAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGcc           >  1:2005454/1‑100 (MQ=255)

                                                                                         cAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGcc           >  1:3378807/1‑100 (MQ=255)

                                                                                           aCCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCgg         >  1:843698/1‑100 (MQ=255)

                                                                                             cTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGtt       >  1:1348618/1‑100 (MQ=255)

                                                                                              ttAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg      <  1:3764092/100‑1 (MQ=255)

                                                                                              ttAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg      <  2:3858173/100‑1 (MQ=255)

                                                                                              ttAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg      <  1:1781719/100‑1 (MQ=255)

                                                                                              ttAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg      >  1:2759747/1‑100 (MQ=255)

                                                                                                 aGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa   <  2:3574871/100‑1 (MQ=255)

                                                                                                  gTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAt  >  1:3063110/1‑100 (MQ=255)

ATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAT  >  NC_000913/4335555‑4335752

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: `‑`