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.6%	708.7 / ‑5.3	221	V289A (GTG→GCG)	ymdD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (123/97); total (124/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 234 total aligned reads displayed.

TGAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCACATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGT  >  NC_000913/1107209‑1107388

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

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

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

 gAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCAt                                                                                        >  1:2144747/1‑93 (MQ=255)

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

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

      agatagatAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCAcc                                                                                   >  2:1895571/1‑93 (MQ=255)

      agatagatAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTAAcc                                                                                   >  1:2103555/1‑93 (MQ=255)

       gatagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCa                                                                                  >  1:984780/1‑93 (MQ=255)

         tagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATc                                                                                <  1:1409592/93‑1 (MQ=255)

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

            ataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCcac                                                                             >  1:3307652/1‑93 (MQ=255)

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

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

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

               aaCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGa                                                                          >  1:194041/1‑93 (MQ=255)

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

                 cAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGAcc                                                                        <  2:5104662/93‑1 (MQ=255)

                 cAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGAcc                                                                        >  2:2716089/1‑93 (MQ=255)

                  aGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCg                                                                       <  2:3854602/93‑1 (MQ=255)

                  aGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCg                                                                       >  2:5473950/1‑93 (MQ=255)

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

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

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

                       tGCGTTAACAAAATAAGTCTCCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGAc                                                                  >  1:4211703/1‑93 (MQ=255)

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

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

                        gCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGAcc                                                                 <  1:1160213/93‑1 (MQ=255)

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

                          gTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCAt                                                               <  2:3996003/93‑1 (MQ=255)

                          gTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCAt                                                               <  1:2034509/93‑1 (MQ=255)

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

                             aaCAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGt                                                            >  1:628915/1‑93 (MQ=255)

                              aCAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTg                                                           >  1:100987/1‑93 (MQ=255)

                              aCAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTg                                                           >  1:5149701/1‑93 (MQ=255)

                                aaaaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGAt                                                         <  1:52065/93‑1 (MQ=255)

                                 aaaTAAGTCACCCGCGCTGGCTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATc                                                        <  1:5006187/93‑1 (MQ=255)

                                 aaaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATc                                                        <  2:3734478/93‑1 (MQ=255)

                                  aaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCa                                                       >  2:4869907/1‑93 (MQ=255)

                                    tAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCAcc                                                     <  2:1293085/93‑1 (MQ=255)

                                     aaGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCg                                                    >  1:2284298/1‑93 (MQ=255)

                                      aGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGa                                                   <  2:572686/93‑1 (MQ=255)

                                       gTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGAc                                                  <  1:4714362/93‑1 (MQ=255)

                                       gTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGAc                                                  <  2:2391261/93‑1 (MQ=255)

                                         cACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTc                                                <  2:2368002/93‑1 (MQ=255)

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

                                          aCCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCg                                               >  2:4470932/1‑93 (MQ=255)

                                          aCCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCACATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCg                                               >  2:116968/1‑93 (MQ=255)

                                           cccGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCgg                                              >  2:4026041/1‑93 (MQ=255)

                                           cccGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCgg                                              >  1:1663498/1‑93 (MQ=255)

                                            ccGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGt                                             <  1:3521924/93‑1 (MQ=255)

                                              gcgcTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGttt                                           <  2:2499437/93‑1 (MQ=255)

                                              gcgcTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGttt                                           >  1:3969555/1‑93 (MQ=255)

                                              gcgcTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGctt                                           >  2:2572768/1‑93 (MQ=255)

                                              gcgcTGACTAGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGttt                                           <  2:422130/93‑1 (MQ=255)

                                                gcTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCg                                         <  2:591186/93‑1 (MQ=255)

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

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

                                                 cTGACTGGAAGATAAGCAAACGGTGACCAAAGGAGAAGGCCGCCTTAACCATCCACAGACCGAGGACCATGGTGATCACCGGCTCGGTTTCGt                                        <  1:3234257/93‑1 (MQ=255)

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

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

                                                       ggAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATcc                                  >  2:1360956/1‑93 (MQ=255)

                                                       ggAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATcc                                  <  2:1140395/93‑1 (MQ=255)

                                                       ggAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATcc                                  <  1:6015300/93‑1 (MQ=255)

                                                        gAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCa                                 >  1:2012837/1‑93 (MQ=255)

                                                        gAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCa                                 >  2:3756149/1‑93 (MQ=255)

                                                        gAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGCGATCACCGACTCGGTTTCGTACATCCa                                 >  2:277497/1‑93 (MQ=255)

                                                           gTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGc                              >  1:5311595/1‑93 (MQ=255)

                                                             tAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCAt                            >  1:1369373/1‑93 (MQ=255)

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

                                                              aaGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATc                           <  1:439545/93‑1 (MQ=255)

                                                               aGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCg                          >  1:5729502/1‑93 (MQ=255)

                                                                 cAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGcc                        <  2:538423/93‑1 (MQ=255)

                                                                 cAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGcc                        >  2:3074901/1‑93 (MQ=255)

                                                                  aaaCGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCa                       >  2:1737085/1‑93 (MQ=255)

                                                                   aaCGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATGGCCAc                      <  1:1197075/93‑1 (MQ=255)

                                                                   aaCGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCAc                      >  2:3307743/1‑93 (MQ=255)

                                                                    aCGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACt                     >  2:3169690/1‑93 (MQ=255)

                                                                    aCGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACt                     <  2:3692340/93‑1 (MQ=255)

                                                                       gTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGcc                  <  1:2404793/93‑1 (MQ=255)

                                                                       gTGGCCAAAGGAGAAGACCGCATACACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGcc                  >  2:5072622/1‑93 (MQ=255)

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

                                                                         ggCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCAt                >  2:3280999/1‑93 (MQ=255)

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

                                                                          gCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATa               <  2:4333100/93‑1 (MQ=255)

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

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

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

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

                                                                               aGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCt          >  1:4111452/1‑93 (MQ=255)

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

                                                                                 gagaAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTgg        <  1:6042032/93‑1 (MQ=255)

                                                                                   gaagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGtt      >  2:4666242/1‑93 (MQ=255)

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

                                                                                    aagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTg     >  1:3595821/1‑93 (MQ=255)

                                                                                     agaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGa    <  1:4178887/93‑1 (MQ=255)

                                                                                     agaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGa    >  2:1655690/1‑93 (MQ=255)

                                                                                     agaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGa    <  2:1396375/93‑1 (MQ=255)

                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTGCATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg   >  2:1685505/1‑93 (MQ=255)

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

TGAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCACATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGT  >  NC_000913/1107209‑1107388

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