BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,973,386	A→G	S298S (AGT→AGC)	motB ←	MotB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,973,386	0	A	G	86.6%	790.4 / 102.4	321	S298S (AGT→AGC)	motB	MotB
Reads supporting (aligned to +/- strand): ref base A (25/18); new base G (147/131); total (172/149)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 321 total aligned reads displayed.

TCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGG  >  NC_000913/1973296‑1973474

  gTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGa                                                                                      <  1:5267926/93‑1 (MQ=255)

  gTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGa                                                                                      <  2:1413405/93‑1 (MQ=255)

  gTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTgg                                                                                      >  1:5738635/1‑92 (MQ=255)

   tCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGAc                                                                                     >  2:820416/1‑93 (MQ=255)

   tCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGAc                                                                                     >  1:2436472/1‑93 (MQ=255)

    ccGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGAcc                                                                                    <  1:4425676/93‑1 (MQ=255)

     cGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCt                                                                                   >  2:5926767/1‑93 (MQ=255)

      gCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCtg                                                                                  <  1:5046932/93‑1 (MQ=255)

      gCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCtg                                                                                  >  1:2935036/1‑93 (MQ=255)

       cTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCtgt                                                                                 >  2:2504880/1‑93 (MQ=255)

       cTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCtgt                                                                                 >  1:3362416/1‑93 (MQ=255)

        ttCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCtgtg                                                                                >  2:6679485/1‑93 (MQ=255)

        ttCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCtgtg                                                                                >  2:6928352/1‑93 (MQ=255)

        ttCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCtgtg                                                                                >  2:3886600/1‑93 (MQ=255)

         tcatcaAAAAATGTCTGATAAAAATCGCTTATATACATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTgtgg                                                                               <  2:3751417/93‑1 (MQ=255)

          catcaAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGAGAACGCTGACCTgtggt                                                                              <  2:3208625/93‑1 (MQ=255)

           atcaAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTgtggtg                                                                             >  2:1672851/1‑93 (MQ=255)

            tcaAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGc                                                                            >  2:3680620/1‑93 (MQ=255)

                aaaaTGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAAcc                                                                        >  1:1318367/1‑93 (MQ=255)

                   aTGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCa                                                                     <  2:4741359/93‑1 (MQ=255)

                   aTGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCa                                                                     >  1:5460686/1‑93 (MQ=255)

                   aTGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCa                                                                     <  2:3335854/93‑1 (MQ=255)

                    tGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAg                                                                    >  2:4117149/1‑93 (MQ=255)

                     gTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAgg                                                                   >  1:7095808/1‑93 (MQ=255)

                     gTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAgg                                                                   <  1:1680571/93‑1 (MQ=255)

                     gTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCann                                                                   >  2:2325986/1‑91 (MQ=255)

                     gTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAgg                                                                   >  1:4881250/1‑93 (MQ=255)

                       cTGATAAAAATCGCTTATATCCATGCTCACGCTGTCGCCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGtt                                                                 <  1:5001232/93‑1 (MQ=255)

                        tGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGttt                                                                <  2:4572865/93‑1 (MQ=255)

                        tGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGttt                                                                >  2:4959274/1‑93 (MQ=255)

                        tGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGttt                                                                <  2:2103896/93‑1 (MQ=255)

                        tGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGttt                                                                >  1:3527564/1‑93 (MQ=255)

                         gATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGtttt                                                               <  2:2701100/93‑1 (MQ=255)

                         gATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGtttt                                                               >  1:1810246/1‑93 (MQ=255)

                         gATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGtttt                                                               >  2:4784988/1‑93 (MQ=255)

                          aTAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGttttt                                                              >  2:6262318/1‑93 (MQ=255)

                          aTAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGttttt                                                              <  1:2164617/93‑1 (MQ=255)

                           tAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTc                                                             <  1:7236763/93‑1 (MQ=255)

                           tAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTc                                                             >  1:3726813/1‑93 (MQ=255)

                           tAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTc                                                             >  2:3119769/1‑93 (MQ=255)

                            aaaaaTCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTcc                                                            <  1:3341848/93‑1 (MQ=255)

                            aaaaaTCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTcc                                                            <  2:7352400/93‑1 (MQ=255)

                             aaaaTCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGTACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCa                                                           <  1:968685/93‑1 (MQ=255)

                             aaaaTCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCa                                                           <  1:1751929/93‑1 (MQ=255)

                             aaaaTCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCa                                                           >  2:3543353/1‑93 (MQ=255)

                               aaTCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAgg                                                         >  2:1820277/1‑93 (MQ=255)

                                 tCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGgc                                                       >  1:5610432/1‑93 (MQ=255)

                                 tCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGgc                                                       <  2:3979049/93‑1 (MQ=255)

                                   gCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGgcgc                                                     <  1:2344065/93‑1 (MQ=255)

                                   gCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGgcgc                                                     >  2:3271806/1‑93 (MQ=255)

                                   gCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGgcgc                                                     <  1:4284679/93‑1 (MQ=255)

                                   gCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGCTTTTTCCAGGgcgc                                                     <  2:6435522/93‑1 (MQ=255)

                                    cTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCt                                                    >  1:6791132/1‑93 (MQ=255)

                                    cTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCt                                                    >  2:4192176/1‑93 (MQ=255)

                                     ttATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCtt                                                   >  1:3005508/1‑93 (MQ=255)

                                     ttATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCtt                                                   >  1:2490920/1‑93 (MQ=255)

                                     ttATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCtt                                                   >  1:6805784/1‑93 (MQ=255)

                                       atatCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTAc                                                 <  2:5620711/93‑1 (MQ=255)

                                        tatCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACt                                                >  2:1547981/1‑93 (MQ=255)

                                        tatCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACt                                                >  2:3011392/1‑93 (MQ=255)

                                            cATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCt                                            >  2:5060483/1‑93 (MQ=255)

                                             aTGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTc                                           <  1:1122570/93‑1 (MQ=255)

                                              tGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCa                                          >  1:1241506/1‑93 (MQ=255)

                                              tGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCa                                          >  1:6029900/1‑93 (MQ=255)

                                               gCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCAt                                         <  1:1555146/93‑1 (MQ=255)

                                                 tCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTc                                       <  2:3495457/93‑1 (MQ=255)

                                                 tCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTc                                       <  2:6350714/93‑1 (MQ=255)

                                                   aCGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTg                                     >  1:3150180/1‑93 (MQ=255)

                                                   aCGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTg                                     <  2:4462140/93‑1 (MQ=255)

                                                    cGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTgg                                    >  1:5122626/1‑93 (MQ=255)

                                                      cTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCt                                  >  1:6210029/1‑93 (MQ=255)

                                                      cTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCt                                  <  1:4205032/93‑1 (MQ=255)

                                                        gTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCttt                                >  1:1050458/1‑93 (MQ=255)

                                                        gTCACCTCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCttt                                >  2:3429640/1‑93 (MQ=255)

                                                        gTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCttt                                >  1:4809446/1‑93 (MQ=255)

                                                              tCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGt                          <  1:6621979/93‑1 (MQ=255)

                                                              tCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGt                          >  2:5560566/1‑93 (MQ=255)

                                                              tCGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGt                          >  2:5085097/1‑93 (MQ=255)

                                                               cGGTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGtt                         <  2:7142374/93‑1 (MQ=255)

                                                                 gTTCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGtttt                       <  1:6449510/93‑1 (MQ=255)

                                                                  ttCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTc                      <  1:3865506/93‑1 (MQ=255)

                                                                  ttCGGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTc                      >  1:7486170/1‑93 (MQ=255)

                                                                  ttCAGCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTc                      <  2:4642509/93‑1 (MQ=255)

                                                                      gCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTTATTCTGGCTTTCGGCGTTTTCATGc                  >  1:1962073/1‑93 (MQ=255)

                                                                      gCTGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGc                  <  1:1430632/93‑1 (MQ=255)

                                                                        tGATGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCaa                >  2:3166956/1‑93 (MQ=255)

                                                                           tGGCATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAAt             >  2:6767148/1‑93 (MQ=255)

                                                                             gCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATgg           <  1:6941459/93‑1 (MQ=255)

                                                                              cATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGc          <  2:5721433/93‑1 (MQ=255)

                                                                              cATTGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGc          <  2:1319187/93‑1 (MQ=255)

                                                                              cATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATACTGGCTTTCGGCGTTTTCATGCAAAATGGc          >  1:2547475/1‑93 (MQ=255)

                                                                               aTTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGcc         >  1:2246658/1‑93 (MQ=255)

                                                                                ttGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCt        >  2:28965/1‑93 (MQ=255)

                                                                                ttGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCt        <  2:5254285/93‑1 (MQ=255)

                                                                                ttGTGGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCt        <  1:3411773/93‑1 (MQ=255)

                                                                                   tgGGAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGtt     <  2:6073734/93‑1 (MQ=255)

                                                                                    gggAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTc    <  1:6287536/93‑1 (MQ=255)

                                                                                     ggAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCg   >  2:1986511/1‑93 (MQ=255)

                                                                                     ggAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCg   >  1:4521555/1‑93 (MQ=255)

                                                                                      gAACGCTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCgg  <  2:4337576/93‑1 (MQ=255)

TCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGG  >  NC_000913/1973296‑1973474

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