BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	603,393	0	T	A	20.0%	356.8 / 89.4	160	intergenic (+58/+23)	pheP/ybdG	phenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand): ref base T (57/71); new base A (18/14); total (75/85)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 163 total aligned reads displayed.

GTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTC  >  NC_000913/603294‑603490

gTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAt                                                                                                   >  1:1218097/1‑100 (MQ=255)

gTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAt                                                                                                   <  1:2786889/100‑1 (MQ=255)

  aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTc                                                                                                 <  2:1665500/100‑1 (MQ=255)

   ttCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg                                                                                                <  1:2338310/100‑1 (MQ=255)

    tCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa                                                                                               >  2:3964519/1‑100 (MQ=255)

        gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                           <  1:1446739/100‑1 (MQ=255)

         tttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgaga                                                                                          <  1:731598/100‑1 (MQ=255)

         tttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgaga                                                                                          <  1:1679459/100‑1 (MQ=255)

          ttATGGGATGTAAAACGCTGCGTCGGAAATAAGGGATTCACGCTACACCCGTCAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGagag                                                                                         <  2:225542/100‑1 (MQ=255)

           tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagagn                                                                                        <  1:328086/100‑2 (MQ=39)

           tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                        <  2:1922673/100‑1 (MQ=39)

           tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                        <  1:2089718/100‑1 (MQ=39)

           tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                        <  1:2225583/100‑1 (MQ=39)

           tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGc                                                                                        >  2:245674/1‑100 (MQ=255)

            aTGGTATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                       <  2:1258638/100‑1 (MQ=38)

            aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                       >  2:508181/1‑100 (MQ=255)

             tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCgg                                                                                      <  1:1530176/100‑1 (MQ=38)

             tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                      >  1:1338270/1‑100 (MQ=255)

             gatcttttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCgg                                                                                      <  2:3865589/95‑1 (MQ=38)

              tctaTTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGt                                                                                     <  1:3787416/97‑1 (MQ=38)

              ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                     >  1:1209071/1‑100 (MQ=255)

               gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAATCGCTCTCATCCATTCGATGAGAGCGGtt                                                                                    >  1:3704005/1‑100 (MQ=255)

                 aTTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGtaga                                                                                  >  2:3707287/1‑97 (MQ=38)

                 aTTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttt                                                                                  <  1:295688/100‑1 (MQ=255)

                 aTTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttt                                                                                  >  1:3871591/1‑100 (MQ=255)

                  tttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttt                                                                                 <  1:3973434/100‑1 (MQ=255)

                  tttAAAACGCTGCGGCGGAAATAAGGCATTCACGCTACATCCGAAAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttt                                                                                 <  2:389640/100‑1 (MQ=255)

                    tAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                               <  2:59457/100‑1 (MQ=255)

                     aaaaCGCGGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                              <  1:3503829/100‑1 (MQ=255)

                      aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGAGAAAACGATGTCAACCATCCGAAACCGCTCTCGTCCATTCGATGAGTGCGGTTTTTTTaa                                                                             <  1:3597022/100‑1 (MQ=255)

                      aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTaa                                                                             <  1:1861022/100‑1 (MQ=255)

                      aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTaa                                                                             >  2:1057595/1‑100 (MQ=255)

                      aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTaa                                                                             <  2:237760/100‑1 (MQ=255)

                        aCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAtt                                                                           <  2:3116477/100‑1 (MQ=255)

                        aCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAtt                                                                           <  2:2364395/100‑1 (MQ=255)

                             gCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGc                                                                      <  2:2052111/100‑1 (MQ=255)

                             gCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGc                                                                      >  2:118300/1‑100 (MQ=255)

                              cGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCt                                                                     <  1:3458200/100‑1 (MQ=255)

                               gTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCtt                                                                    >  1:1956595/1‑100 (MQ=255)

                               gTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCtt                                                                    <  2:1381154/100‑1 (MQ=255)

                                 cGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaa                                                                  >  2:2328047/1‑100 (MQ=255)

                                  ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                 >  1:86976/1‑100 (MQ=255)

                                  ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                 >  2:3156505/1‑100 (MQ=255)

                                  ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                 <  2:2676815/100‑1 (MQ=255)

                                    aaaTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                               <  1:2882527/100‑1 (MQ=255)

                                      aTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCa                                                             <  2:1468627/100‑1 (MQ=255)

                                        aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAcc                                                           <  2:382856/100‑1 (MQ=255)

                                        aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAcc                                                           <  2:2775920/100‑1 (MQ=255)

                                          ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                         >  1:1149681/1‑100 (MQ=255)

                                            cATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGcc                                                       >  1:2761881/1‑100 (MQ=255)

                                             aTTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCa                                                      >  2:858932/1‑100 (MQ=255)

                                             aTTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCa                                                      >  2:1468465/1‑100 (MQ=255)

                                                cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                   >  2:3551195/1‑100 (MQ=255)

                                                  cGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagagag                                                 >  1:868487/1‑100 (MQ=255)

                                                   gCTACATCCGACAAAACGATGTCAACCCTCCGAAACCGCTCTCATCCATGCTATGAGGGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGc                                                <  2:2662658/100‑1 (MQ=255)

                                                   gCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGGGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGc                                                <  1:3523973/100‑1 (MQ=255)

                                                   gCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGc                                                >  2:23105/1‑100 (MQ=255)

                                                      aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGaa                                             >  2:2752518/1‑100 (MQ=255)

                                                              cAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg                                     <  1:1971639/100‑1 (MQ=255)

                                                                aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                   >  1:3735030/1‑100 (MQ=255)

                                                                  aCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCgg                                 <  1:3857710/100‑1 (MQ=255)

                                                                   cGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGt                                <  1:96722/100‑1 (MQ=255)

                                                                     aTGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTCCTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCAGGTTg                              <  1:126723/100‑1 (MQ=255)

                                                                     aTGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTg                              <  1:3825826/100‑1 (MQ=255)

                                                                     aTGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTg                              >  1:2580591/1‑100 (MQ=255)

                                                                       gTCAACCATCCGAAACCGCGCTCATCCATTCGATGATAGCGTTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTggn                            <  2:3203948/100‑2 (MQ=255)

                                                                        tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                           <  1:3753144/100‑1 (MQ=255)

                                                                         cAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGa                          >  1:2690588/1‑100 (MQ=255)

                                                                             cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                      <  2:1535830/100‑1 (MQ=255)

                                                                              aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                     <  1:233137/100‑1 (MQ=255)

                                                                               tCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                    <  2:1732936/100‑1 (MQ=255)

                                                                                ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTTCTTAAATGCACCCGCCAGAAAGCGAATATCATTGCCGGTTGGCGACtgatga                   >  2:3580739/1‑100 (MQ=255)

                                                                                ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatga                   <  2:3665303/100‑1 (MQ=255)

                                                                                  gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAg                 <  1:453251/100‑1 (MQ=255)

                                                                                  aaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAg                 >  1:2898023/2‑100 (MQ=38)

                                                                                   aaaCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGt                <  2:278484/100‑1 (MQ=255)

                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGCGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc             >  2:931541/1‑100 (MQ=38)

                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc             >  1:2075782/1‑100 (MQ=38)

                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc             >  1:3394463/1‑100 (MQ=38)

                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc             >  2:3193514/1‑100 (MQ=38)

                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc             >  1:918661/1‑100 (MQ=38)

                                                                                      ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc             >  1:3984308/1‑100 (MQ=255)

                                                                                       cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa            >  2:3276972/1‑100 (MQ=38)

                                                                                       cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACACGCCAGAGAGGGAATATCATTGTCGGTTGGCGACTGATGAAGTCGcc            >  2:837541/1‑99 (MQ=255)

                                                                                        gCTCTCGTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGGAAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           <  2:3163223/100‑1 (MQ=255)

                                                                                        gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           >  2:2550564/1‑100 (MQ=39)

                                                                                        gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           >  1:3649046/1‑100 (MQ=39)

                                                                                        gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           >  2:177028/1‑100 (MQ=39)

                                                                                         ctctCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa          >  1:2538939/1‑100 (MQ=255)

                                                                                          tctcATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGAc         >  1:3016497/1‑100 (MQ=39)

                                                                                          tctcATCCATTTGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGGGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGAc         <  2:204015/100‑1 (MQ=255)

                                                                                           ctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGAcc        >  1:1021415/1‑100 (MQ=255)

                                                                                            tcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCa       <  2:2602797/100‑1 (MQ=255)

                                                                                             cATCCATTCGATGAGAGCGGTTTGTTTAATTACTGCTTAAATTCACCCGCCAGATAGCGAATATCATTGTCGGTTGGCGACTGACGAAGTCGCAGACCaa      <  1:3548381/100‑1 (MQ=255)

                                                                                              aTCCATTCGATGAGAGCGGTTTTTTTTATTACTTCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaaa     >  2:2380108/1‑100 (MQ=255)

                                                                                               tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc    <  1:601581/100‑1 (MQ=255)

                                                                                               tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc    <  2:1930844/100‑1 (MQ=255)

                                                                                               tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc    <  1:1660430/100‑1 (MQ=255)

                                                                                                 naTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActc  >  1:3234809/2‑100 (MQ=255)

                                                                                                 cATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActc  >  2:707143/1‑100 (MQ=255)

GTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTC  >  NC_000913/603294‑603490

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