BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	857,539	T→C	D249G (GAC→GGC)	ybiU ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	857,539	0	T	C	98.7%	693.9 / ‑2.2	227	D249G (GAC→GGC)	ybiU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (109/115); total (110/117)
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.62e-01

Only 100 of 244 total aligned reads displayed.

GACGTGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACG  >  NC_000913/857451‑857628

   gTGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCgt                                                                                    >  2:1605466/1‑93 (MQ=255)

    tGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCgtg                                                                                   >  2:1353328/1‑93 (MQ=255)

    tGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCgtg                                                                                   >  1:5828569/1‑93 (MQ=255)

        gcagcCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTAct                                                                               >  2:968143/1‑93 (MQ=255)

         cagcCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTActc                                                                              <  2:1587389/93‑1 (MQ=255)

          cgcCCCTGGCCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTActct                                                                             <  1:6859921/92‑1 (MQ=255)

           gcCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCtt                                                                            >  1:1338886/1‑93 (MQ=255)

           gcCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCtt                                                                            >  2:7758900/1‑93 (MQ=255)

              ccTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCac                                                                         >  1:2155038/1‑92 (MQ=255)

              ccTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCaa                                                                         <  2:2421336/93‑1 (MQ=255)

                tGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACt                                                                       >  1:4581452/1‑93 (MQ=255)

                   ccAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTcc                                                                    >  2:1045962/1‑93 (MQ=255)

                    cAGGCAGCATATCAGAGAGCGCTGTCGAGCCCTGGAATGTCCGAAACCCGGAATATTTGGTGGTGTAGGGCACCGTGTGCTCTTACACTTCCg                                                                   >  1:4622432/1‑93 (MQ=255)

                    cAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCg                                                                   <  1:2476891/93‑1 (MQ=255)

                    cAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCAAAAAACGGAAAATTGGGGGGGTTTGCCCCCCGGGT‑CTCTTTAACTTccct                                                                  >  1:1606879/1‑91 (MQ=255)

                      ggCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTa                                                                 >  2:3000941/1‑93 (MQ=255)

                      ggCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTa                                                                 <  1:4372916/93‑1 (MQ=255)

                       gcagcaTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTAc                                                                <  2:6988339/93‑1 (MQ=255)

                         agcaTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGa                                                              <  2:2694617/93‑1 (MQ=255)

                          gcaTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAt                                                             >  2:1733857/1‑93 (MQ=255)

                           caTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCACCTTCCGTACGAtg                                                            <  2:1930512/93‑1 (MQ=255)

                           caTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtg                                                            <  1:3404931/93‑1 (MQ=255)

                           caTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtg                                                            <  1:7159514/93‑1 (MQ=255)

                           caTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtg                                                            <  2:1245789/93‑1 (MQ=255)

                            atatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgt                                                           <  2:6868718/93‑1 (MQ=255)

                            atatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgt                                                           >  2:6322177/1‑93 (MQ=255)

                            atatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgt                                                           <  1:4501055/93‑1 (MQ=255)

                            atatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgt                                                           >  1:5479426/1‑93 (MQ=255)

                             tatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgtg                                                          >  1:5204436/1‑93 (MQ=255)

                             tatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgtg                                                          >  2:790231/1‑93 (MQ=255)

                             tatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgtg                                                          >  2:5164946/1‑93 (MQ=255)

                              atCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGc                                                         >  2:2936928/1‑93 (MQ=255)

                                cAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCg                                                       >  1:3525007/1‑93 (MQ=255)

                                  gcgagCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCa                                                     <  1:6147871/91‑1 (MQ=255)

                                  gagagCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCa                                                     <  1:1371217/93‑1 (MQ=255)

                                  gagagCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCa                                                     <  2:60887/93‑1 (MQ=255)

                                    gagCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATg                                                   >  2:3328445/1‑93 (MQ=255)

                                     agCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGc                                                  >  1:7270390/1‑93 (MQ=255)

                                       cgcTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCa                                                <  2:5345267/93‑1 (MQ=255)

                                        gcTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAg                                               >  1:2948840/1‑93 (MQ=255)

                                        gcTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAg                                               <  2:6160617/93‑1 (MQ=255)

                                         cTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAgg                                              >  1:6633020/1‑93 (MQ=255)

                                         cTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAgg                                              >  2:7056867/1‑93 (MQ=255)

                                          tGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAggg                                             >  2:5519233/1‑93 (MQ=255)

                                           gtacAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGa                                            <  2:2363105/90‑1 (MQ=255)

                                           gTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTCCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGa                                            <  2:3294874/93‑1 (MQ=255)

                                            tCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGAt                                           >  1:1966660/1‑93 (MQ=255)

                                               aGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCat                                        <  2:4366325/93‑1 (MQ=255)

                                                gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCata                                       <  2:3624720/93‑1 (MQ=255)

                                                gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCata                                       >  2:4872660/1‑93 (MQ=255)

                                                  ccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATAtt                                     <  2:3824388/93‑1 (MQ=255)

                                                  ccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATAtt                                     >  1:2998364/1‑93 (MQ=255)

                                                  ccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATAtt                                     <  2:4823118/93‑1 (MQ=255)

                                                  ccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATAtt                                     <  2:7503003/93‑1 (MQ=255)

                                                   cTGTAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTg                                    <  1:1229716/93‑1 (MQ=255)

                                                   cTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTg                                    >  1:932070/1‑93 (MQ=255)

                                                   cTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCATACGATGTGCCGCATGCCAGGGATCATATTg                                    >  1:7000751/1‑93 (MQ=255)

                                                    tGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTgc                                   >  1:443111/1‑93 (MQ=255)

                                                    tGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTATTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTgc                                   >  1:6979885/1‑93 (MQ=255)

                                                     ggAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTgcg                                  <  2:6189547/93‑1 (MQ=255)

                                                       aaTGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGcc                                >  2:5274096/1‑93 (MQ=255)

                                                       aaTGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGcc                                >  1:2886592/1‑93 (MQ=255)

                                                       aaTGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGcc                                >  1:3680729/1‑93 (MQ=255)

                                                        aTGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCAGCATGCCAGGGATCATATTGCGCCa                               <  2:457250/93‑1 (MQ=255)

                                                         tGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAg                              >  2:5596525/1‑93 (MQ=255)

                                                           tCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAt                            <  2:7636751/93‑1 (MQ=255)

                                                           tCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAt                            <  1:7812287/93‑1 (MQ=255)

                                                           tCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAt                            >  1:1895370/1‑93 (MQ=255)

                                                           tCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAt                            <  1:6440405/93‑1 (MQ=255)

                                                           tCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAt                            <  1:6435394/93‑1 (MQ=255)

                                                             cGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTg                          >  2:1871072/1‑93 (MQ=255)

                                                             cGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTg                          >  1:1879969/1‑93 (MQ=255)

                                                               aaaCACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGcc                        >  2:6722624/1‑93 (MQ=255)

                                                                aaCACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCa                       <  1:2868873/93‑1 (MQ=255)

                                                                aaCACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCa                       <  1:4853397/93‑1 (MQ=255)

                                                                 acacGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGTTGCATGCTAGGGATCATATTGCGCCAGATTGCCAt                      >  1:5393670/1‑93 (MQ=255)

                                                                 acacGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCAt                      >  1:2554867/1‑93 (MQ=255)

                                                                 acacGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCAt                      >  1:5131131/1‑93 (MQ=255)

                                                                  cacGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCAtt                     >  1:419777/1‑93 (MQ=255)

                                                                   acGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTa                    >  1:2589844/1‑93 (MQ=255)

                                                                   acGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTa                    >  1:1538101/1‑93 (MQ=255)

                                                                   acGGAACATTTGGTGGTGTTGCCCACCGTGTACTC‑TCAACTTTCCTACGATGTGCCGCATGCCAGGGATCATATTGGGCCAGATTGCCCTTaa                   >  2:2719032/1‑93 (MQ=255)

                                                                     ggAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTaaa                  <  1:4062071/93‑1 (MQ=255)

                                                                     ggAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTaaa                  <  2:316510/93‑1 (MQ=255)

                                                                      gAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAg                 <  1:7750193/93‑1 (MQ=255)

                                                                      gAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAg                 <  2:4152262/93‑1 (MQ=255)

                                                                       aaCATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGa                <  2:4399613/93‑1 (MQ=255)

                                                                       aaCATTTGGTGGTGTTGCCCACAGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCCGATTCCCATTAAAGa                <  2:5393459/93‑1 (MQ=255)

                                                                        aCATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGAc               <  2:7377817/93‑1 (MQ=255)

                                                                         cATTTGGTGGTGTTGTCCACCGTGTACTCTTCAGCTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACg              <  1:7608871/93‑1 (MQ=255)

                                                                          aTTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGt             >  1:7587118/1‑93 (MQ=255)

                                                                           tttGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGtt            <  1:4942062/93‑1 (MQ=255)

                                                                           tttGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCAAGGGATCATATTGCGCCAGATTGCCATTAAAGACGtt            >  1:5873008/1‑93 (MQ=255)

                                                                                tggtgTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGa       <  2:5924350/93‑1 (MQ=255)

                                                                                 ggtgTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGaa      >  2:1446237/1‑93 (MQ=255)

                                                                                  gtgtTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGaaa     <  2:4076914/93‑1 (MQ=255)

                                                                                   tgttgcCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGaaaa    >  2:572833/1‑93 (MQ=255)

                                                                                   tgttgcCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGaaaa    <  1:386352/93‑1 (MQ=255)

                                                                                   tgttgcCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGaaaa    <  1:362155/93‑1 (MQ=255)

                                                                                    gttgcCCACCGTGTACTCTTCAACTGCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAAc   <  1:5765892/93‑1 (MQ=255)

GACGTGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACG  >  NC_000913/857451‑857628

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