BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,115,116	+C	coding (2013/2028 nt)	b2974 ←	putative endoglucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,115,115	1	.	C	99.2%	790.8 / ‑1.4	243	coding (2014/2028 nt)	b2974	putative endoglucanase
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (123/117); total (124/119)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 267 total aligned reads displayed.

GTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115023‑3115202

gTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAttt                                                                                          >  2:2029490/1‑93 (MQ=255)

 tttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑c                                                                                        <  1:4029841/93‑1 (MQ=255)

  ttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cc                                                                                       <  2:3057659/93‑2 (MQ=255)

   tGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cct                                                                                      <  1:1239443/93‑3 (MQ=255)

    gCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cctg                                                                                     >  2:6337001/1‑90 (MQ=255)

     cTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cctgc                                                                                    <  2:1316473/93‑5 (MQ=255)

        gTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGa                                                                                  <  2:2534602/93‑1 (MQ=255)

         tATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGAt                                                                                 <  2:5801145/93‑1 (MQ=255)

          aTCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTGCTGCGGATc                                                                                >  1:6267574/1‑93 (MQ=255)

           tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                                                                               >  1:1046811/1‑93 (MQ=255)

           tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                                                                               >  1:6400677/1‑93 (MQ=255)

           tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                                                                               >  2:3379250/1‑93 (MQ=255)

           tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                                                                               >  2:4938563/1‑93 (MQ=255)

             tGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgag                                                                             >  2:2661215/1‑93 (MQ=255)

             tGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgag                                                                             >  2:4276170/1‑93 (MQ=255)

              gCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgaga                                                                            <  2:1092823/93‑1 (MQ=255)

               cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGAt                                                                           <  2:3794799/93‑1 (MQ=255)

               cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGAt                                                                           <  1:5477178/93‑1 (MQ=255)

                 cGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCg                                                                         <  1:456858/93‑1 (MQ=255)

                   cAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGcc                                                                       >  2:4039312/1‑93 (MQ=255)

                        ggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc                                                                  <  1:735499/93‑1 (MQ=255)

                         gATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCt                                                                 >  1:3109356/1‑93 (MQ=255)

                          aTTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTa                                                                <  1:5698623/93‑1 (MQ=255)

                           ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAc                                                               <  2:138401/93‑1 (MQ=255)

                           ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAc                                                               <  1:2957679/93‑1 (MQ=255)

                            tGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACt                                                              >  2:4777046/1‑93 (MQ=255)

                            tGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACt                                                              <  1:3483447/93‑1 (MQ=255)

                             gCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACtt                                                             >  2:6195833/1‑93 (MQ=255)

                              cATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTa                                                            >  1:2767325/1‑93 (MQ=255)

                              cATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTa                                                            >  2:2931461/1‑93 (MQ=255)

                               aTAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTaa                                                           >  2:695323/1‑93 (MQ=255)

                                tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          <  2:2971604/93‑1 (MQ=255)

                                tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          >  1:3267143/1‑93 (MQ=255)

                                tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          >  1:4329219/1‑93 (MQ=255)

                                tAAGGATAGTTAGCCACCTGAGCCAGATATTCAAAGCAGTTAAGGATCAGCAGCGGCATTTCCTGCGCATCGAGATCGCCATAGCTACTTaaa                                                          <  1:1626450/93‑1 (MQ=255)

                                 aaGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAAt                                                         <  1:985651/93‑1 (MQ=255)

                                   ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGc                                                       >  2:5295127/1‑93 (MQ=255)

                                   ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGc                                                       <  1:3462057/93‑1 (MQ=255)

                                     aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                                     <  1:4295493/93‑1 (MQ=255)

                                     aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                                     >  1:5928427/1‑93 (MQ=255)

                                     aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                                     >  2:4892379/1‑93 (MQ=255)

                                        gTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAAc                                                  <  2:1754678/93‑1 (MQ=255)

                                        gTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAAc                                                  >  1:3908216/1‑93 (MQ=255)

                                         ttACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATTAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACa                                                 <  2:3643807/93‑1 (MQ=255)

                                         ttACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACa                                                 <  2:6393262/93‑1 (MQ=255)

                                         ttACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACa                                                 <  1:1850416/93‑1 (MQ=255)

                                           aCCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACaga                                               >  1:5095489/1‑93 (MQ=255)

                                           aCCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACaga                                               >  2:5056094/1‑93 (MQ=255)

                                            cccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACagag                                              >  1:1138350/1‑93 (MQ=255)

                                             ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAgaga                                             <  1:6146675/93‑1 (MQ=255)

                                             ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAgaga                                             <  2:4485233/93‑1 (MQ=255)

                                              cACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                                            <  2:2338115/93‑1 (MQ=255)

                                               aCCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATg                                           <  1:6484479/93‑1 (MQ=255)

                                                ccTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGc                                          <  2:6620332/93‑1 (MQ=255)

                                                 cTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGcc                                         >  2:6102086/1‑93 (MQ=255)

                                                 cTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGcc                                         <  1:100875/93‑1 (MQ=255)

                                                     gTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCa                                     >  2:1931926/1‑93 (MQ=255)

                                                         caTATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAAt                                 <  2:1198068/93‑1 (MQ=255)

                                                           tatTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc                               >  1:3375683/1‑93 (MQ=255)

                                                           tatTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc                               >  2:3335755/1‑93 (MQ=255)

                                                               cAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGAt                           <  2:4069417/93‑1 (MQ=255)

                                                                aaaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATg                          <  2:4385011/93‑1 (MQ=255)

                                                                aaaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATg                          >  2:5123290/1‑93 (MQ=255)

                                                                 aaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg                         >  1:6702735/1‑93 (MQ=255)

                                                                 aaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg                         >  1:4046577/1‑93 (MQ=255)

                                                                  aGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGa                        >  2:2223611/1‑93 (MQ=255)

                                                                  aGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGa                        <  1:1516221/93‑1 (MQ=255)

                                                                   gCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGaa                       <  2:5749805/93‑1 (MQ=255)

                                                                     cGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGt                     >  1:1938228/1‑93 (MQ=255)

                                                                     cGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGt                     <  2:215373/93‑1 (MQ=255)

                                                                     cGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGt                     <  1:5191422/93‑1 (MQ=255)

                                                                     cGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGt                     <  1:6667871/93‑1 (MQ=255)

                                                                      gTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc                    <  1:3007245/93‑1 (MQ=255)

                                                                       ttAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCg                   >  2:3835779/1‑93 (MQ=255)

                                                                         aaGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGaa                 <  1:5766504/93‑1 (MQ=255)

                                                                          aGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAAc                >  1:1109378/1‑93 (MQ=255)

                                                                           ggATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACg               >  1:6361232/1‑93 (MQ=255)

                                                                            gATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGc              <  2:2819416/93‑1 (MQ=255)

                                                                            gATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGc              <  1:4997785/93‑1 (MQ=255)

                                                                            gATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGc              <  2:1215896/93‑1 (MQ=255)

                                                                             aTCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCa             >  2:6030317/1‑93 (MQ=255)

                                                                             aTCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCa             >  2:4614967/1‑93 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            >  1:1162550/1‑93 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            >  1:5384849/1‑93 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            >  1:3292604/1‑93 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            <  2:1558452/93‑1 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            <  2:359073/93‑1 (MQ=255)

                                                                                agcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcg          >  1:2102396/1‑93 (MQ=255)

                                                                                 gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcgc         >  2:3900664/1‑93 (MQ=255)

                                                                                 gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcgc         >  1:2501697/1‑93 (MQ=255)

                                                                                 gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcgc         <  2:551575/93‑1 (MQ=255)

                                                                                 gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcgc         >  1:2232262/1‑93 (MQ=255)

                                                                                  cagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCt        <  1:242081/93‑1 (MQ=255)

                                                                                    gcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      <  1:2272694/93‑1 (MQ=255)

                                                                                    gcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      <  1:3682854/93‑1 (MQ=255)

                                                                                    gcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      <  1:4706622/93‑1 (MQ=255)

                                                                                      ggCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGttt    <  1:4424910/93‑1 (MQ=255)

                                                                                       gCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTc   >  1:1570808/1‑93 (MQ=255)

                                                                                       gCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTc   <  1:2243474/93‑1 (MQ=255)

                                                                                        cATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTcc  <  1:62083/93‑1 (MQ=255)

GTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115023‑3115202

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