BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,400,388	A→G	S298P (TCC→CCC)	b1337 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,400,388	0	A	G	99.4%	555.8 / ‑5.3	178	S298P (TCC→CCC)	b1337	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (90/87); total (91/87)
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.91e-01

Only 100 of 185 total aligned reads displayed.

CCAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCT  >  NC_000913/1400298‑1400477

ccAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGAc                                                                                         >  1:6902144/1‑93 (MQ=255)

 cAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACa                                                                                        >  2:557518/1‑93 (MQ=255)

 cAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACa                                                                                        <  1:4824511/93‑1 (MQ=255)

 cAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACa                                                                                        <  1:453986/93‑1 (MQ=255)

  aGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAg                                                                                       >  2:925725/1‑93 (MQ=255)

  aGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAg                                                                                       <  2:4338419/93‑1 (MQ=255)

  aGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAg                                                                                       >  2:3136697/1‑93 (MQ=255)

  aGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAg                                                                                       <  2:1964550/93‑1 (MQ=255)

  aGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAg                                                                                       <  1:5860470/93‑1 (MQ=255)

   gTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAgt                                                                                      >  1:2386174/1‑92 (MQ=255)

   gTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAgg                                                                                      >  2:6510664/1‑93 (MQ=255)

   gTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAgg                                                                                      <  1:25383/93‑1 (MQ=255)

    ttcttcGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTGCATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGc                                                                                     >  1:1185933/1‑93 (MQ=255)

     tcttcGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCt                                                                                    >  1:6386381/1‑93 (MQ=255)

      cttcGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCtt                                                                                   >  2:5077037/1‑93 (MQ=255)

       ttcGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCttt                                                                                  <  2:5527139/93‑1 (MQ=255)

       ttcGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCttt                                                                                  <  2:2003800/93‑1 (MQ=255)

       ttcGGAGTTCCATTCCGGAGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCttt                                                                                  >  2:7904977/1‑93 (MQ=255)

         cGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGt                                                                                <  2:3984960/93‑1 (MQ=255)

          ggAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTc                                                                               <  2:4495518/93‑1 (MQ=255)

             gTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATCGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGaa                                                                            >  2:4792749/1‑93 (MQ=255)

              ttCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAg                                                                           >  2:5790457/1‑93 (MQ=255)

                 cATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCgg                                                                        >  2:1985524/1‑93 (MQ=255)

                   ttCCGGGGTACCCAACTGGGATAGGGCCTGGTACATGGCATTTTCTAGGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGct                                                                      <  2:7024824/93‑2 (MQ=255)

                   ttCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCa                                                                      <  1:3577433/93‑1 (MQ=255)

                   ttCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCa                                                                      <  2:2098921/93‑1 (MQ=255)

                     ccGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCAtt                                                                    <  2:7667865/93‑1 (MQ=255)

                      cGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTc                                                                   <  1:4313973/93‑1 (MQ=255)

                      cGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTc                                                                   <  2:7285306/93‑1 (MQ=255)

                       ggggTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCa                                                                  >  2:7800727/1‑93 (MQ=255)

                       ggggTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCa                                                                  <  2:300339/93‑1 (MQ=255)

                       ggggTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCa                                                                  >  2:5602279/1‑93 (MQ=255)

                       ggggTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCTGCATTCa                                                                  >  2:5193018/1‑93 (MQ=255)

                           tACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGGCAGGCTTTGTCGAAGCGGCATTCAACCg                                                              <  2:4723457/93‑1 (MQ=255)

                           tACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCg                                                              >  1:5728574/1‑93 (MQ=255)

                           tACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCg                                                              <  2:7508992/93‑1 (MQ=255)

                            aCCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCgg                                                             >  1:1531575/1‑92 (MQ=255)

                             ccAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTg                                                            <  1:219403/93‑1 (MQ=255)

                             ccAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTg                                                            >  1:1916173/1‑93 (MQ=255)

                             ccAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTg                                                            <  1:3436852/93‑1 (MQ=255)

                               aaaaTGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGCGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGt                                                          <  1:5615373/93‑1 (MQ=255)

                                aaaTGGGATAGGGCCTGGGACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGtt                                                         >  2:6288164/1‑93 (MQ=255)

                                 aaTGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGttt                                                        >  2:3687550/1‑93 (MQ=255)

                                  aTGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTc                                                       >  1:1276262/1‑93 (MQ=255)

                                  aTGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTc                                                       <  2:1475352/93‑1 (MQ=255)

                                     ggATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGt                                                    >  2:6057088/1‑93 (MQ=255)

                                       aTAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTGGGtca                                                  >  1:2270359/1‑93 (MQ=255)

                                       aTAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGtca                                                  >  2:6271013/1‑93 (MQ=255)

                                         aGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGtcatc                                                >  1:7211041/1‑93 (MQ=255)

                                              cTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGc                                           <  1:1925213/93‑1 (MQ=255)

                                              cTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGc                                           >  2:2409816/1‑93 (MQ=255)

                                               tGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGcc                                          >  2:7480933/1‑93 (MQ=255)

                                               tGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGcc                                          <  1:293823/93‑1 (MQ=255)

                                                ggTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCg                                         <  1:415011/93‑1 (MQ=255)

                                                 gTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGc                                        >  1:5488037/1‑93 (MQ=255)

                                                    cATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCAcc                                     <  1:13791/93‑1 (MQ=255)

                                                     aTGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCt                                    >  1:4996948/1‑93 (MQ=255)

                                                     aTGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCt                                    >  2:1322464/1‑93 (MQ=255)

                                                     aTGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCt                                    >  1:5551843/1‑93 (MQ=255)

                                                      tGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCtt                                   >  2:26545/1‑93 (MQ=255)

                                                       ggCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTc                                  >  1:4387616/1‑93 (MQ=255)

                                                       ggCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTc                                  <  1:7521454/93‑1 (MQ=255)

                                                       ggCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTc                                  <  1:6337090/93‑1 (MQ=255)

                                                       ggCATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTc                                  <  2:3616581/93‑1 (MQ=255)

                                                         cATTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCgg                                <  2:4879651/93‑1 (MQ=255)

                                                          aTTTTCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGTATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGc                               >  2:3084395/1‑93 (MQ=255)

                                                           ttttCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCg                              >  2:6824164/1‑93 (MQ=255)

                                                           ttttCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCg                              >  1:1612594/1‑93 (MQ=255)

                                                             ttCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAt                            >  2:2459537/1‑93 (MQ=255)

                                                             ttCTAAGGTGCGATTCGGCAGATCACTGGGACAGGCTTTGTCGAAGCGGCATTCACCCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAt                            <  2:3418615/93‑1 (MQ=255)

                                                              tCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAtt                           <  2:6780290/93‑1 (MQ=255)

                                                              tCTAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAtt                           >  1:3230121/1‑93 (MQ=255)

                                                                tAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGTATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAtttt                         <  1:4940574/93‑1 (MQ=255)

                                                                tAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAtttt                         >  1:3134605/1‑93 (MQ=255)

                                                                tAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAtttt                         <  1:4244784/93‑1 (MQ=255)

                                                                tAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAtttt                         >  2:2624408/1‑93 (MQ=255)

                                                                tAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAtttt                         >  2:1791331/1‑93 (MQ=255)

                                                                tAAGGTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAtttt                         >  2:1339724/1‑93 (MQ=255)

                                                                    gTGCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCg                     >  1:2318046/1‑93 (MQ=255)

                                                                      gCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGAc                   >  2:1543337/1‑93 (MQ=255)

                                                                      gCGATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGCCATCAATGCCGCACCTTCGGCGATTTTGGCGAc                   >  1:7705970/1‑93 (MQ=255)

                                                                        gATTCGGGAGATAACTGGGCCAGGCTTTGTCGAGGCGGCATTCACCCGTGGTTTCGGGCATCACTGGCGCACTTTCGGCGATTTTGGCGcccc                 >  1:2790907/1‑93 (MQ=255)

                                                                        gATTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGAccc                 >  1:2790908/1‑93 (MQ=255)

                                                                         aTTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCg                >  1:1106344/1‑93 (MQ=255)

                                                                         aTTCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCg                <  2:4110717/93‑1 (MQ=255)

                                                                          ttCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGa               >  1:3665559/1‑93 (MQ=255)

                                                                          ttCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGa               >  1:3665560/1‑93 (MQ=255)

                                                                          ttCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGa               >  1:6676824/1‑93 (MQ=255)

                                                                          ttCGGGAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGGTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCgt               >  2:4493320/1‑92 (MQ=255)

                                                                              ggAGATAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCAt           >  2:6231803/1‑93 (MQ=255)

                                                                               gagaTAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATa          >  2:2258030/1‑93 (MQ=255)

                                                                                agaTAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATaa         >  1:456866/1‑93 (MQ=255)

                                                                                  aTAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAat       <  1:6276499/93‑1 (MQ=255)

                                                                                   tAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAata      >  1:7110734/1‑93 (MQ=255)

                                                                                   tAACTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAata      <  2:3340067/93‑1 (MQ=255)

                                                                                    aaCTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAatat     <  2:2909755/93‑1 (MQ=255)

                                                                                      cTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATAtgc   <  2:773227/93‑1 (MQ=255)

                                                                                      cTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATAtgc   <  2:3901398/93‑1 (MQ=255)

                                                                                      cTGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATAtgc   >  2:5683672/1‑93 (MQ=255)

                                                                                       tGGGACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATAtgct  <  2:1916173/93‑1 (MQ=255)

CCAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCT  >  NC_000913/1400298‑1400477

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