BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,729,644	T→C	M845T (ATG→ACG)	lhr →	Lhr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,729,644	0	T	C	99.1%	709.6 / ‑3.5	226	M845T (ATG→ACG)	lhr	Lhr
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (118/106); total (120/106)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 234 total aligned reads displayed.

TCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATG  >  NC_000913/1729555‑1729733

tCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGAcgcg                                                                                        >  1:4192163/1‑93 (MQ=255)

tCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGAcgcg                                                                                        >  2:5476296/1‑93 (MQ=255)

 cAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGt                                                                                       <  1:4606644/93‑1 (MQ=255)

  aGGGATAGCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTc                                                                                      <  2:3335277/93‑1 (MQ=255)

   gggATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCg                                                                                     >  1:6784430/1‑93 (MQ=255)

   gggATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCg                                                                                     >  1:1809593/1‑93 (MQ=255)

    ggATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGc                                                                                    >  2:2920213/1‑93 (MQ=255)

       tatCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTg                                                                                 <  1:2703616/93‑1 (MQ=255)

        atCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGa                                                                                >  1:4860484/1‑93 (MQ=255)

         tCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGaa                                                                               <  1:3828567/93‑1 (MQ=255)

          ccGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAAc                                                                              <  1:6225017/93‑1 (MQ=255)

           cGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGTTGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACg                                                                             <  2:3322320/93‑1 (MQ=255)

           cGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACg                                                                             <  2:5571734/93‑1 (MQ=255)

            ggATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACgg                                                                            <  1:2387881/93‑1 (MQ=255)

            ggATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACgg                                                                            >  2:6143541/1‑93 (MQ=255)

             gATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGt                                                                           <  1:2526177/93‑1 (MQ=255)

              aTTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTg                                                                          <  1:3443776/93‑1 (MQ=255)

                tttCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGc                                                                        <  2:4307606/93‑1 (MQ=255)

                tttCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGc                                                                        <  2:1151486/93‑1 (MQ=255)

                    cGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGaaa                                                                    <  2:2994929/93‑1 (MQ=255)

                    cGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGaaa                                                                    <  1:426741/93‑1 (MQ=255)

                     ggTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAAt                                                                   <  2:6571173/93‑1 (MQ=255)

                      gTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAAtt                                                                  >  1:3644899/1‑93 (MQ=255)

                      gTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAAtt                                                                  <  2:7473680/93‑1 (MQ=255)

                             tcGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAatat                                                           >  1:7846538/1‑93 (MQ=255)

                               gAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATcc                                                         >  2:2835568/1‑93 (MQ=255)

                                 aaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGa                                                       >  2:1793403/1‑93 (MQ=255)

                                 aaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGa                                                       >  1:5585459/1‑93 (MQ=255)

                                 aaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGa                                                       <  2:2296268/93‑1 (MQ=255)

                                  aCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGAt                                                      <  2:4657370/93‑1 (MQ=255)

                                   cccTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATg                                                     >  2:734703/1‑93 (MQ=255)

                                    ccTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGt                                                    <  2:5034779/93‑1 (MQ=255)

                                    ccTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGt                                                    >  1:6065660/1‑93 (MQ=255)

                                    ccTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGt                                                    >  2:5768810/1‑93 (MQ=255)

                                         gcgAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGa                                               >  2:7709716/1‑93 (MQ=255)

                                         gcgAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGa                                               <  1:3221321/93‑1 (MQ=255)

                                            aaTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCa                                            <  1:4706003/93‑1 (MQ=255)

                                            aaTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCa                                            >  1:4029456/1‑93 (MQ=255)

                                             aTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCAc                                           >  2:6823807/1‑93 (MQ=255)

                                             aTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCAc                                           >  1:7677289/1‑93 (MQ=255)

                                             aTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCAc                                           >  1:1329323/1‑93 (MQ=255)

                                             aTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCAc                                           <  1:6534862/93‑1 (MQ=255)

                                              tGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACt                                          >  2:1888172/1‑93 (MQ=255)

                                              tGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACt                                          >  1:559456/1‑93 (MQ=255)

                                               gTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTa                                         <  1:1112562/93‑1 (MQ=255)

                                               gTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTa                                         >  1:3259178/1‑93 (MQ=255)

                                               gTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTa                                         >  1:1697830/1‑93 (MQ=255)

                                               gTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTa                                         <  1:7511504/93‑1 (MQ=255)

                                                tCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTAc                                        >  2:3490150/1‑93 (MQ=255)

                                                  tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGc                                      >  1:1233672/1‑93 (MQ=255)

                                                  tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGc                                      >  2:7403941/1‑93 (MQ=255)

                                                   gCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGcc                                     <  1:6012487/93‑1 (MQ=255)

                                                   gCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGcc                                     <  1:200424/93‑1 (MQ=255)

                                                   gCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGcc                                     <  1:7076336/93‑1 (MQ=255)

                                                    cAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGccc                                    >  1:1410485/1‑93 (MQ=255)

                                                      aGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTc                                  >  2:4484959/1‑93 (MQ=255)

                                                       gATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCg                                 >  2:1101296/1‑93 (MQ=255)

                                                        aTGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGc                                >  1:3673178/1‑93 (MQ=255)

                                                        aTGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATGCAAATATCCGATGTAACGACCACTACGCCCTCGc                                >  1:4834220/1‑93 (MQ=255)

                                                         tGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGcc                               <  2:6951354/93‑1 (MQ=255)

                                                         tGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGcc                               <  2:1466555/93‑1 (MQ=255)

                                                         tGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGcc                               >  1:2785659/1‑93 (MQ=255)

                                                           tttATGATCTCCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCtt                             >  1:7032618/1‑93 (MQ=255)

                                                            ttATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCttt                            <  1:2211697/93‑1 (MQ=255)

                                                            ttATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCttt                            <  1:3555/93‑1 (MQ=255)

                                                            ttATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCttt                            <  2:6050862/93‑1 (MQ=255)

                                                             tATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAACTATCCGATGTAACGACCACTACGCCCTCGCCTTTc                           <  2:7831463/93‑1 (MQ=255)

                                                                 aTCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCca                       >  2:3168946/1‑93 (MQ=255)

                                                                  tCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCcac                      >  2:7209856/1‑93 (MQ=255)

                                                                   cTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCcaca                     >  1:6313974/1‑93 (MQ=255)

                                                                    ttCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACaa                    >  2:797729/1‑93 (MQ=255)

                                                                    ttCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACaa                    <  1:5306216/93‑1 (MQ=255)

                                                                    ttCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACaa                    <  1:6925472/93‑1 (MQ=255)

                                                                    ctCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACaa                    <  1:1009605/92‑1 (MQ=255)

                                                                    cgcccGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACaa                    <  1:7278434/91‑1 (MQ=255)

                                                                       ccGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGtt                 >  2:3977831/1‑93 (MQ=255)

                                                                       ccGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGtt                 <  2:4192162/93‑1 (MQ=255)

                                                                       ccGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGtt                 <  2:1757513/93‑1 (MQ=255)

                                                                        cGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGttt                <  1:5164173/93‑1 (MQ=255)

                                                                          cACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTAt              <  2:5632155/93‑1 (MQ=255)

                                                                          cACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTAt              >  1:4567485/1‑93 (MQ=255)

                                                                           aCTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTAtt             <  2:1052031/93‑1 (MQ=255)

                                                                           aCTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTAtt             <  2:2081121/93‑1 (MQ=255)

                                                                            cTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTg            >  2:2000428/1‑93 (MQ=255)

                                                                             tGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCATCGCCTTTCGCCACAAGTTTATTGt           >  1:7749329/1‑93 (MQ=255)

                                                                              ggAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGtt          <  1:842834/93‑1 (MQ=255)

                                                                              ggAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGtt          <  1:5817099/93‑1 (MQ=255)

                                                                              ggAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGta          >  2:2618293/1‑92 (MQ=255)

                                                                               gAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTc         <  2:2438154/93‑1 (MQ=255)

                                                                               gAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTc         <  1:2286457/93‑1 (MQ=255)

                                                                               gAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTc         >  2:1281018/1‑93 (MQ=255)

                                                                                aaCGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCg        >  2:1937434/1‑93 (MQ=255)

                                                                                aaCGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCg        >  1:6945331/1‑93 (MQ=255)

                                                                                aaCGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCg        >  2:1387612/1‑93 (MQ=255)

                                                                                 aCGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCgg       >  1:1337556/1‑93 (MQ=255)

                                                                                 aCGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCgg       >  2:4827308/1‑93 (MQ=255)

                                                                                  cGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGc      <  1:4362277/93‑1 (MQ=255)

                                                                                     ttGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTAt   <  2:7306768/93‑1 (MQ=255)

                                                                                     ttGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTAt   <  1:3542959/93‑1 (MQ=255)

                                                                                      tGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATg  <  2:3846243/93‑1 (MQ=255)

TCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATG  >  NC_000913/1729555‑1729733

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