BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,997,285	(T)_6→7	coding (629/756 nt)	b2865 ←	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,997,279	1	.	T	98.4%	681.7 / 2.2	189	coding (635/756 nt)	b2865	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base . (1/2); new base T (80/106); total (81/108)
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.06e-01

Only 100 of 220 total aligned reads displayed.

GCAAGTAACGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑TTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTT  >  NC_000913/2997187‑2997366

gcAAGTAACGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGa                                                                                          >  2:4031155/1‑93 (MQ=255)

gcAAGTAACGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGa                                                                                          >  2:3720848/1‑93 (MQ=255)

 cAAGTAACGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑t                                                                                        >  1:1179488/1‑93 (MQ=255)

   aGTAACGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑ttt                                                                                      >  2:2976025/1‑93 (MQ=255)

   aGTAACGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑ttt                                                                                      <  2:739727/93‑1 (MQ=255)

    gTAACGTAGCGGATCAATTGCTGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑tttt                                                                                     <  2:3599657/93‑1 (MQ=255)

    gTAACGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑tttt                                                                                     >  1:4237291/1‑93 (MQ=255)

     tAACGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑ttttt                                                                                    <  2:1436644/93‑1 (MQ=255)

       aCGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑TTTTTTg                                                                                  >  2:2096924/1‑93 (MQ=255)

        cGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑tttttttg                                                                                 <  1:2922964/93‑3 (MQ=255)

         gTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑tttttttgc                                                                                >  1:3272660/1‑90 (MQ=255)

          tAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑tttttttgcc                                                                               >  1:248926/1‑89 (MQ=255)

            gCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCa                                                                              <  1:3798951/93‑1 (MQ=255)

            gCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCa                                                                              >  1:2989619/1‑93 (MQ=255)

            gCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCa                                                                              <  1:3380460/93‑1 (MQ=255)

             cGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAg                                                                             <  1:4288956/93‑1 (MQ=255)

              ggATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGc                                                                            >  2:2231181/1‑93 (MQ=255)

               gATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGcc                                                                           >  2:3336298/1‑93 (MQ=255)

                aTCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCt                                                                          <  2:1661146/93‑1 (MQ=255)

                 tCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCGCCATAGTGGCGATTTTTTTGCCCAGCCtt                                                                         >  1:2026569/1‑93 (MQ=255)

                 tCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCtt                                                                         <  1:1000846/93‑1 (MQ=255)

                 tCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCtt                                                                         >  2:1124698/1‑93 (MQ=255)

                  cAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTc                                                                        <  2:2632464/93‑1 (MQ=255)

                   aaTTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCa                                                                       <  1:161013/93‑1 (MQ=255)

                   aaTTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCa                                                                       <  1:50923/93‑1 (MQ=255)

                    aTTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCAc                                                                      >  1:4206019/1‑93 (MQ=255)

                      tGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGc                                                                    <  2:2422179/93‑1 (MQ=255)

                      tGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGc                                                                    <  2:4195110/93‑1 (MQ=255)

                       gCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCt                                                                   >  1:2077987/1‑93 (MQ=255)

                       gCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCt                                                                   <  1:778136/93‑1 (MQ=255)

                       gCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCt                                                                   <  1:4445636/93‑1 (MQ=255)

                        ccGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCtt                                                                  <  1:2269177/93‑1 (MQ=255)

                           ttGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGc                                                               >  1:2354035/1‑93 (MQ=255)

                            tGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGcc                                                              >  2:3498208/1‑93 (MQ=255)

                              cACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCa                                                            <  2:3693579/93‑1 (MQ=255)

                              cACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCa                                                            <  2:864052/93‑1 (MQ=255)

                              cACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCa                                                            >  1:922656/1‑93 (MQ=255)

                               aCGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCat                                                           >  2:2858447/1‑93 (MQ=255)

                               aCGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCat                                                           >  2:1538889/1‑93 (MQ=255)

                               aCGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCat                                                           >  1:1563903/1‑93 (MQ=255)

                               aCGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCat                                                           >  2:4466876/1‑93 (MQ=255)

                                 ggTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCatta                                                         <  1:4125737/93‑1 (MQ=255)

                                 ggTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCatta                                                         >  2:3041391/1‑93 (MQ=255)

                                 ggTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCatta                                                         <  2:4243121/93‑1 (MQ=255)

                                  gTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATATTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCattat                                                        >  1:4504739/1‑93 (MQ=255)

                                  gTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCattat                                                        <  1:3059279/93‑1 (MQ=255)

                                  gTAACGAAGCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCattat                                                        <  2:160115/93‑1 (MQ=255)

                                   tAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCattatt                                                       >  2:4649420/1‑93 (MQ=255)

                                    aaCGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTAttt                                                      <  2:1956851/93‑1 (MQ=255)

                                     aCGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTa                                                     >  1:1287950/1‑93 (MQ=255)

                                     aCGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTa                                                     >  1:1212757/1‑93 (MQ=255)

                                     aCGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTa                                                     >  1:1600978/1‑93 (MQ=255)

                                        aaTCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCa                                                  >  2:2552049/1‑93 (MQ=255)

                                        aaTCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCa                                                  <  1:2646583/93‑1 (MQ=255)

                                        aaTCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCa                                                  <  1:1714253/93‑1 (MQ=255)

                                        aaTCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCGCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCa                                                  <  1:45832/93‑1 (MQ=255)

                                         aTCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAg                                                 >  2:1857228/1‑93 (MQ=255)

                                           cTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCa                                               >  1:4141098/1‑93 (MQ=255)

                                            tGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCAt                                              <  2:3740844/93‑1 (MQ=255)

                                             ggAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATc                                             <  1:4478481/93‑1 (MQ=255)

                                             ggAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATc                                             >  2:458750/1‑93 (MQ=255)

                                              gAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCg                                            <  1:3759177/93‑1 (MQ=255)

                                              gAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCg                                            >  2:1726139/1‑93 (MQ=255)

                                               aaaTGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCgt                                           <  1:1791391/93‑1 (MQ=255)

                                                aaTGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCgtg                                          >  2:660308/1‑93 (MQ=255)

                                                 aTGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCgtgt                                         <  1:1066597/93‑1 (MQ=255)

                                                  tGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTc                                        <  2:2396403/93‑1 (MQ=255)

                                                      ggCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTa                                    <  2:1176548/93‑1 (MQ=255)

                                                        cGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATg                                  <  2:3763094/93‑1 (MQ=255)

                                                         gAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATgg                                 >  2:3922161/1‑93 (MQ=255)

                                                           aCAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGc                               <  2:3565185/93‑1 (MQ=255)

                                                            cAGATGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCg                              >  1:4375914/1‑93 (MQ=255)

                                                               aggCCCCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTaa                           <  2:1397764/91‑1 (MQ=255)

                                                               aTGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTaa                           <  1:4036108/93‑1 (MQ=255)

                                                               aTGCCGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTaa                           <  2:2719500/93‑1 (MQ=255)

                                                                  ccGCATCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGcc                        <  2:4099158/93‑1 (MQ=255)

                                                                    gCATCCGTTCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCgt                      >  2:4555815/1‑93 (MQ=255)

                                                                      aTCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCgtaa                    <  2:3531063/93‑1 (MQ=255)

                                                                      aTCCGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCgtaa                    <  1:1959564/93‑1 (MQ=255)

                                                                        ccGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCgtaatg                  <  1:2031667/93‑1 (MQ=255)

                                                                         cGTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCgtaatgt                 >  1:2728151/1‑93 (MQ=255)

                                                                          gTGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCgtaatgta                >  1:3887578/1‑93 (MQ=255)

                                                                           tGCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCgtaatgtaa               >  1:64464/1‑93 (MQ=255)

                                                                            gCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCgtaatgtaat              <  1:244847/93‑1 (MQ=255)

                                                                            gCTCCCCATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCgtaatgtaat              <  1:3099077/93‑1 (MQ=255)

                                                                               ccccATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCtt           <  1:745283/93‑1 (MQ=255)

                                                                                cccATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTc          <  1:3139636/93‑1 (MQ=255)

                                                                                 ccATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCa         <  1:3829228/93‑1 (MQ=255)

                                                                                 ccATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCa         <  2:3121126/93‑1 (MQ=255)

                                                                                 ccATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCa         >  2:1698891/1‑93 (MQ=255)

                                                                                 ccATAGTGGCGATTTTTTTGCCAAGCCTTCACGCTTTGCCCATTATTTATTAGCATCGTGTCATTTAGTGCGTAAGTCGTAATGTAATCTTCa         >  2:2946733/1‑93 (MQ=255)

                                                                                  cATAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCAc        <  2:557610/93‑1 (MQ=255)

                                                                                    tAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACtg      >  1:2795867/1‑93 (MQ=255)

                                                                                    tAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACtg      >  2:678757/1‑93 (MQ=255)

                                                                                    tAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACtg      <  1:4362350/93‑1 (MQ=255)

                                                                                    tAGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCACATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACtg      >  1:1743769/1‑93 (MQ=255)

                                                                                     aGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACtgt     <  1:2566065/93‑1 (MQ=255)

                                                                                     aGTGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACtgt     <  2:2055555/93‑1 (MQ=255)

                                                                                       tGGCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACTgtgt   >  2:782044/1‑93 (MQ=255)

                                                                                        ggCGATTTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGtt  <  2:437852/93‑1 (MQ=255)

GCAAGTAACGTAGCGGATCAATTGCCGTTGCACGGTAACGAATCTGGAAATGCAGGCGAACAGATGCCGCATCCGTGCTCCCCATAGTGGCGA‑TTTTTTGCCCAGCCTTCACGCTTTGCCCATTATTTACCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTT  >  NC_000913/2997187‑2997366

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