BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,462,362	T→C	I267T (ATT→ACT)	b1399 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,462,362	0	T	C	99.6%	770.2 / ‑4.8	235	I267T (ATT→ACT)	b1399	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (110/124); total (111/124)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

Only 100 of 246 total aligned reads displayed.

GCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTA  >  NC_000913/1462272‑1462452

gCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCCCACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACActt                                                                                          <  2:1506804/93‑1 (MQ=255)

gCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACActt                                                                                          >  1:2616930/1‑93 (MQ=255)

   gtgtCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATc                                                                                       <  1:6024951/93‑1 (MQ=255)

   gtgtCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATc                                                                                       <  1:6422848/93‑1 (MQ=255)

      tcgtcCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAgc                                                                                    <  2:4996473/93‑1 (MQ=255)

      tcgtcCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAgc                                                                                    <  1:2350380/93‑1 (MQ=255)

        gtcCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAgcgc                                                                                  >  1:4148354/1‑93 (MQ=255)

        gtcCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAgcgc                                                                                  >  1:4672352/1‑93 (MQ=255)

        gtcCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAgcgc                                                                                  >  2:4230603/1‑93 (MQ=255)

          ccTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGt                                                                                <  1:2029081/93‑1 (MQ=255)

          ccTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGt                                                                                <  2:6707771/93‑1 (MQ=255)

           cTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTa                                                                               >  1:808533/1‑93 (MQ=255)

           cTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTa                                                                               <  2:4863226/93‑1 (MQ=255)

           cTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTa                                                                               >  2:285719/1‑93 (MQ=255)

            ttAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAg                                                                              >  1:731877/1‑93 (MQ=255)

            ttAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGACATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAg                                                                              >  1:6623448/1‑93 (MQ=255)

             tAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAgc                                                                             <  1:5136228/93‑1 (MQ=255)

                aGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGcc                                                                          >  1:6449479/1‑93 (MQ=255)

                 gACCCATTGTTGCCGGAGGGGTTGCTTCCGGCACACTGGGCGAGGCATACGGCGCGTCAGCTGGGCATCAACACTTACAAGCGCGTAGGGccg                                                                         >  2:4704911/1‑92 (MQ=255)

                 gACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCct                                                                         <  2:6163290/93‑1 (MQ=255)

                   cccATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCctgc                                                                       <  2:2908209/93‑1 (MQ=255)

                    ccATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCctgct                                                                      >  2:5545419/1‑93 (MQ=255)

                    ccATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCctgct                                                                      >  2:3448470/1‑93 (MQ=255)

                      aTTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCtgctgc                                                                    >  1:3992614/1‑93 (MQ=255)

                       ttgttgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTgctgct                                                                   >  2:2646716/1‑93 (MQ=255)

                       ttgttgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTgctgct                                                                   >  1:3102512/1‑93 (MQ=255)

                        tgttgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCtt                                                                  <  2:327788/93‑1 (MQ=255)

                        tgttgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCtt                                                                  <  2:4354877/93‑1 (MQ=255)

                          ttgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTa                                                                <  2:3125042/93‑1 (MQ=255)

                          ttgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTa                                                                <  1:5920300/93‑1 (MQ=255)

                          ttgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGCAGCGCCTGCTGCtttg                                                                >  1:5443391/1‑92 (MQ=255)

                           tgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAg                                                               <  1:6104589/93‑1 (MQ=255)

                           tgCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAg                                                               <  2:1129962/93‑1 (MQ=255)

                            gCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGc                                                              <  2:6593499/93‑1 (MQ=255)

                            gCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGc                                                              >  1:2262891/1‑93 (MQ=255)

                            gCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGc                                                              <  2:3597191/93‑1 (MQ=255)

                             ccGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGcg                                                             <  1:5792727/93‑1 (MQ=255)

                               ggaggaGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGcgtt                                                           >  2:1461639/1‑93 (MQ=255)

                               ggaggaGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGcgtt                                                           <  1:5209370/93‑1 (MQ=255)

                               ggaggaGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGcgtt                                                           >  2:248842/1‑93 (MQ=255)

                                gaggagTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGcgttc                                                          <  2:2085101/93‑1 (MQ=255)

                                gaggagTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGcgttc                                                          >  1:1017263/1‑93 (MQ=255)

                                  ggagTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCCGCGCGTAGCGCCTGCTGCTTTAGcgttcgt                                                        >  1:3366601/1‑93 (MQ=255)

                                  ggagTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGcgttcgt                                                        >  2:3646698/1‑93 (MQ=255)

                                  ggagTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGcgttcgt                                                        <  1:2556718/93‑1 (MQ=255)

                                     gTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAg                                                     <  2:475272/93‑1 (MQ=255)

                                       tGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTg                                                   >  1:884222/1‑93 (MQ=255)

                                               gCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGt                                           <  2:5591386/93‑1 (MQ=255)

                                               gCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGt                                           <  2:2370719/93‑1 (MQ=255)

                                                 acacTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGa                                         >  2:40333/1‑93 (MQ=255)

                                                 acacTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGa                                         <  2:5842767/93‑1 (MQ=255)

                                                  cacTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGaa                                        >  2:5210406/1‑93 (MQ=255)

                                                    cTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAAc                                      <  1:2324762/93‑1 (MQ=255)

                                                      gggcagggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCt                                    <  1:6295836/93‑1 (MQ=255)

                                                      gggcagggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCt                                    >  1:3658308/1‑93 (MQ=255)

                                                       ggcagggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtc                                   <  1:1616935/93‑1 (MQ=255)

                                                       ggcagggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtc                                   <  2:4512724/93‑1 (MQ=255)

                                                        gcagggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtcg                                  >  2:4152039/1‑93 (MQ=255)

                                                        gcagggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtcg                                  <  2:3316171/93‑1 (MQ=255)

                                                         cagggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtcgg                                 <  2:1947531/93‑1 (MQ=255)

                                                         cagggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtcgg                                 <  2:6291973/93‑1 (MQ=255)

                                                         cagggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtcgg                                 >  1:3783433/1‑93 (MQ=255)

                                                          agggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtcggt                                <  2:3730720/93‑1 (MQ=255)

                                                          agggcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtcggt                                <  1:1819224/93‑1 (MQ=255)

                                                             gcaTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCtcggtcgg                             <  1:5694784/93‑1 (MQ=255)

                                                              caTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTcggtcggt                            <  1:1119347/93‑1 (MQ=255)

                                                              caTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTcggtcggt                            <  1:3953185/93‑1 (MQ=255)

                                                               aTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTg                           <  2:425956/93‑1 (MQ=255)

                                                               aTACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAAGAGAAAACTCGGTCGGTg                           >  1:3201382/1‑93 (MQ=255)

                                                                tACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGa                          >  2:1210349/1‑93 (MQ=255)

                                                                tACGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGa                          <  1:525493/93‑1 (MQ=255)

                                                                  cGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAAc                        <  2:3381017/93‑1 (MQ=255)

                                                                  cGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAAc                        <  1:2764047/93‑1 (MQ=255)

                                                                  cGGCGCGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAAc                        >  2:1729624/1‑93 (MQ=255)

                                                                  cGGCGATTCAGCTGTGTATCAACTCTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAAc                        <  1:4207393/93‑1 (MQ=255)

                                                                    gcgcGTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAActg                      <  1:1289096/93‑1 (MQ=255)

                                                                      gcgTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAActgcc                    <  1:5398597/93‑1 (MQ=255)

                                                                      gcgTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAActgcc                    <  2:5680915/93‑1 (MQ=255)

                                                                       cgTCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAActgcct                   <  2:2290317/93‑1 (MQ=255)

                                                                         tCAGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAActgcctgc                 <  1:3793603/93‑1 (MQ=255)

                                                                           aGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTgcgc               <  2:52200/93‑1 (MQ=255)

                                                                           aGCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTgcgc               >  1:6150295/1‑93 (MQ=255)

                                                                            gCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGcc              >  1:4389442/1‑93 (MQ=255)

                                                                            gCTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGcc              <  2:838290/93‑1 (MQ=255)

                                                                             cTGTGTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCg             >  1:3159931/1‑93 (MQ=255)

                                                                              tgtgTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCgg            >  1:3574346/1‑93 (MQ=255)

                                                                              tgtgTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCgg            <  1:1463918/93‑1 (MQ=255)

                                                                              tgtgTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCgg            >  2:3806753/1‑93 (MQ=255)

                                                                              tgtgTATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAAACTCGGTCGGTGAACTGCCTGCGCCgg            >  1:6228076/1‑93 (MQ=255)

                                                                               gtgtATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCggg           >  2:3228936/1‑93 (MQ=255)

                                                                                 gtATCAACACTTATTAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGaa         <  2:1619333/93‑1 (MQ=255)

                                                                                 gtATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGaa         <  2:5390570/93‑1 (MQ=255)

                                                                                 gtATCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGCGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGaa         >  1:5907996/1‑93 (MQ=255)

                                                                                   aTCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGc       <  1:5035993/93‑1 (MQ=255)

                                                                                    tCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGcc      <  2:3837498/93‑1 (MQ=255)

                                                                                    tCAACACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGcc      <  1:4862972/93‑1 (MQ=255)

                                                                                      aaCACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTg    >  2:5969427/1‑93 (MQ=255)

                                                                                      aaCACTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGGTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTg    >  2:1119166/1‑93 (MQ=255)

                                                                                      aaCACTTATCAGCGCGCAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTg    >  1:104881/1‑93 (MQ=255)

                                                                                        cATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTa  >  2:2568736/1‑93 (MQ=255)

GCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGCAGGGCATACGGCGCGTCAGCTGTGTATCAACATTTATCAGCGCGTAGCGCCTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTA  >  NC_000913/1462272‑1462452

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