BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,456,531	A→G	M82V (ATG→GTG)	b1394 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,456,531	0	A	G	98.0%	310.8 / ‑4.4	101	M82V (ATG→GTG)	b1394	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (49/50); total (49/52)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 105 total aligned reads displayed.

TTACTTACAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCA  >  NC_000913/1456441‑1456620

ttacttacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTg                                                                                         >  2:581537/1‑93 (MQ=255)

   cttacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCa                                                                                      <  1:244293/93‑1 (MQ=255)

    ttacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAg                                                                                     <  2:3070896/93‑1 (MQ=255)

    ttacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAg                                                                                     <  2:1827459/93‑1 (MQ=255)

     tacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGt                                                                                    <  1:2347580/93‑1 (MQ=255)

      acAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGtt                                                                                   >  1:3482586/1‑93 (MQ=255)

         ggTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGaa                                                                                <  1:4684491/93‑1 (MQ=255)

         ggTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGaa                                                                                >  2:4401542/1‑93 (MQ=255)

          gTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAAc                                                                               <  1:1326866/93‑1 (MQ=255)

          gTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAAc                                                                               <  2:4059526/93‑1 (MQ=255)

               gggCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCAAAGGCACCGCACCGGATTTAGGGGGGTCAGTTGAACGTTTc                                                                          >  1:38629/1‑93 (MQ=255)

                ggCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCt                                                                         >  1:323021/1‑93 (MQ=255)

                 gcgcGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCta                                                                        >  1:3061810/1‑93 (MQ=255)

                  cgcgGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCtat                                                                       <  2:16049/93‑1 (MQ=255)

                  cgcgGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCtat                                                                       >  2:3265410/1‑93 (MQ=255)

                     ggATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAc                                                                    <  1:2437846/93‑1 (MQ=255)

                     ggATTTTGTGCTGGGCAGGATATTAACGATCGTAACGTTGATCCCAAAGGCTCCGCACCGGATTTAGGGGTGTCAGTTGCACGTTTCTATAAc                                                                    <  1:2437847/93‑1 (MQ=255)

                      gATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAcc                                                                   >  1:254235/1‑93 (MQ=255)

                      gATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAcc                                                                   >  2:4329387/1‑93 (MQ=255)

                       aTTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAccc                                                                  >  2:74769/1‑93 (MQ=255)

                       aTTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAccc                                                                  >  1:1792549/1‑93 (MQ=255)

                        ttttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCa                                                                 <  1:2523792/93‑1 (MQ=255)

                        ttttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCa                                                                 >  1:4745511/1‑93 (MQ=255)

                        ttttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCa                                                                 <  1:1866401/93‑1 (MQ=255)

                        ttttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCa                                                                 >  1:4245786/1‑93 (MQ=255)

                         tttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCAc                                                                <  2:1063702/93‑1 (MQ=255)

                           tgtgCTGGTGAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTg                                                              <  2:1249734/93‑1 (MQ=255)

                            gtgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCGCTgg                                                             <  2:3029401/93‑1 (MQ=255)

                             tgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGt                                                            >  1:743115/1‑93 (MQ=255)

                              gCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTa                                                           >  2:1746558/1‑93 (MQ=255)

                              gCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTa                                                           <  1:3023195/93‑1 (MQ=255)

                               cTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTAc                                                          >  2:1333645/1‑93 (MQ=255)

                               cTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTAc                                                          <  1:3095451/93‑1 (MQ=255)

                                tGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACg                                                         >  2:3442834/1‑93 (MQ=255)

                                tGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACg                                                         >  2:1165169/1‑93 (MQ=255)

                                 ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGt                                                        <  2:1463034/93‑1 (MQ=255)

                                     aGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGcc                                                    >  2:2803572/1‑93 (MQ=255)

                                      ggATCTTAACGATCGTACCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCt                                                   <  1:4217013/93‑1 (MQ=255)

                                      ggATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCt                                                   <  2:3584746/93‑1 (MQ=255)

                                        aTCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTgg                                                 >  2:3509612/1‑93 (MQ=255)

                                        aTCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTgg                                                 >  2:905938/1‑93 (MQ=255)

                                         tCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGc                                                >  2:3457858/1‑93 (MQ=255)

                                          cTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCa                                               >  2:4642298/1‑93 (MQ=255)

                                           ttAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaa                                              >  2:841379/1‑93 (MQ=255)

                                             aaCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaaaa                                            >  2:2228544/1‑93 (MQ=255)

                                             aaCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaaaa                                            >  2:4068577/1‑93 (MQ=255)

                                              acaaTCGTCGCGTTGAGCACACAGGCACCGCAACGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTAGCATGGCAAAAc                                           <  2:2772728/90‑1 (MQ=255)

                                              aCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAAc                                           <  1:934710/93‑1 (MQ=255)

                                               cGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACt                                          <  1:3805254/93‑1 (MQ=255)

                                                gATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTCGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTg                                         <  2:2055893/93‑1 (MQ=255)

                                                 aTCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGc                                        >  2:3704982/1‑93 (MQ=255)

                                                 aTCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGc                                        <  2:3463627/93‑1 (MQ=255)

                                                 aTCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGATGAACGATTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGc                                        >  2:3704981/1‑93 (MQ=255)

                                                 aTCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGc                                        <  2:4326549/93‑1 (MQ=255)

                                                  tCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGcc                                       >  2:4189311/1‑93 (MQ=255)

                                                  tCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGcc                                       <  2:2842490/93‑1 (MQ=255)

                                                  tCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGcc                                       >  2:2285415/1‑93 (MQ=255)

                                                  tCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGcc                                       >  2:4229691/1‑93 (MQ=255)

                                                   cGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCa                                      <  1:2700503/93‑1 (MQ=255)

                                                     tAACGTTGATCCCACAGGCCCCACACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCaaa                                    >  1:4727818/1‑93 (MQ=255)

                                                       aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCACAAc                                  >  1:4289619/1‑93 (MQ=255)

                                                       aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAc                                  <  1:4544431/93‑1 (MQ=255)

                                                       aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAc                                  <  1:86089/93‑1 (MQ=255)

                                                        cGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAcc                                 <  1:4174658/93‑1 (MQ=255)

                                                         gTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCg                                <  1:3377888/93‑1 (MQ=255)

                                                          ttGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                               <  1:4286367/93‑1 (MQ=255)

                                                            gATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTg                             >  2:1927034/1‑93 (MQ=255)

                                                             aTCCCACAGGCCCCGCACCGGATCTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGa                            <  2:2619241/93‑1 (MQ=255)

                                                               cccACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATc                          <  2:2706148/93‑1 (MQ=255)

                                                               cccACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATc                          >  1:607997/1‑93 (MQ=255)

                                                                ccACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTATCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGAGAtat                         <  1:4221719/93‑3 (MQ=255)

                                                                  acaGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgt                       <  2:3269501/93‑1 (MQ=255)

                                                                   caGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgtg                      <  2:4103470/93‑1 (MQ=255)

                                                                   caGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgtg                      <  1:4620999/93‑1 (MQ=255)

                                                                    agaccccGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGc                     <  2:486758/90‑1 (MQ=255)

                                                                    aGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAACACCGGTGATCTGTGc                     >  2:677866/1‑93 (MQ=255)

                                                                    aGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGc                     <  1:4303514/93‑1 (MQ=255)

                                                                     ggCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCa                    <  1:1492994/93‑1 (MQ=255)

                                                                     ggCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCa                    >  2:1255765/1‑93 (MQ=255)

                                                                      gCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAg                   >  1:1918965/1‑93 (MQ=255)

                                                                      gCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAg                   >  1:969055/1‑93 (MQ=255)

                                                                        cccGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTc                 <  1:4146830/93‑1 (MQ=255)

                                                                         ccGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCa                >  1:3390624/1‑93 (MQ=255)

                                                                           gCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAt              >  2:3327163/1‑93 (MQ=255)

                                                                           gCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAt              <  1:956611/93‑1 (MQ=255)

                                                                            cACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtg             >  2:3970216/1‑93 (MQ=255)

                                                                            cACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtg             <  2:267874/93‑1 (MQ=255)

                                                                            cACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtg             <  2:641866/93‑1 (MQ=255)

                                                                            cACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtg             <  2:3278866/93‑1 (MQ=255)

                                                                             aCCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtgg            <  1:4447428/93‑1 (MQ=255)

                                                                             aCCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtgg            >  2:3488022/1‑93 (MQ=255)

                                                                              ccGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggc           <  1:1793880/93‑1 (MQ=255)

                                                                                ggATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgt         <  2:265404/93‑1 (MQ=255)

                                                                                ggATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgt         <  2:3061594/93‑1 (MQ=255)

                                                                                   tttAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgtggc      >  2:1057102/1‑93 (MQ=255)

                                                                                   tttAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgtggc      <  2:2531674/93‑1 (MQ=255)

                                                                                    ttAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgtggcg     >  1:518909/1‑93 (MQ=255)

                                                                                      aGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTggcggc   >  2:1883092/1‑93 (MQ=255)

                                                                                      aGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTggcggc   >  1:4748867/1‑93 (MQ=255)

                                                                                       ggggTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggca  <  1:538549/93‑1 (MQ=255)

TTACTTACAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCA  >  NC_000913/1456441‑1456620

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