BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,049,140	A→G	I72V (ATT→GTT)	yihI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,049,140	0	A	G	99.6%	778.2 / NA	242	I72V (ATT→GTT)	yihI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (126/115); minor base C (1/0); total (127/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

CCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTA  >  NC_000913/4049050‑4049229

 cGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTc                                                                                        <  2:3102105/93‑1 (MQ=255)

  gggCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTcc                                                                                       >  1:4388530/1‑93 (MQ=255)

  gggCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTcc                                                                                       >  2:4499404/1‑93 (MQ=255)

   ggCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCa                                                                                      <  2:1089923/93‑1 (MQ=255)

   ggCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCa                                                                                      >  2:3972874/1‑93 (MQ=255)

    gCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAATCCCCTGTTCCAt                                                                                     <  2:1601652/93‑1 (MQ=255)

    gCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCAt                                                                                     <  1:2392812/93‑1 (MQ=255)

    gCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCAt                                                                                     <  2:2491405/93‑1 (MQ=255)

        ccGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTgag                                                                                 >  2:1462944/1‑91 (MQ=255)

        ccGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTgag                                                                                 <  2:358102/93‑3 (MQ=255)

         cgcgCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGc                                                                                <  1:461988/93‑1 (MQ=255)

         cgcgCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGc                                                                                <  2:3394546/93‑1 (MQ=255)

            gcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTg                                                                             >  2:240868/1‑93 (MQ=255)

              aGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGAc                                                                           <  1:2990347/93‑1 (MQ=255)

               gcgggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACt                                                                          >  1:225550/1‑93 (MQ=255)

                 gggcggAAACACCACGTCTGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                        >  1:3375206/1‑93 (MQ=255)

                 gggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                        >  1:1027453/1‑93 (MQ=255)

                 gggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                        <  1:2541141/93‑1 (MQ=255)

                     ggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaaa                                                                    >  1:1927526/1‑93 (MQ=255)

                     ggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaaa                                                                    <  2:185609/93‑1 (MQ=255)

                     ggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaaa                                                                    <  1:3799834/93‑1 (MQ=255)

                     ggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaaa                                                                    <  1:4221461/93‑1 (MQ=255)

                      gAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAg                                                                   >  2:2570935/1‑93 (MQ=255)

                        aaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTc                                                                 <  1:659936/93‑1 (MQ=255)

                        aaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTc                                                                 >  2:4380726/1‑93 (MQ=255)

                            ccacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCa                                                             >  1:3063472/1‑93 (MQ=255)

                             cacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCaa                                                            <  2:924787/93‑1 (MQ=255)

                              acGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCaaa                                                           <  1:2463092/93‑1 (MQ=255)

                               cGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAAc                                                          >  2:728659/1‑93 (MQ=255)

                                 tCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAg                                                        >  2:2414172/1‑93 (MQ=255)

                                 tCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAg                                                        >  2:4214125/1‑93 (MQ=255)

                                   aGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGca                                                      <  1:61693/93‑1 (MQ=255)

                                   aGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGca                                                      >  2:2884879/1‑93 (MQ=255)

                                    ggCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGcac                                                     <  1:3599078/93‑1 (MQ=255)

                                      cagcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACaa                                                   <  2:943262/93‑1 (MQ=255)

                                        gcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAc                                                 >  2:4411526/1‑93 (MQ=255)

                                         caAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAcc                                                <  1:4008187/93‑1 (MQ=255)

                                          aaaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCg                                               >  1:1316596/1‑93 (MQ=255)

                                           aaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGa                                              <  1:2066891/93‑1 (MQ=255)

                                            aGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                                             >  2:352930/1‑93 (MQ=255)

                                            aGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGCAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                                             <  1:2866028/93‑1 (MQ=255)

                                             ggCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAag                                            >  2:548829/1‑93 (MQ=255)

                                              gCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAaga                                           <  1:728212/93‑1 (MQ=255)

                                               ccAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAagag                                          <  1:3151072/93‑1 (MQ=255)

                                                cAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGt                                         <  1:2607253/93‑1 (MQ=255)

                                                 aGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTg                                        >  1:350365/1‑93 (MQ=255)

                                                 aGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTg                                        >  1:3456246/1‑93 (MQ=255)

                                                  gAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTga                                       >  1:4280881/1‑93 (MQ=255)

                                                  gAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTga                                       <  1:3214100/93‑1 (MQ=255)

                                                    aCGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTgaga                                     <  2:224582/93‑1 (MQ=255)

                                                        aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                 >  2:3668802/1‑93 (MQ=255)

                                                        aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                 >  1:4444680/1‑93 (MQ=255)

                                                        aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                 >  2:4477323/1‑93 (MQ=255)

                                                         ccAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCt                                >  1:1261109/1‑93 (MQ=255)

                                                          cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                               >  2:3159436/1‑93 (MQ=255)

                                                          cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                               >  1:177820/1‑93 (MQ=255)

                                                          cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                               <  2:3797468/93‑1 (MQ=255)

                                                          cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                               >  2:3706876/1‑93 (MQ=255)

                                                          cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                               >  1:2727507/1‑93 (MQ=255)

                                                          cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                               <  1:1152289/93‑1 (MQ=255)

                                                           aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                              <  1:2622395/93‑1 (MQ=255)

                                                           aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                              <  2:1632576/93‑1 (MQ=255)

                                                           aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                              >  1:3939618/1‑93 (MQ=255)

                                                            aaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATg                             >  2:3543826/1‑93 (MQ=255)

                                                             aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGc                            >  1:1182702/1‑93 (MQ=255)

                                                             aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGc                            <  2:2087176/93‑1 (MQ=255)

                                                             aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGc                            >  1:3359180/1‑93 (MQ=255)

                                                              aGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCt                           <  1:4645697/93‑1 (MQ=255)

                                                               gATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCtt                          >  1:1344867/1‑93 (MQ=255)

                                                               gATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCtt                          >  1:1481269/1‑93 (MQ=255)

                                                               gATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCtt                          >  2:2122817/1‑93 (MQ=255)

                                                                aTCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCttt                         <  2:403333/93‑1 (MQ=255)

                                                                aTCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCttt                         <  1:4513254/93‑1 (MQ=255)

                                                                  ccACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCa                       <  2:2900129/93‑1 (MQ=255)

                                                                  ccACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCa                       >  2:3343381/1‑93 (MQ=255)

                                                                  ccACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCa                       >  2:1816471/1‑93 (MQ=255)

                                                                    aCGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCAcc                     >  1:2587326/1‑93 (MQ=255)

                                                                       tATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCa                  <  1:2831884/93‑1 (MQ=255)

                                                                       tATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCa                  >  1:1098197/1‑93 (MQ=255)

                                                                        aTTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAg                 <  1:3875495/93‑1 (MQ=255)

                                                                         ttGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAgg                >  1:99696/1‑93 (MQ=255)

                                                                         ttGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAgg                <  1:2895200/93‑1 (MQ=255)

                                                                         ttGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAgg                >  1:2088473/1‑93 (MQ=255)

                                                                         ctgacAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAgg                <  2:4100409/89‑1 (MQ=255)

                                                                          tGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGc               <  2:1399595/93‑1 (MQ=255)

                                                                           ggCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCg              >  1:92532/1‑93 (MQ=255)

                                                                           ggCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCg              <  1:1150121/93‑1 (MQ=255)

                                                                             cAGTAAACCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGa            <  1:1409284/93‑1 (MQ=255)

                                                                              aGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAg           >  2:491330/1‑93 (MQ=255)

                                                                              aGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAg           <  1:784087/93‑1 (MQ=255)

                                                                               gTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGt          <  1:1117135/93‑1 (MQ=255)

                                                                                 aaaaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTg        >  1:2713517/1‑93 (MQ=255)

                                                                                 aaaaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTg        >  1:2433780/1‑93 (MQ=255)

                                                                                 aaaaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTg        <  1:1692368/93‑1 (MQ=255)

                                                                                 aaaaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTg        <  1:342905/93‑1 (MQ=255)

                                                                                   aaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGa      >  1:1940420/1‑93 (MQ=255)

                                                                                     ccccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAt    >  1:2925082/1‑93 (MQ=255)

                                                                                     ccccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAt    <  2:2289008/93‑1 (MQ=255)

                                                                                       ccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTa  >  1:880493/1‑93 (MQ=255)

                                                                                       ccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTa  >  2:865176/1‑93 (MQ=255)

CCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTA  >  NC_000913/4049050‑4049229

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