BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,605,226	G→A	T68T (ACG→ACA)	b1518 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,605,226	0	G	A	78.3%	463.0 / 114.8	230	T68T (ACG→ACA)	b1518	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (24/26); new base A (98/82); total (122/108)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-04

Only 100 of 232 total aligned reads displayed.

TCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGG  >  NC_000913/1605135‑1605315

tCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGc                                                                                          <  1:4149772/93‑1 (MQ=255)

  gATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACAccc                                                                                        >  2:160246/1‑93 (MQ=255)

  gATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACAccc                                                                                        >  1:372052/1‑93 (MQ=255)

    tGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCa                                                                                      >  2:394536/1‑93 (MQ=255)

       cTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTa                                                                                   <  2:4607954/93‑1 (MQ=255)

        ttACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTAc                                                                                  >  2:30720/1‑93 (MQ=255)

           cAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACaaa                                                                               >  1:5149967/1‑93 (MQ=255)

             ggACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAAc                                                                             <  2:3711161/93‑1 (MQ=255)

              gACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAAcc                                                                            >  1:928947/1‑93 (MQ=255)

              gACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAAcc                                                                            <  2:7306789/93‑1 (MQ=255)

               aCCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCt                                                                           <  1:4350458/93‑1 (MQ=255)

               aCCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCt                                                                           <  2:6967201/93‑1 (MQ=255)

                cccGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCtg                                                                          <  2:2841683/93‑1 (MQ=255)

                  cGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCtgtg                                                                        >  1:4702314/1‑93 (MQ=255)

                  cGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCtgtg                                                                        <  2:6055001/93‑1 (MQ=255)

                   ggAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTgtgt                                                                       >  1:1151465/1‑93 (MQ=255)

                   ggAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTgtgt                                                                       <  2:1892705/93‑1 (MQ=255)

                   ggAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTgtgt                                                                       >  2:3129116/1‑93 (MQ=255)

                    gAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTc                                                                      <  2:2332990/93‑1 (MQ=255)

                    gAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTc                                                                      >  1:228139/1‑93 (MQ=255)

                    gAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTCGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTc                                                                      >  2:3156627/1‑93 (MQ=255)

                     aaGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTcg                                                                     >  2:4451010/1‑93 (MQ=255)

                      aGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTcgc                                                                    <  1:6308355/93‑1 (MQ=255)

                      aGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTcgc                                                                    >  2:1500236/1‑93 (MQ=255)

                        tGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGa                                                                  >  1:2847416/1‑93 (MQ=255)

                          aaTTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGaaa                                                                >  2:3668668/1‑93 (MQ=255)

                           aTTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAAc                                                               <  1:7505448/93‑1 (MQ=255)

                           aTTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAAc                                                               >  2:6602103/1‑93 (MQ=255)

                            ttCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACt                                                              >  1:1318032/1‑93 (MQ=255)

                              cgcgCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTgg                                                            <  2:3513671/93‑1 (MQ=255)

                                cgcTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGaa                                                          <  2:6758933/93‑1 (MQ=255)

                                cgcTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGaa                                                          <  2:2830538/93‑1 (MQ=255)

                                cgcTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGaa                                                          <  2:2620386/93‑1 (MQ=255)

                                 gcTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAAt                                                         <  1:4429409/93‑1 (MQ=255)

                                 gcTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAAt                                                         >  2:112988/1‑93 (MQ=255)

                                   ttttATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCt                                                       >  1:6940734/1‑93 (MQ=255)

                                   ttttATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCt                                                       <  1:7217139/93‑1 (MQ=255)

                                    tttATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCtt                                                      >  2:7192418/1‑93 (MQ=255)

                                    tttATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCtt                                                      <  1:5965078/93‑1 (MQ=255)

                                     ttATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCttt                                                     <  1:2425822/93‑1 (MQ=255)

                                     ttATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCttt                                                     <  1:3280319/93‑1 (MQ=255)

                                     ttATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCttt                                                     <  1:1038860/93‑1 (MQ=255)

                                       atatCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTaa                                                   <  2:5168746/93‑1 (MQ=255)

                                          tCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGa                                                <  2:4114510/93‑1 (MQ=255)

                                          tCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGa                                                >  1:1055791/1‑93 (MQ=255)

                                          tCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGa                                                >  1:346736/1‑93 (MQ=255)

                                           cTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGAc                                               <  1:6430833/93‑1 (MQ=255)

                                           cTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGAc                                               <  1:5903237/93‑1 (MQ=255)

                                            tACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGAcc                                              >  2:3856494/1‑93 (MQ=255)

                                             aCGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCg                                             >  1:121081/1‑93 (MQ=255)

                                                aaGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGc                                          >  1:6500838/1‑93 (MQ=255)

                                                aaGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGc                                          <  2:2949258/93‑1 (MQ=255)

                                                 aGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGCGTCGCGAAACTGGAATCTTTAATGACCGGGcc                                         <  2:3919551/93‑1 (MQ=255)

                                                 aGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGcc                                         <  1:6740882/93‑1 (MQ=255)

                                                 aGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGcc                                         >  1:5318841/1‑93 (MQ=255)

                                                  gCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCcg                                        >  1:2059168/1‑93 (MQ=255)

                                                  gCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCcg                                        >  2:7064181/1‑93 (MQ=255)

                                                  gCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCcg                                        >  1:6630955/1‑93 (MQ=255)

                                                    cTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTTGCGAAACTGGAATCTTTAATGACCGGGCcgcg                                      >  2:3455797/1‑93 (MQ=255)

                                                      ataAAGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTa                                    <  2:3354517/93‑1 (MQ=255)

                                                        aaaGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaa                                  >  1:1547591/1‑93 (MQ=255)

                                                        aaaGATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaa                                  <  2:6165137/93‑1 (MQ=255)

                                                         aaGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTaaaa                                 >  2:7422658/1‑93 (MQ=255)

                                                           gATGAAGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGTAATCTTTAATGACCGGGCCGCGTaaaaaa                               <  2:6705837/93‑1 (MQ=255)

                                                             tGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACg                             >  2:5024166/1‑93 (MQ=255)

                                                               aaGACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTc                           >  2:6369213/1‑93 (MQ=255)

                                                                aGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCt                          <  2:867461/93‑1 (MQ=255)

                                                                 gACGCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTg                         >  1:1658401/1‑93 (MQ=255)

                                                                  aCGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGt                        >  2:6289585/1‑93 (MQ=255)

                                                                    gCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTc                      <  1:2700759/93‑1 (MQ=255)

                                                                    gCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTc                      <  2:5552457/93‑1 (MQ=255)

                                                                    gCAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTc                      >  2:4315689/1‑93 (MQ=255)

                                                                     cAGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCa                     >  1:5257863/1‑93 (MQ=255)

                                                                      aGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTATAAAACGTCTGTTCaa                    <  1:1463547/93‑1 (MQ=255)

                                                                      aGTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCaa                    <  2:4723728/93‑1 (MQ=255)

                                                                       gTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAAt                   <  1:3207522/93‑1 (MQ=255)

                                                                       gTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAAt                   >  2:6145547/1‑93 (MQ=255)

                                                                       gTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAAt                   >  2:4107106/1‑93 (MQ=255)

                                                                       gTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAAt                   >  2:4209636/1‑93 (MQ=255)

                                                                       gTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAAt                   >  1:5045801/1‑93 (MQ=255)

                                                                       gTGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAAt                   <  2:6253581/93‑1 (MQ=255)

                                                                        tGGCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATg                  >  1:3806951/1‑93 (MQ=255)

                                                                         ggCGTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATgg                 <  1:764810/93‑1 (MQ=255)

                                                                            gTTCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGttt              <  2:4275622/93‑1 (MQ=255)

                                                                             ttCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTg             >  2:4347954/1‑93 (MQ=255)

                                                                              tCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGCCCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGa            >  1:1703515/1‑93 (MQ=255)

                                                                              tCCATAAAACCACACCCTACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGa            <  1:6200325/93‑1 (MQ=255)

                                                                              tCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGa            <  2:4792077/93‑1 (MQ=255)

                                                                              tCCATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGa            >  2:1317713/1‑93 (MQ=255)

                                                                               ccATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGAt           >  2:5920399/1‑93 (MQ=255)

                                                                               ccATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGAt           <  1:3834863/93‑1 (MQ=255)

                                                                                cATAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATg          >  1:3036111/1‑93 (MQ=255)

                                                                                 aTAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGc         <  1:4862442/93‑1 (MQ=255)

                                                                                  tAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGcc        <  2:5717464/93‑1 (MQ=255)

                                                                                  tAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGcc        <  2:6507903/93‑1 (MQ=255)

                                                                                  tAAAACCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGcc        >  2:2217536/1‑93 (MQ=255)

                                                                                    aaaCCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGt      <  1:5554025/93‑1 (MQ=255)

                                                                                      aCCACACCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGa    <  2:2097088/93‑1 (MQ=255)

                                                                                       ccACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAg   <  2:1016289/93‑1 (MQ=255)

                                                                                        cacaCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAgg  <  1:1903974/93‑1 (MQ=255)

TCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGG  >  NC_000913/1605135‑1605315

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