BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	342,821	C→T	F238F (TTC→TTT)	yahK →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	342,821	0	C	T	99.3%	1085.3 / ‑4.4	300	F238F (TTC→TTT)	yahK	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (153/145); total (155/145)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-01

Only 100 of 311 total aligned reads displayed.

GGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCTGGTTGG  >  NC_000913/342731‑342910

ggCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGa                                                                                         >  1:3600650/1‑93 (MQ=255)

 gCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGAt                                                                                        >  2:4884799/1‑93 (MQ=255)

   aaaaCGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGAttt                                                                                      <  1:4920267/93‑1 (MQ=255)

   aaaaCGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGAttt                                                                                      <  1:4657781/93‑1 (MQ=255)

   aaaaCGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGAttt                                                                                      <  1:2069693/93‑1 (MQ=255)

    aaaCGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTc                                                                                     <  1:1881823/93‑1 (MQ=255)

    aaaCGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTc                                                                                     <  1:445687/93‑1 (MQ=255)

    aaaCGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTc                                                                                     >  2:3646568/1‑93 (MQ=255)

     aaCGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTct                                                                                    >  1:6852780/1‑92 (MQ=255)

      aCGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAATTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCAt                                                                                   >  1:6759130/1‑93 (MQ=255)

       cgcgAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCAtt                                                                                  >  1:3702177/1‑93 (MQ=255)

        gcgAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCAttt                                                                                 <  1:601043/93‑1 (MQ=255)

          gAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTg                                                                               <  1:5886888/93‑1 (MQ=255)

           aaGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGa                                                                              >  1:5502441/1‑93 (MQ=255)

            aGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGaa                                                                             <  2:5767077/93‑1 (MQ=255)

              cGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATa                                                                           <  2:3922545/93‑1 (MQ=255)

              cGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATa                                                                           <  2:1206769/93‑1 (MQ=255)

               ggCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATAc                                                                          <  1:7513792/93‑1 (MQ=255)

                 cAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAg                                                                        <  2:4063377/93‑1 (MQ=255)

                  aaaaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGt                                                                       <  2:6929333/93‑1 (MQ=255)

                  aaaaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGt                                                                       <  1:7879587/93‑1 (MQ=255)

                   aaaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTa                                                                      <  2:3561977/93‑1 (MQ=255)

                   aaaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTa                                                                      <  2:1480418/93‑1 (MQ=255)

                   aaaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTa                                                                      >  2:731626/1‑93 (MQ=255)

                   aaaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTa                                                                      >  1:3889552/1‑93 (MQ=255)

                   aaaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTa                                                                      >  2:3256485/1‑93 (MQ=255)

                   aaaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTa                                                                      >  2:3363460/1‑93 (MQ=255)

                    aaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAg                                                                     <  2:535833/93‑1 (MQ=255)

                    aaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAg                                                                     <  2:5248364/93‑1 (MQ=255)

                    aaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAg                                                                     <  2:2556698/93‑1 (MQ=255)

                      gCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCt                                                                   >  1:7225289/1‑93 (MQ=255)

                      gCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCt                                                                   >  2:1967994/1‑93 (MQ=255)

                        ccTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTgc                                                                 <  1:5331535/93‑1 (MQ=255)

                        ccTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTgc                                                                 <  2:2276630/93‑1 (MQ=255)

                        ccTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTgc                                                                 <  1:4078326/93‑1 (MQ=255)

                         cTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTgcg                                                                <  1:3029043/93‑1 (MQ=255)

                         cTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTgcg                                                                <  2:802661/93‑1 (MQ=255)

                           gggggCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGcc                                                              <  1:5236013/93‑1 (MQ=255)

                             gggCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCcac                                                            >  1:5010813/1‑93 (MQ=255)

                              ggCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCcaca                                                           >  1:1009272/1‑93 (MQ=255)

                                ccGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATa                                                         >  1:1131669/1‑93 (MQ=255)

                                  gATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAt                                                       >  1:7592355/1‑93 (MQ=255)

                                   aTGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtc                                                      >  2:7753988/1‑93 (MQ=255)

                                    tGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtct                                                     >  1:905376/1‑93 (MQ=255)

                                      aaGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcg                                                   <  2:4420665/93‑1 (MQ=255)

                                        gttgttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgac                                                 >  1:7333240/1‑93 (MQ=255)

                                         ttgttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgacg                                                >  2:7090871/1‑93 (MQ=255)

                                         ttgttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgacg                                                >  2:1871532/1‑93 (MQ=255)

                                           gttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAt                                              >  2:2025873/1‑93 (MQ=255)

                                            ttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAtt                                             >  2:5579251/1‑93 (MQ=255)

                                              aaCTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAtttt                                           >  1:695341/1‑93 (MQ=255)

                                              aaCTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAtttt                                           <  1:355555/93‑1 (MQ=255)

                                                 tCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacc                                        <  1:2639813/93‑1 (MQ=255)

                                                 tCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacc                                        >  2:31271/1‑93 (MQ=255)

                                                  cACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacaa                                       >  2:1124214/1‑91 (MQ=255)

                                                   aCGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTagcac                                      >  2:2987009/1‑93 (MQ=255)

                                                   aCGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTaccac                                      <  2:4750137/93‑1 (MQ=255)

                                                    ctcAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTaccacc                                     <  1:798685/91‑1 (MQ=255)

                                                     gCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCt                                    <  1:1680625/93‑1 (MQ=255)

                                                     gCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCt                                    >  2:1705997/1‑93 (MQ=255)

                                                     gCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCt                                    <  1:1289171/93‑1 (MQ=255)

                                                       aaTGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTg                                  >  2:130453/1‑93 (MQ=255)

                                                        aTGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGc                                 >  1:1618919/1‑93 (MQ=255)

                                                        aTGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGc                                 >  1:5652483/1‑93 (MQ=255)

                                                        aTGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGc                                 >  1:6444383/1‑93 (MQ=255)

                                                             gATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAg                            >  1:7113307/1‑93 (MQ=255)

                                                              aTGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGc                           >  2:2466810/1‑93 (MQ=255)

                                                               tGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCg                          <  1:972720/93‑1 (MQ=255)

                                                                gagaTGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGt                         <  2:6684829/93‑1 (MQ=255)

                                                                gagaTGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGt                         >  2:4670573/1‑93 (MQ=255)

                                                                  gaTGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGa                       <  1:3487580/93‑1 (MQ=255)

                                                                  gaTGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGa                       >  1:2137862/1‑93 (MQ=255)

                                                                   aTGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGAt                      <  1:4200482/93‑1 (MQ=255)

                                                                    tGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATg                     >  2:6007438/1‑93 (MQ=255)

                                                                         gcTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCAcc                <  2:5520653/93‑1 (MQ=255)

                                                                         gcTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCAcc                >  1:3388092/1‑93 (MQ=255)

                                                                         gcTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCAcc                >  1:2896357/1‑93 (MQ=255)

                                                                          cTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCa               <  2:5132298/93‑1 (MQ=255)

                                                                           tCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCAt              >  1:5962257/1‑93 (MQ=255)

                                                                           tCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCAt              >  1:3191175/1‑93 (MQ=255)

                                                                           tCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCAt              >  2:3076009/1‑93 (MQ=255)

                                                                             aTCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGa            <  1:6191058/93‑1 (MQ=255)

                                                                             aTCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGa            >  1:7894309/1‑93 (MQ=255)

                                                                              tCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGAc           >  2:157301/1‑93 (MQ=255)

                                                                              tCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGAc           >  2:2889461/1‑93 (MQ=255)

                                                                                 gAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCt        >  2:2699996/1‑93 (MQ=255)

                                                                                 gAAGAGTGTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCt        <  2:1346453/93‑1 (MQ=255)

                                                                                   agagTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtgg      >  2:3424125/1‑93 (MQ=255)

                                                                                   agagTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtgg      <  2:6587054/93‑1 (MQ=255)

                                                                                    gagTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggt     >  1:1809465/1‑93 (MQ=255)

                                                                                    gagTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggt     <  2:4196634/93‑1 (MQ=255)

                                                                                      gTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggttg   <  2:6764346/93‑1 (MQ=255)

                                                                                      gTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggttg   >  2:7560152/1‑93 (MQ=255)

                                                                                      gTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggttg   <  2:4259204/93‑1 (MQ=255)

                                                                                      gTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTCTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggttg   <  1:2869574/93‑1 (MQ=255)

                                                                                       ttttGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGCTGGCACCATGACGCtggttgg  <  1:5420834/93‑1 (MQ=255)

                                                                                       ttttGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggttgg  >  2:6356300/1‑93 (MQ=255)

                                                                                       ttttGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggttgg  >  1:271266/1‑93 (MQ=255)

                                                                                       ttttGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggttgg  >  1:3937299/1‑93 (MQ=255)

                                                                                       ttttGATTTCATTTTGAATACAGGAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggttgg  <  1:7710027/93‑1 (MQ=255)

GGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCTGGTTGG  >  NC_000913/342731‑342910

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: `‑`