BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,107,299	A→G	V289A (GTG→GCG)	ymdD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,107,299	0	A	G	99.5%	601.4 / ‑5.3	192	V289A (GTG→GCG)	ymdD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (108/83); total (108/84)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 208 total aligned reads displayed.

TGAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCACATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGT  >  NC_000913/1107209‑1107388

tgAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACcgca                                                                                         <  1:6622478/93‑1 (MQ=255)

tgAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACcgca                                                                                         >  2:5444243/1‑93 (MQ=255)

tgAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACcgca                                                                                         >  1:7187019/1‑93 (MQ=255)

tgAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACcgca                                                                                         >  2:3574627/1‑93 (MQ=255)

ggAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACcgca                                                                                         <  2:5721260/92‑1 (MQ=255)

ggAACCACATAGATAAACAGCGATGCGTTAACAAAATAAATCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACcgca                                                                                         <  2:3502230/92‑1 (MQ=255)

 gAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCAt                                                                                        <  1:7625334/93‑1 (MQ=255)

  aaCCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCAtt                                                                                       <  1:2993832/93‑1 (MQ=255)

   aCCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTc                                                                                      >  2:1380662/1‑93 (MQ=255)

    ccAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCa                                                                                     >  2:1020549/1‑93 (MQ=255)

    ccAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCa                                                                                     <  2:7745807/93‑1 (MQ=255)

     cAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCAc                                                                                    <  2:3140888/93‑1 (MQ=255)

     cAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCAc                                                                                    <  2:5690639/93‑1 (MQ=255)

     cAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCAc                                                                                    >  1:7555886/1‑93 (MQ=255)

       gatagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCa                                                                                  <  2:4323466/93‑1 (MQ=255)

           gataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCCTCca                                                                              >  1:3143851/1‑93 (MQ=255)

           gataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCca                                                                              >  1:3229781/1‑93 (MQ=255)

           gataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCca                                                                              <  1:4351300/93‑1 (MQ=255)

             taAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCcaca                                                                            >  2:7867232/1‑93 (MQ=255)

              aaaCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAg                                                                           <  2:7562584/93‑1 (MQ=255)

               aaCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGa                                                                          <  2:83106/93‑1 (MQ=255)

               aaCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGa                                                                          <  2:5972498/93‑1 (MQ=255)

                aCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGAc                                                                         <  2:6992462/93‑1 (MQ=255)

                aCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGAc                                                                         <  2:1103212/93‑1 (MQ=255)

                   gCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCCCCATCCACAGACCGa                                                                      <  2:3986557/93‑1 (MQ=255)

                   gCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGa                                                                      >  1:6439552/1‑93 (MQ=255)

                    cGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAg                                                                     <  1:351596/93‑1 (MQ=255)

                      aTGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGa                                                                   <  1:2862030/93‑1 (MQ=255)

                       tGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGAc                                                                  >  2:6677905/1‑93 (MQ=255)

                       tGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGAc                                                                  >  2:4804946/1‑93 (MQ=255)

                       tGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGAc                                                                  >  2:5044517/1‑93 (MQ=255)

                        gCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGAcc                                                                 >  1:4745877/1‑93 (MQ=255)

                          gTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCAt                                                               >  2:249542/1‑93 (MQ=255)

                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTGACCATCCACAGACCGAGGACCATg                                                              >  1:5664914/1‑93 (MQ=255)

                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                              >  2:4443002/1‑93 (MQ=255)

                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                              <  1:3993633/93‑1 (MQ=255)

                             aaCAAAATAAGTCTCCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGt                                                            <  2:3153937/93‑1 (MQ=255)

                               cAAAATAAGTCCCCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGa                                                          <  1:6101153/93‑1 (MQ=255)

                               cAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGa                                                          <  1:4636954/93‑1 (MQ=255)

                                  aaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCa                                                       <  1:2742280/93‑1 (MQ=255)

                                   aTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCAc                                                      >  1:7087921/1‑93 (MQ=255)

                                   aTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCAc                                                      <  2:4537057/93‑1 (MQ=255)

                                       gTCCCCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGAc                                                  <  1:1102689/93‑1 (MQ=255)

                                          aCCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCg                                               <  1:7677687/93‑1 (MQ=255)

                                             cgcgCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGtt                                            <  2:3798168/93‑1 (MQ=255)

                                             cgcgCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGtt                                            <  2:2694848/93‑1 (MQ=255)

                                             cgcgCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGtt                                            >  1:5252474/1‑93 (MQ=255)

                                                gcTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCg                                         >  2:4115968/1‑93 (MQ=255)

                                                 cTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGt                                        >  1:4471940/1‑93 (MQ=255)

                                                 cTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGt                                        >  1:4679856/1‑93 (MQ=255)

                                                  tGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTa                                       >  2:2999196/1‑93 (MQ=255)

                                                    aCTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACa                                     <  1:682446/93‑1 (MQ=255)

                                                    aCTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACa                                     >  2:3649183/1‑93 (MQ=255)

                                                     cTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACAt                                    >  2:4916606/1‑93 (MQ=255)

                                                     cTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACAt                                    <  1:3878525/93‑1 (MQ=255)

                                                      tGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                   >  1:7937036/1‑93 (MQ=255)

                                                      tGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                   <  1:4156869/93‑1 (MQ=255)

                                                      tGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                   <  2:6326021/93‑1 (MQ=255)

                                                         aaGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAg                                <  1:1577528/93‑1 (MQ=255)

                                                         aaGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAg                                >  1:2309948/1‑93 (MQ=255)

                                                         aaGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAg                                <  2:4724427/93‑1 (MQ=255)

                                                         aaGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAg                                <  2:1213528/93‑1 (MQ=255)

                                                          aGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAgg                               >  2:863204/1‑93 (MQ=255)

                                                            ttAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCa                             >  1:1800031/1‑93 (MQ=255)

                                                            ttAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCa                             <  1:4941057/93‑1 (MQ=255)

                                                            ttAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCa                             <  2:2285945/93‑1 (MQ=255)

                                                            ttAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCa                             <  2:1751212/93‑1 (MQ=255)

                                                            ttAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCa                             <  2:5896235/93‑1 (MQ=255)

                                                             tAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCAt                            <  1:5986544/93‑1 (MQ=255)

                                                                     cGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTg                    >  1:2135401/1‑93 (MQ=255)

                                                                       gTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGcc                  >  2:2802478/1‑93 (MQ=255)

                                                                       gTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGcc                  >  2:2552408/1‑93 (MQ=255)

                                                                        gggCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCa                 <  2:7202705/92‑1 (MQ=255)

                                                                         ggCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCAt                <  2:5178259/93‑1 (MQ=255)

                                                                          gCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATa               >  1:2815234/1‑93 (MQ=255)

                                                                           gcAAAGGAGAAGCCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAg              <  2:5764349/92‑1 (MQ=255)

                                                                           ctaaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAg              <  2:5268198/91‑1 (MQ=255)

                                                                           ccAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAg              >  1:6291171/1‑93 (MQ=255)

                                                                            cAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgc             >  2:2900391/1‑93 (MQ=255)

                                                                            cAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgc             <  1:7382667/93‑1 (MQ=255)

                                                                            cAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgc             >  2:4855215/1‑93 (MQ=255)

                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg            >  2:1540195/1‑93 (MQ=255)

                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg            >  1:167039/1‑93 (MQ=255)

                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg            >  2:4835121/1‑93 (MQ=255)

                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg            >  1:6411314/1‑93 (MQ=255)

                                                                             aaaGGAGAAGACCACATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg            <  1:6642620/93‑1 (MQ=255)

                                                                              aaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcgc           <  1:174166/93‑1 (MQ=255)

                                                                               aGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCt          <  2:3411257/93‑1 (MQ=255)

                                                                                ggAGAAGCGCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGAGTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCTCTg         <  2:164662/93‑1 (MQ=255)

                                                                                ggAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTg         <  1:5083551/93‑1 (MQ=255)

                                                                                ggAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTg         <  2:6974410/93‑1 (MQ=255)

                                                                                ggAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTg         >  2:3861620/1‑93 (MQ=255)

                                                                                   gaagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGtt      <  2:1582039/93‑1 (MQ=255)

                                                                                   gaagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGtt      <  2:4634993/93‑1 (MQ=255)

                                                                                   gaagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCCGGtt      >  2:2323115/1‑93 (MQ=255)

                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg   >  1:5499010/1‑93 (MQ=255)

                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg   >  2:6771602/1‑93 (MQ=255)

                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg   >  2:2717780/1‑93 (MQ=255)

                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg   >  2:2832254/1‑93 (MQ=255)

                                                                                       aCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGt  >  2:5490757/1‑93 (MQ=255)

TGAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCACATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGT  >  NC_000913/1107209‑1107388

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