BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,957,451	T→C	L140L (CTA→CTG)	b1875 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,957,451	0	T	C	99.2%	388.5 / ‑4.5	121	L140L (CTA→CTG)	b1875	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (52/68); total (52/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

Only 100 of 126 total aligned reads displayed.

GTTGTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGTAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCA  >  NC_000913/1957362‑1957540

gttgttTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAgg                                                                                        <  2:2967759/93‑1 (MQ=255)

 ttgttTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGc                                                                                       <  2:64930/93‑1 (MQ=255)

 ttgttTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGc                                                                                       >  2:39195/1‑93 (MQ=255)

  tgttTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCt                                                                                      <  1:2104143/93‑1 (MQ=255)

   gttTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTg                                                                                     <  1:4684128/93‑1 (MQ=255)

       tGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTc                                                                                 >  1:1034072/1‑93 (MQ=255)

          cATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCtt                                                                              <  2:4562444/93‑1 (MQ=255)

            tCGGTGACCGCCAGCGTCGGATTCGTCAATCGCTCATACTCACCTTTTGGCGAACTGGTTTTCACTGAGGTCCCCGGCAGGCTGCGTCCTtca                                                                            >  1:3024212/1‑93 (MQ=255)

             cGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTtcat                                                                           <  2:1734212/93‑1 (MQ=255)

              ggTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTtcatc                                                                          >  1:4246701/1‑93 (MQ=255)

               gTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTtcatca                                                                         <  1:3095458/93‑1 (MQ=255)

                ggACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAg                                                                        <  1:1015916/92‑1 (MQ=255)

                   ccGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGCGTCCCTCATCaaaac                                                                     >  1:4364912/1‑92 (MQ=255)

                   ccGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGaaa                                                                     >  2:4719213/1‑93 (MQ=255)

                    cGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAg                                                                    <  1:614819/93‑1 (MQ=255)

                     gCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGc                                                                   >  1:1664453/1‑93 (MQ=255)

                     gCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGc                                                                   <  2:1748843/93‑1 (MQ=255)

                      ccAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGATCTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCa                                                                  <  2:4225515/93‑1 (MQ=255)

                      ccAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCa                                                                  <  1:3237533/93‑1 (MQ=255)

                      ccAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCa                                                                  >  1:4293623/1‑93 (MQ=255)

                         gCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATg                                                               >  1:490141/1‑93 (MQ=255)

                         gCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATg                                                               <  1:2153779/93‑1 (MQ=255)

                           gTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGTAGGCTGAGTCCTTCATCAGAAAGCAATGcc                                                             <  2:4528385/93‑1 (MQ=255)

                           gTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGcc                                                             >  1:4628388/1‑93 (MQ=255)

                           gTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGcc                                                             >  2:1301316/1‑93 (MQ=255)

                           gTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGcc                                                             >  2:1443200/1‑93 (MQ=255)

                           gTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGcc                                                             >  1:429442/1‑93 (MQ=255)

                            tCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCa                                                            >  1:79341/1‑93 (MQ=255)

                             cGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAg                                                           >  1:133143/1‑93 (MQ=255)

                             cGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAg                                                           <  1:3769985/93‑1 (MQ=255)

                              ggATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAgc                                                          <  2:2414564/93‑1 (MQ=255)

                              ggATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAgc                                                          <  1:2900169/93‑1 (MQ=255)

                               gATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAgcg                                                         <  2:4292087/93‑1 (MQ=255)

                               gATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAgcg                                                         >  2:308896/1‑93 (MQ=255)

                                 ttCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCa                                                       <  2:3833406/93‑1 (MQ=255)

                                 ttCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGGTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCa                                                       >  1:837982/1‑93 (MQ=255)

                                   cGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACg                                                     <  2:1605852/93‑1 (MQ=255)

                                   cGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACg                                                     <  2:1992544/93‑1 (MQ=255)

                                    tgTAATCGCTCATGCTCGCCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACgg                                                    <  1:4503829/92‑1 (MQ=255)

                                      tAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGc                                                  >  1:1706820/1‑93 (MQ=255)

                                       aaTCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCa                                                 >  1:3754189/1‑93 (MQ=255)

                                        aTCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAt                                                >  2:217703/1‑93 (MQ=255)

                                        aTCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAt                                                >  1:3755024/1‑93 (MQ=255)

                                        aTCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAt                                                <  2:2605299/93‑1 (MQ=255)

                                         tCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAtt                                               >  2:3518642/1‑93 (MQ=255)

                                         tCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAtt                                               <  1:4539629/93‑1 (MQ=255)

                                         tCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAtt                                               <  1:4456395/93‑1 (MQ=255)

                                          cGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTc                                              <  2:431289/93‑1 (MQ=255)

                                          cGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTc                                              <  1:1028051/93‑1 (MQ=255)

                                            cTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAg                                            >  2:674478/1‑93 (MQ=255)

                                            cTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAg                                            >  2:2876774/1‑93 (MQ=255)

                                            cTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAg                                            >  2:527187/1‑93 (MQ=255)

                                            cTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAg                                            >  2:4284682/1‑93 (MQ=255)

                                            cTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAg                                            >  2:1851113/1‑93 (MQ=255)

                                             tCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGc                                           >  1:95069/1‑93 (MQ=255)

                                             tCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGc                                           >  2:2799564/1‑93 (MQ=255)

                                             tCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGc                                           <  1:2305083/93‑1 (MQ=255)

                                             tCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGc                                           <  2:1404076/93‑1 (MQ=255)

                                             tCATGCTCACCTTATGGCGCACTGATTTTCACTGAGATCCCAGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCaaa                                           <  1:4204240/93‑3 (MQ=255)

                                             ccATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGc                                           <  2:258029/92‑1 (MQ=255)

                                              cATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCg                                          >  2:4245752/1‑93 (MQ=255)

                                               aTGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGt                                         <  1:4474688/93‑1 (MQ=255)

                                                 gCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGttt                                       >  1:3110444/1‑93 (MQ=255)

                                                  cTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTc                                      <  2:2425760/93‑1 (MQ=255)

                                                  cTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTc                                      >  1:2628597/1‑93 (MQ=255)

                                                  cTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTc                                      <  1:3029149/93‑1 (MQ=255)

                                                    cACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTg                                    >  2:1600358/1‑93 (MQ=255)

                                                    cACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTg                                    <  1:178356/93‑1 (MQ=255)

                                                    cACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTg                                    >  1:4447886/1‑93 (MQ=255)

                                                     aCCTTTTGGCGAACTGGTGTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCAGCAGAAAGCAATGCCAGCCCACGGGCCTTCAGCGTTTCTgg                                   <  1:3282186/93‑1 (MQ=255)

                                                      ccTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGa                                  <  1:4127034/93‑1 (MQ=255)

                                                       cTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGAt                                 <  2:1544020/93‑1 (MQ=255)

                                                           tGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCAcc                             <  2:618313/93‑1 (MQ=255)

                                                            ggCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCAccc                            >  2:4729562/1‑93 (MQ=255)

                                                            ggCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCAccc                            <  1:3314390/93‑1 (MQ=255)

                                                              cGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCgg                          >  2:2930549/1‑93 (MQ=255)

                                                              cGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCGTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGGTTCTGGATCACCCgg                          <  2:2655000/93‑1 (MQ=255)

                                                                aaCTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCa                        >  1:19416/1‑93 (MQ=255)

                                                                  cTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCCGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGa                      >  1:904258/1‑93 (MQ=255)

                                                                   tGGTTTTCACTGAGATCCCCCGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGaa                     >  1:2627893/1‑93 (MQ=255)

                                                                     gTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACg                   <  2:4342029/93‑1 (MQ=255)

                                                                     gTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGAGCATTCAGCGTTTCTGGATCACCCGGCAGAACg                   >  1:4574797/1‑93 (MQ=255)

                                                                      ttttCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGa                  >  1:2580459/1‑93 (MQ=255)

                                                                        ttCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGAtt                <  1:39200/93‑1 (MQ=255)

                                                                        ttCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGAtt                <  2:4644577/93‑1 (MQ=255)

                                                                         tCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGAttt               <  2:348641/93‑1 (MQ=255)

                                                                         tCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGAttt               <  2:1941981/93‑1 (MQ=255)

                                                                            cTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCaa            >  2:2347939/1‑93 (MQ=255)

                                                                             tGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGGCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAat           <  1:4492594/93‑1 (MQ=255)

                                                                              gagaTCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAata          <  1:3883157/93‑1 (MQ=255)

                                                                                gaTCCCCGGCAGGCTGGGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATg        <  2:3504463/93‑1 (MQ=255)

                                                                                gaTCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATg        <  1:2148302/93‑1 (MQ=255)

                                                                                gaTCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATg        <  2:2946104/93‑1 (MQ=255)

                                                                                gaTCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATg        >  2:310716/1‑93 (MQ=255)

                                                                                 aTCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGt       >  2:3349517/1‑93 (MQ=255)

                                                                                  tCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGtt      <  1:3359763/93‑1 (MQ=255)

                                                                                  tCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGtt      >  2:4541971/1‑93 (MQ=255)

                                                                                   ccccGGCAGGCTGAGTCCTTCATCAGAAAGCAAACCCCACGCACGGGCCTTCAGCGGTTCTGGATCAGCCCGCAGACAGATTTCAATATGTTc     >  2:1123741/1‑93 (MQ=255)

                                                                                     ccGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTccc   <  2:721553/93‑1 (MQ=255)

                                                                                      cGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCa  <  2:3186250/93‑1 (MQ=255)

GTTGTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGTAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCA  >  NC_000913/1957362‑1957540

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