BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,924,185	A→G	G14G (GGA→GGG)	b1844 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,924,185	0	A	G	96.8%	369.9 / 1.0	123	G14G (GGA→GGG)	b1844	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/3); new base G (68/51); total (69/54)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 134 total aligned reads displayed.

AAGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTATTA  >  NC_000913/1924095‑1924273

aaGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCagggggg                                                                                        <  2:1664295/93‑1 (MQ=255)

aaGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCagggggg                                                                                        <  2:1236193/93‑1 (MQ=255)

aaGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCagggagg                                                                                        >  1:4031916/1‑93 (MQ=255)

 aGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCaggggggg                                                                                       <  1:1089931/93‑1 (MQ=255)

 aGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCaggggggg                                                                                       >  1:2882696/1‑93 (MQ=255)

  ggCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAggggggga                                                                                      >  2:461964/1‑93 (MQ=255)

   gCGATATCATTCCATTACGCTAGGCTTTTGCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGAt                                                                                     <  2:4456223/93‑1 (MQ=255)

    cGATATCATTCCATTACGCTAGGCTTTTTCGGCCTAGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATc                                                                                    <  2:3909902/93‑1 (MQ=255)

     gATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCg                                                                                   >  1:1436004/1‑93 (MQ=255)

      atatCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGt                                                                                  >  1:1287925/1‑93 (MQ=255)

      atatCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGt                                                                                  >  1:3277942/1‑93 (MQ=255)

      atatCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGt                                                                                  >  1:1925652/1‑93 (MQ=255)

       tatCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGtt                                                                                 >  1:2640742/1‑93 (MQ=255)

           aTTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGGTCGTTgaga                                                                             >  2:2547095/1‑93 (MQ=255)

           aTTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTgaga                                                                             <  2:3308234/93‑1 (MQ=255)

           aTTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTgaga                                                                             <  1:827970/93‑1 (MQ=255)

            ttCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGAt                                                                            <  1:3016406/93‑1 (MQ=255)

                aTTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGcc                                                                        >  1:3905075/1‑93 (MQ=255)

                aTTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGcc                                                                        <  1:134711/93‑1 (MQ=255)

                 ttgcGCTAGGCTTGTTCGGCCTGGAGCATGCCTTGTTGCGCAGTATCGATCCAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCct                                                                       <  1:343875/90‑1 (MQ=255)

                  tACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCctc                                                                      <  1:2529167/93‑1 (MQ=255)

                     gCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGt                                                                   >  2:3366914/1‑93 (MQ=255)

                     gCTAGGCTTGTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGt                                                                   <  1:542372/93‑1 (MQ=255)

                      cTAGGCTTTTTCGGCCTGGAGCATTCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGtt                                                                  <  2:483921/93‑1 (MQ=255)

                        aGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGa                                                                <  2:3576137/93‑1 (MQ=255)

                        aGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGa                                                                >  1:4257594/1‑93 (MQ=255)

                        aGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGa                                                                <  1:1576825/93‑1 (MQ=255)

                         ggCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGAt                                                               <  1:623320/93‑1 (MQ=255)

                         ggCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGAt                                                               <  2:3618825/93‑1 (MQ=255)

                         ggCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGAt                                                               >  1:2233136/1‑93 (MQ=255)

                         ggCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGAt                                                               >  1:230634/1‑93 (MQ=255)

                         ggCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGCTGAt                                                               >  2:3064888/1‑93 (MQ=255)

                          gCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATg                                                              <  1:3832569/93‑1 (MQ=255)

                            tttttCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTc                                                            <  1:4434786/93‑1 (MQ=255)

                            tttttCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTc                                                            >  2:332848/1‑93 (MQ=255)

                            tttttCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTc                                                            >  2:1828229/1‑93 (MQ=255)

                               ttCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCAtt                                                         <  1:2584081/93‑1 (MQ=255)

                               ttCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCAtt                                                         <  2:2903321/93‑1 (MQ=255)

                                tCGTCCTTGAGCATGCCATGTTGCGCATTATCGATACAGAAAGCTGCGGTTTGCAGGTATGGATCTTTGAGATTGCCTCTGTTTATGTCATTg                                                        <  2:4334376/93‑1 (MQ=255)

                                  ggCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGAc                                                      <  1:4615282/93‑1 (MQ=255)

                                   gCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACg                                                     <  1:1617508/93‑1 (MQ=255)

                                     cTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGa                                                   >  1:1323806/1‑93 (MQ=255)

                                     cTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGa                                                   >  1:556630/1‑93 (MQ=255)

                                      tGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGaa                                                  >  2:2795498/1‑93 (MQ=255)

                                       ggAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGaaa                                                 >  2:4725725/1‑93 (MQ=255)

                                       ggAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGaaa                                                 >  1:2286098/1‑93 (MQ=255)

                                         aGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGaaaaa                                               >  2:1044626/1‑93 (MQ=255)

                                         aGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGaaaaa                                               >  1:1435621/1‑93 (MQ=255)

                                          gCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAACAAt                                              >  2:1988488/1‑93 (MQ=255)

                                          gCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAAt                                              >  2:2858353/1‑93 (MQ=255)

                                          gCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAAt                                              >  1:1548716/1‑93 (MQ=255)

                                           catgccatgcTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAACAATc                                             >  2:1787110/1‑93 (MQ=255)

                                           cATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATc                                             <  2:3047656/93‑1 (MQ=255)

                                           cATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATc                                             <  1:2865876/93‑1 (MQ=255)

                                            aTGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCg                                            >  1:1806105/1‑93 (MQ=255)

                                            aTGCCATGTTGCGCATTATCGATACAGAAACCTGCGGGTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGCAAAATCg                                            >  2:1937900/1‑93 (MQ=255)

                                             tGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGt                                           >  2:784505/1‑93 (MQ=255)

                                             tGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGt                                           >  1:3740968/1‑93 (MQ=255)

                                                cATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCaa                                        >  1:45218/1‑93 (MQ=255)

                                                cATGTTGCGCATTATCGATACAGAAACCTGCGGTGTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCaa                                        >  2:2978393/1‑93 (MQ=255)

                                                cATGTTGCGCATTATCGATACAGAAACCTGCGGGTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCaa                                        >  1:2860181/1‑93 (MQ=255)

                                                cATGTTGCGCATTATCGATACAGAAAACTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCATCaa                                        >  1:558257/1‑93 (MQ=255)

                                                 aTGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAAc                                       >  1:2625770/1‑93 (MQ=255)

                                                 aTGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAAc                                       >  1:2625896/1‑93 (MQ=255)

                                                   gTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAAccc                                     <  2:1942445/93‑1 (MQ=255)

                                                    ttGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAAcccc                                    <  1:833734/93‑1 (MQ=255)

                                                     tGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCa                                   <  1:360512/93‑1 (MQ=255)

                                                     tGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCa                                   <  2:2247963/93‑1 (MQ=255)

                                                      gcgcATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCAt                                  <  2:762459/93‑1 (MQ=255)

                                                        gcATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGa                                >  2:2774103/1‑93 (MQ=255)

                                                         cATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAg                               >  1:154262/1‑93 (MQ=255)

                                                         cATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAAACCCATgcg                               >  2:1800721/1‑91 (MQ=255)

                                                           ttATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGcc                             <  2:2915836/93‑1 (MQ=255)

                                                             aTCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCAAc                           >  1:630827/1‑93 (MQ=255)

                                                              tCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCAcc                          <  1:3378042/93‑1 (MQ=255)

                                                              tCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCAcc                          >  2:1698971/1‑93 (MQ=255)

                                                              tCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCAcc                          <  2:2054491/93‑1 (MQ=255)

                                                                gATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTg                        <  1:3483758/93‑1 (MQ=255)

                                                                gATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTg                        >  1:2140558/1‑93 (MQ=255)

                                                                  tACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGt                      <  1:4344456/93‑1 (MQ=255)

                                                                  tACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGt                      <  2:2282532/93‑1 (MQ=255)

                                                                   aCAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTTATGTCATTGACGGAAAAATAGTCAACACCATGAGCCACACGGTg                     >  2:4616993/1‑93 (MQ=255)

                                                                      gAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTgcgc                  >  2:4093366/1‑93 (MQ=255)

                                                                       aaaCCTGCGGTGTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGcc                 <  1:1819246/93‑1 (MQ=255)

                                                                        aaCCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGccc                <  1:2045285/93‑1 (MQ=255)

                                                                        aaCCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGccc                >  2:4010628/1‑93 (MQ=255)

                                                                          ccTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACTCCATGAGCCACCTGGTGCGCCCCg              <  2:4555648/93‑1 (MQ=255)

                                                                           cTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGa             >  2:4680276/1‑93 (MQ=255)

                                                                            tGCGGTTTGCAGGGGGGGCTCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCACCGAt            >  2:2266014/1‑93 (MQ=255)

                                                                            tGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGAt            >  2:1024045/1‑93 (MQ=255)

                                                                            tGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTAAACCCCATGAGCCACCTGGTGCGCCCCGAt            >  2:4512940/1‑93 (MQ=255)

                                                                             gCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATc           <  2:3648157/93‑1 (MQ=255)

                                                                                gTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTc        >  1:1029290/1‑93 (MQ=255)

                                                                                 tttGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTcc       >  2:20574/1‑93 (MQ=255)

                                                                                   tGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTa     <  1:800148/93‑1 (MQ=255)

                                                                                     cAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTAtt   <  2:1836307/93‑1 (MQ=255)

                                                                                     cAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTAtt   <  1:2532203/93‑1 (MQ=255)

                                                                                      agggggggATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTATTa  >  2:1703340/1‑93 (MQ=255)

                                                                                      agggggggATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTATTa  >  1:1434787/1‑93 (MQ=255)

                                                                                      agggggggATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTATTa  >  1:1056744/1‑93 (MQ=255)

AAGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTATTA  >  NC_000913/1924095‑1924273

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