BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,210,781	T→C	M68V (ATG→GTG)	b2125 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,210,781	0	T	C	99.6%	700.8 / ‑5.3	222	M68V (ATG→GTG)	b2125	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (124/96); total (125/97)
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.90e-01

Only 100 of 235 total aligned reads displayed.

AATCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCATCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATCCCTT  >  NC_000913/2210691‑2210870

aaTCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCAccc                                                                                         <  2:1425393/93‑1 (MQ=255)

aaTCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCAccc                                                                                         >  2:7131076/1‑93 (MQ=255)

 aTCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCAcccc                                                                                        <  1:1251210/93‑1 (MQ=255)

  tCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACcccc                                                                                       >  2:2814212/1‑93 (MQ=255)

  tCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGGAAAGCACcccc                                                                                       >  1:1655937/1‑93 (MQ=255)

   cAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCcca                                                                                      >  1:1094102/1‑93 (MQ=255)

   cAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCcca                                                                                      >  1:1483616/1‑93 (MQ=255)

   cAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCcca                                                                                      <  2:4105373/93‑1 (MQ=255)

    aaaGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccac                                                                                     <  2:3918993/93‑1 (MQ=255)

    aaaGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccac                                                                                     >  1:3313558/1‑93 (MQ=255)

     aaGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacc                                                                                    >  2:756136/1‑93 (MQ=255)

     aaGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacc                                                                                    >  1:3660850/1‑93 (MQ=255)

     aaGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacc                                                                                    <  1:3475552/93‑1 (MQ=255)

     aaGGCATGTTCTTCAAAGGCTTAAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacc                                                                                    >  2:4981722/1‑93 (MQ=255)

       ggCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCAt                                                                                  >  2:2424355/1‑93 (MQ=255)

       ggCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCAt                                                                                  >  2:2175893/1‑93 (MQ=255)

         cATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCAttt                                                                                <  1:6617050/93‑1 (MQ=255)

           tGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTcc                                                                              <  1:5500555/93‑1 (MQ=255)

           tGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTcc                                                                              >  1:4179698/1‑93 (MQ=255)

             ttcttcAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAg                                                                            <  2:4791268/93‑1 (MQ=255)

             ttcttcAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAg                                                                            >  2:1074532/1‑93 (MQ=255)

              tcttcAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGa                                                                           <  1:5735253/93‑1 (MQ=255)

                ttcAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGAcc                                                                         >  1:1330003/1‑93 (MQ=255)

                  cAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCAc                                                                       <  2:3741654/93‑1 (MQ=255)

                   aaaGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACt                                                                      >  2:6025571/1‑93 (MQ=255)

                   aaaGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACt                                                                      >  1:257951/1‑93 (MQ=255)

                   aaaGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACt                                                                      >  1:4732910/1‑93 (MQ=255)

                   aaaGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACt                                                                      <  1:5956756/93‑1 (MQ=255)

                    aaGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTg                                                                     <  2:4646861/93‑1 (MQ=255)

                     aGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGa                                                                    >  1:2809129/1‑93 (MQ=255)

                      ggCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGAt                                                                   <  1:6801696/93‑1 (MQ=255)

                       gCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATg                                                                  <  2:7898567/93‑1 (MQ=255)

                        cTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgc                                                                 >  1:952275/1‑93 (MQ=255)

                        cTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgc                                                                 >  1:6531574/1‑93 (MQ=255)

                        cTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgc                                                                 >  2:618963/1‑93 (MQ=255)

                         tttAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgcg                                                                >  2:4439557/1‑93 (MQ=255)

                            aaTTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCgg                                                             >  1:3957946/1‑93 (MQ=255)

                              ttacGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCa                                                           <  2:1854501/90‑1 (MQ=255)

                              ttGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCGCCACCATTTCCAGACCACTGATGCGCGGct                                                           >  1:4504260/1‑92 (MQ=255)

                               tGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCAt                                                          >  2:2128854/1‑93 (MQ=255)

                                gCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATc                                                         >  1:6407104/1‑93 (MQ=255)

                                gCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATc                                                         <  1:3071818/93‑1 (MQ=255)

                                 cGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCTCCCCCACCATTTCCAGACCACTGATGCGCGGCATCt                                                        <  2:557649/93‑1 (MQ=255)

                                  gTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTg                                                       >  2:2608959/1‑93 (MQ=255)

                                  gTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTg                                                       <  2:7473584/93‑1 (MQ=255)

                                   tATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCGAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTgg                                                      >  1:3839399/1‑93 (MQ=255)

                                     ttCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCCCCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGat                                                    <  1:5128758/93‑1 (MQ=255)

                                     ttCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGat                                                    >  2:676453/1‑93 (MQ=255)

                                      tCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGata                                                   <  2:3442024/93‑1 (MQ=255)

                                       cGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGatat                                                  <  2:5836909/93‑1 (MQ=255)

                                        gTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGatata                                                 >  2:2555061/1‑92 (MQ=255)

                                        gTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATc                                                 >  2:5364313/1‑93 (MQ=255)

                                        gTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATc                                                 <  2:309175/93‑1 (MQ=255)

                                         tCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCg                                                >  1:950438/1‑93 (MQ=255)

                                         tCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCg                                                <  1:1315402/93‑1 (MQ=255)

                                          cAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCga                                               <  1:5988505/93‑1 (MQ=255)

                                              cGCAGTGATAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGaaa                                           <  2:4501068/93‑1 (MQ=255)

                                              cGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGaaa                                           >  1:3568805/1‑93 (MQ=255)

                                                cAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACa                                         <  2:7188742/93‑1 (MQ=255)

                                                 aGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAg                                        <  1:6800373/93‑1 (MQ=255)

                                                  gTGAGAAAAACAATATACGGGCGATGTTCCGGGGCAAGCACGCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGc                                       >  1:7679353/1‑93 (MQ=255)

                                                   tGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGCTGCGCGGCATCTGGATATCGAGAAACAGca                                      >  2:4060581/1‑93 (MQ=255)

                                                   tGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGca                                      <  1:1116135/93‑1 (MQ=255)

                                                   tGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGca                                      >  2:2860186/1‑93 (MQ=255)

                                                    gagaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGcac                                     <  1:5921887/93‑1 (MQ=255)

                                                    gagaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGcac                                     >  2:3468728/1‑93 (MQ=255)

                                                     agaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGcaca                                    <  2:7279483/93‑1 (MQ=255)

                                                     agaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGcaca                                    <  2:5022122/93‑1 (MQ=255)

                                                      gaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGcacac                                   >  1:6092700/1‑92 (MQ=255)

                                                      gaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACAt                                   >  1:1350110/1‑93 (MQ=255)

                                                      gaAAAACAATATACGGGCGATGTTCCGGGTAAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGcacac                                   >  1:3711601/1‑92 (MQ=255)

                                                       aaaaaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATc                                  <  2:1675722/93‑1 (MQ=255)

                                                        aaaaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATcc                                 >  2:4201174/1‑93 (MQ=255)

                                                        aaaaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATcc                                 >  1:7848630/1‑93 (MQ=255)

                                                        aaaaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATcc                                 <  1:5714853/93‑1 (MQ=255)

                                                         aaaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCg                                <  1:4363705/93‑1 (MQ=255)

                                                          aaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCgg                               >  1:2376438/1‑93 (MQ=255)

                                                           aCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCggg                              <  1:2276694/93‑1 (MQ=255)

                                                           aCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCggg                              >  1:4570294/1‑93 (MQ=255)

                                                            cAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGgc                             >  2:1555937/1‑93 (MQ=255)

                                                            cAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGgc                             >  2:7780389/1‑93 (MQ=255)

                                                             aaTATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGgcg                            >  2:6663186/1‑93 (MQ=255)

                                                                ataCGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAg                         >  1:5346835/1‑93 (MQ=255)

                                                                 taCGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGt                        <  2:6427319/93‑1 (MQ=255)

                                                                    gggCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTa                     <  2:7683987/93‑1 (MQ=255)

                                                                    gggCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTa                     <  1:235606/93‑1 (MQ=255)

                                                                     ggCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTAt                    <  2:5926430/93‑1 (MQ=255)

                                                                       cGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGc                  >  1:7207377/1‑93 (MQ=255)

                                                                        gATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCa                 >  2:4332691/1‑93 (MQ=255)

                                                                        gATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCa                 >  1:4914962/1‑93 (MQ=255)

                                                                         aTGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCAc                <  1:527735/93‑1 (MQ=255)

                                                                              ccGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCgcgc           >  1:6881046/1‑93 (MQ=255)

                                                                              ccGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCgcgc           <  2:5293231/93‑1 (MQ=255)

                                                                                gggTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCg         >  1:4818160/1‑93 (MQ=255)

                                                                                  gTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGAt       >  2:6124856/1‑93 (MQ=255)

                                                                                    cAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATcc     <  1:5664188/93‑1 (MQ=255)

                                                                                    cAAGCAACCCCACCATTTCCAGTCCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATcc     <  2:5148134/93‑1 (MQ=255)

                                                                                     aaGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATccc    >  1:1529960/1‑93 (MQ=255)

                                                                                      aGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATCCCt   >  2:1826087/1‑93 (MQ=255)

                                                                                       gCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATCCCtt  >  1:7557198/1‑93 (MQ=255)

AATCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCATCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATCCCTT  >  NC_000913/2210691‑2210870

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