BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,460,324	A→G	D59G (GAC→GGC)	b1398 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,460,324	0	A	G	86.1%	615.6 / 81.7	259	D59G (GAC→GGC)	b1398	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (17/19); new base G (111/112); total (128/131)
Fisher's exact test for biased strand distribution p-value = 8.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 259 total aligned reads displayed.

ACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGC  >  NC_000913/1460233‑1460413

  gCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCt                                                                                        >  2:5689353/1‑93 (MQ=255)

   cTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTg                                                                                       <  2:3337779/93‑1 (MQ=255)

    tCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGc                                                                                      >  2:5515702/1‑93 (MQ=255)

     cAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCg                                                                                     >  1:442899/1‑93 (MQ=255)

     cAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCg                                                                                     >  2:589760/1‑93 (MQ=255)

     cAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCg                                                                                     <  1:7084318/93‑1 (MQ=255)

     cAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCg                                                                                     <  1:2032753/93‑1 (MQ=255)

     cAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCg                                                                                     >  1:2874555/1‑93 (MQ=255)

       aaCACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTa                                                                                   <  1:716827/93‑1 (MQ=255)

        acacGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTaa                                                                                  >  2:6113556/1‑93 (MQ=255)

         cacGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTaaa                                                                                 >  2:5739589/1‑93 (MQ=255)

         cacGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTaaa                                                                                 <  1:5063120/93‑1 (MQ=255)

         cacGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTaaa                                                                                 <  1:2793875/93‑1 (MQ=255)

          acGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAAt                                                                                <  2:1575467/93‑1 (MQ=255)

          acGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAAt                                                                                <  1:3598774/93‑1 (MQ=255)

             cgTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTc                                                                             >  2:5952673/1‑93 (MQ=255)

              gTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTcc                                                                            <  1:4593881/93‑1 (MQ=255)

               tATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTccc                                                                           <  2:5555325/93‑1 (MQ=255)

               tATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTccc                                                                           >  2:292777/1‑93 (MQ=255)

                aTGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCt                                                                          <  1:2571016/93‑1 (MQ=255)

                aTGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCt                                                                          <  1:6499726/93‑1 (MQ=255)

                 tGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTg                                                                         >  2:6232949/1‑93 (MQ=255)

                  gAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGt                                                                        <  2:2511414/93‑1 (MQ=255)

                    aaaTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTac                                                                      >  1:3372732/1‑93 (MQ=255)

                    aaaTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTac                                                                      >  2:2660721/1‑93 (MQ=255)

                     aaTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTacc                                                                     <  2:1153782/93‑1 (MQ=255)

                     aaTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTacc                                                                     >  1:4203192/1‑93 (MQ=255)

                     aaTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTacc                                                                     >  1:440400/1‑93 (MQ=255)

                      aTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTacca                                                                    >  2:4410450/1‑93 (MQ=255)

                        gtgCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTaccacc                                                                  <  1:111475/93‑1 (MQ=255)

                        gtgCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTaccacc                                                                  >  1:3202363/1‑93 (MQ=255)

                            cGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAAc                                                              >  2:1727180/1‑93 (MQ=255)

                            cGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAAc                                                              <  2:6014279/93‑1 (MQ=255)

                            cGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAAc                                                              <  2:393639/93‑1 (MQ=255)

                             gATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGGACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACa                                                             >  1:2199326/1‑93 (MQ=255)

                             gATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACa                                                             <  1:730531/93‑1 (MQ=255)

                             gATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACa                                                             >  2:5281688/1‑93 (MQ=255)

                                 tATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGAt                                                         <  1:6738524/93‑1 (MQ=255)

                                 tATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGAt                                                         <  1:1807056/93‑1 (MQ=255)

                                 tATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACGGGAt                                                         <  1:6005143/93‑1 (MQ=255)

                                  aTCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATc                                                        <  2:6882884/93‑1 (MQ=255)

                                     ggCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTgc                                                     >  1:2189734/1‑93 (MQ=255)

                                     ggCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTAGGGGACTTTTAAGGCCTGCGTAATTTTCCCTGTCCCACAAAACAGAACCTgc                                                     >  1:6141551/1‑93 (MQ=255)

                                      gcgcAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTgcg                                                    >  2:4127233/1‑93 (MQ=255)

                                      gcgcAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTgcg                                                    >  1:313006/1‑93 (MQ=255)

                                      gcgcAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTgcg                                                    >  1:2402768/1‑93 (MQ=255)

                                       cgcAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTgcgc                                                   <  1:975741/93‑1 (MQ=255)

                                        gcAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGcgcg                                                  >  2:1382228/1‑93 (MQ=255)

                                          aaaTTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGAc                                                <  2:2175646/93‑1 (MQ=255)

                                           aaTTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACa                                               >  2:6868225/1‑93 (MQ=255)

                                           aaTTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACa                                               >  1:1032160/1‑93 (MQ=255)

                                             ttCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGATAAc                                             >  2:6397280/1‑93 (MQ=255)

                                             ttCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAAc                                             >  1:1971940/1‑93 (MQ=255)

                                               cGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTa                                           <  2:5122768/93‑1 (MQ=255)

                                               cGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTa                                           >  2:1218056/1‑93 (MQ=255)

                                                gACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTAt                                          >  2:3563740/1‑93 (MQ=255)

                                                gACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTAt                                          >  1:3569448/1‑93 (MQ=255)

                                                 aCGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATc                                         <  2:6638928/93‑1 (MQ=255)

                                                   gcagcaGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATccc                                       <  1:1354762/93‑1 (MQ=255)

                                                    cagcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCt                                      <  2:7504392/93‑1 (MQ=255)

                                                     agcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCtt                                     >  1:1492003/1‑93 (MQ=255)

                                                      gcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCttt                                    <  1:5230820/93‑1 (MQ=255)

                                                        agGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGa                                  >  1:3630872/1‑93 (MQ=255)

                                                          gCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGaca                                >  2:5947159/1‑93 (MQ=255)

                                                           cGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGacac                               <  2:1931463/93‑1 (MQ=255)

                                                            gTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACAcc                              >  2:5609054/1‑93 (MQ=255)

                                                             tACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCt                             >  1:4863065/1‑93 (MQ=255)

                                                             tACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCt                             <  1:6986891/93‑1 (MQ=255)

                                                             tACATCCTGATGATTTCAGGGAACTTTCAGACCTGTGTAAATTTCCCTGTACCACCAAACGGGATCTGCGCGACAACTATCCCTTTGACACCt                             >  2:364608/1‑93 (MQ=255)

                                                              aCATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCtt                            <  1:910396/93‑1 (MQ=255)

                                                              aCATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCtt                            <  1:5913215/93‑1 (MQ=255)

                                                               cATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCttt                           >  2:1651126/1‑93 (MQ=255)

                                                                 tCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGc                         <  2:6760738/93‑1 (MQ=255)

                                                                 tCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGc                         >  2:5835539/1‑93 (MQ=255)

                                                                    tgatgaTTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGt                      <  1:6089663/93‑1 (MQ=255)

                                                                    tgatgaTTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGt                      <  1:4967776/93‑1 (MQ=255)

                                                                    tgatgaTTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGt                      >  1:571571/1‑93 (MQ=255)

                                                                     gatgatTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTg                     >  1:3863639/1‑93 (MQ=255)

                                                                     gatgatTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTg                     >  2:7279611/1‑93 (MQ=255)

                                                                      atgatTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGc                    <  2:1262103/93‑1 (MQ=255)

                                                                       tgatATCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGcc                   >  1:435169/1‑93 (MQ=255)

                                                                         atTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCaa                 <  2:703618/93‑1 (MQ=255)

                                                                          tttCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAAt                >  1:5889675/1‑93 (MQ=255)

                                                                          tttCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAAt                <  2:5486859/93‑1 (MQ=255)

                                                                           ttCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATg               <  2:1899787/93‑1 (MQ=255)

                                                                              aGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGaa            >  1:5521800/1‑93 (MQ=255)

                                                                              aGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGaa            >  1:5820490/1‑93 (MQ=255)

                                                                              aGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGaa            <  2:6571943/93‑1 (MQ=255)

                                                                               gggAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAAc           <  2:383807/93‑1 (MQ=255)

                                                                               gggAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAAc           >  1:6553763/1‑93 (MQ=255)

                                                                               gggAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAAc           >  1:1300428/1‑93 (MQ=255)

                                                                                ggAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACa          >  2:3969545/1‑93 (MQ=255)

                                                                                 gAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACaa         <  2:4727872/93‑1 (MQ=255)

                                                                                 gAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACaa         <  2:1506968/93‑1 (MQ=255)

                                                                                  aaCTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAg        >  1:2639555/1‑93 (MQ=255)

                                                                                     tttCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAgtgg     <  1:4031637/93‑1 (MQ=255)

                                                                                      ttCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAgtggt    <  1:1920777/93‑1 (MQ=255)

                                                                                      ttCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAgtggt    <  1:367260/93‑1 (MQ=255)

                                                                                        cAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTgc  <  2:4066824/93‑1 (MQ=255)

                                                                                        cAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTgc  <  1:630426/93‑1 (MQ=255)

ACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGC  >  NC_000913/1460233‑1460413

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